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Abstract 

It is remarkable how the Frohhch polaron, one of the simplest examples of a Quan- 
tum Field Theoretical problem, as it basically consists of a single fcrmion interacting 
with a scalar Bose field of ion displacements, has resisted full analytical or numerical 
solution at all coupling since ~ 1950, when its Hamiltonian was first written. The 
field has been a testing ground for analytical, semi-analytical, and numerical tech- 
niques, such as path integrals, strong-coupling perturbation expansion, advanced 
variational, exact diagonalisation (ED), and quantum Monte Carlo (QMC) tech- 
niques. This article reviews recent developments in the field of continuum and 
discrete (lattice) Frohlich (bi)polarons starting with the basics and covering a num- 
ber of active directions of research. 
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I. INTRODUCTION 

Charge carriers in inorganic and organic matter interact with ion vibrations. The cor- 
responding electron-phonon interaction (EPI) causes phase transformations, including su- 
perconductivity, and dominates the transport properties of many metals and semiconduc- 
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tors. EPIs have been shown to be relevant in cuprate and other hig h-temperature su- 



perconductors through, 



2004 



Zhao and Morris 



mission (ARPES) flGweon et al. 



or example, isotope sub stitution experiments (IKhasanov et al. 



1995 



Zhao et al. 



number of earlier optical ( 



2004 



19971). high resolution angle resolved phot oe- 



Lanzara et al. 



Calvani et al. 



19891). neutron - scatt ering (ISendvka et al. 



( iReznik et al. 



1994 



2001 



Meevasana et al. 



Mihailovic et al. 



1995!) and more recent inelastic scatt ering 



1990 



2006 La 



Zamboni et al. 



20061). pump- probe (iGadermaier et al. 



2009 



Radovic et al. 



20081) and 



tunnelling (jShim et al. 1 . 12008| ) measuremeii t s. In colossal magnetoresistance (CMR) man- 



scopies ( ICampbell et al. 



2001 



ganites, isotope substitutions ( Zhao et al. I . Il996l ). X-ray and neutron scattering spectro- 



2OO3I ) and a number of other expe riments also s how a sig- 



nificant effect of EPI on the physical properties (for review see 



fore it has been suggested that th e long-range (IFrohlich I . Il954l ) and/or the molecular- 



Tokura 



2000)). There- 



type (e.g. (jjahn and Teller 



perconductors 



Devreese 



1996 



see flAlexandrov 



2005 



Miiller 



( lAlexandrov and Bratkovsky 



1937)) E 



1996 



s play sign if icant r ole in high-temperature su 



Alexandrov 



1998 



Bednorz and Miiller 



1988 



20001). and references th e rein), and in CMR manganites (see 



1999a 



Edwards 



2002 



Milhs et al 



I995I ) and references 



therein). Very recent experimenta. 



SrTiOa ( Ivan Mechelen et al. 



2008 



observations of the optical conductivity in the Nb-doped 
) reveal the evidence of the mid-infrared optical conduc- 
tivity band provided by the polaron mechanism like in many other oxides. The effective 
mass of the charge carriers is obtained by analyzing the Drude spectral weight. Defining 
the mass renormalization of the charge carriers as the ratio of the total electronic spectral 
weight and the Dr ude spectral weight, a tw ofold mass enhancement is obtained, which is 



attributed in Ref. ( van Mechelen et al. 



20081 ) to th e electron-phonon coup ling. The missing 



spectral weight is recovered according the sum rule (IDevreese et al. 



19771 ) in a mid-infrared 



optical conductivity band. This band results from the electron-phonon coupling interaction, 
traditionally associated with the polaronic nature of the charge carriers. The effective mass 
obtained from the optical spectral weights yields an intermediat e electron-phonon coupling 
strength, 3 < a < 4. Therefore it has been suggested in Ref. ( Ivan Mechelen et all l2008l ) 
that the charge transport in the Nb-doped SrTiOs is carried by large polarons. 

When EPI is sufficiently strong, electron Bloch states are affected even in the nor- 
mal phase. Phonons are also affected by conduction electrons. In doped insulators hare 
phonons are well defined in insulating parent compounds, but microscopic separation of 
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electrons and ph o nons is not so straightforward in metals and heavily doped insulators 



(Maximov et al. 



( iHohenberg and Kohn 



1997 ) , where the Born and Oppen 



1964 



Kohn and Sham 



leimer (1927) and density functional 
19651 ) methods are used. Here we have to 



start with the first principle Hamiltonian describing conduction electrons and ions coupled 
by the Coulomb forces. 

One cannot solve the corresponding Schrodinger equation perturbatively because the 
Coulomb interaction is strong. The ratio of the characteristic Coulomb energy to the kinetic 
energy is = mee^/(47r?T,e/3)^/^ ~ 1 for the electron density n^, = ZN = lO^^cm"^ (here 
and further we take the volume of the system a.s V = 1, unless specified otherwise, and 
h = c = = I). However, one can take advantage of the small value of the electron 
to ion mass ratio, rrie/M < 10~^. Ions are heavy and the amplitudes of their vibrations, 
(|u|) ~ a/I/Mu;/), near equilibrium positions are much smaller than the lattice constant 



(a = N ^^'^), (|u|)/a ~ (nie/MrsY^^ <^ 1. In this estimate we take the characteristic 



vibration frequency ujd of the order of the ion plasma frequency ujq = ^-kN Z'^e^ jM . 
Hence one can expand the Hamiltonian in powers of |u|. 

Any further progress requires a simplifying physical idea, which commonly is to approach 
the ground state of the many-electron system via a one-electron picture. In the framework 
of the local density approximation (LDA), where the Coulomb electron-electron interaction 
is replaced by an effective one-body potential, the Hamiltonian is written as 



' e—ph 



where 



(1) 

(2) 
(3) 



describe independent Bloch electrons and phonons, created (annihilated) by c]^^^ (c„ks) and 
by c/^^ (dc^u), respectively, ^nks = Enks — /i is the band energy spectrum with respect to the 
chemical potential /i, (k, q) are quasi-momenta of electrons and phonons, respectively, n is 
the electron band index, u is the phonon mode index, and s is the electron spin. The part 
of the electron-phonon interaction, which is linear in the phonon operators, is written as 



e~ph 



'2N 



E 



In 



nk— qs'^qi^ ~l~ H.C., 



(4) 



where 7nn'(Q> u) is the dimensionless matrix element. If we restrict the summations over q 
and k to the first Brillouin zone of the crystal, then He-ph should also include the summation 
over reciprocal lattice vectors G of umklapp scattering contributions where q is replaced 
by q + G. The terms of H^-ph which are quadratic and of higher orders in the phonon 
operators, c/q;^, are usually small. They play a role only for those phonons which are not 
coupled with electrons by the linear interaction, Eq.(j4]). 

The electron-electron correlation energy of a homogeneous electron system is often written 

as 

He-e = lj2^c{ci)p[p^, (5) 

q 

where 14 (q) is a matrix element, which is zero for q = because of electroneutrality and 
= Xlk s '^L'^k+qs is the density fluctuation operator. H should also include a random 
potential in doped semiconductor s and amorphous metals, which could affect the EPI matrix 
element ( iBelitz and Kirkpatrick I . Il994l ) . 

For the purpose of this review we mostly confine our discussions to a single band approx- 
imation with the EPI matrix element 7nn(q, k, u) = 7(q) depending only on the momentum 
transfer q. The approximation allows for qualitative and in many cases quantitative descrip- 
tions of essential polaronic effects in advanced materials. Nevertheless there are might be 
degenerate atomic orbitals in solids coupled to local molecular-type Jahn- Teller distortions, 
where one has to consider multi-band electron energy structures. 

Quantitative calculations of th e matrix element in the whole region of moment a can be 



performed from pseudopotentials (IBaryakhtar et al. 



1999 



Maximov et al. 



19971). On the 



other hand one can paramet erize EP I rathe r than to compute it from first principles in many 
physically important cases (IMahan I . Il990l ). There are three most important interactions in 
doped semiconductors, which are polar coupling to optical phonons (the Frohlich EPI), 
deformation potential coupling to acoustical phonons, and the local (Holstein) EPI with 
molecular type vibrations in complex lattices. While the matrix element is ill defined in 
metals, the bare phonons cjq^ and the electron band structure -Enks are well defined in 
doped semiconductors, which have their parent dielectric compounds. Here the effect of 
carriers on the crystal field and on the dynamic matrix is small while the carrier density 
is much less than the atomic one. Hence one can use the band structure and the crystal 
field of parent insulators to calculate the matrix element in doped semiconductors. The 
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interaction constant 7(q) has different g-dependence for different phonon branches. In the 
long wavelength limit (g ^ 7r/a), 7(q) oc g", where n = —1, and n = —1/2 for polar optical, 
molecular (cjq = uq)) and acoustic (cjq oc q) phonons, respectively. Not only q dependence 
is known but also the absolute values of 7(q) are well parameterized in this limit. For 
example in polar semiconductors |7(q)P = 4:ne^ / KUoq"^ , where k = (e^^ — Eq^)^^, and e and 
Eq are high-frequency and static dielectric constants, respectively. If the crystal lacks an 
inversion center to be piezoelectric, there is EPI with piezoelectric (acoustic) phonons with 
an anysotropic matrix elem ent, which also contribute to a polaron effect and a Coulomb-like 
attraction of two polarons (IMahan 1 . Il972l ). 

To get a better insight into physical constraints of the above approximation let us trans- 
form the Bloch states to the real space or Wannier states using the canonical linear transfor- 
mation of the electron operators, q = A^~^ Ylk^^^'^'^^s, where i = (m, s) includes both site 
m and spin quantum numbers. In this site (Wannier) representation the electron kinetic 
energy takes the following form 



He = ^ t(m - n)Sss'clcj, 



(6) 



where t(m) = ^^^^^ "^ "bare" hopping integral, j = (n, s'), and E]^ is the 

Bloch band dispersion in the rigid lattice. 

The electron-phonon interaction and the Coulomb correlations acquire simple forms in the 
Wannier representation, if their matrix elements i n the momentum repre s entat ion depend 
only on the momentum transfer q (here we follow (lAlexandrov and Mott I . Il995l ) ) , 



He-e = ^ ^c(m - n)ninj, 



where 



and 



^2N 



7(q)e^'i-' 



(7) 
(8) 

(9) 
(10) 



are the matrix elements of the electron-phonon and Coulomb interactions, respectively, in 
the Wannier representation for electrons, and hi = clct is the density operator. 



We see that taking the interaction matrix element depending only on the momentum 
transfer one neglects the terms in the electron-phonon and Coulomb interactions, which are 
proportional to the overlap integrals of the Wannier orbitals on different sites. This approx- 
imation is justified for narrow band materials, where the electron bandwidth is less than the 
characteristic magnitude of the crystal field potential. In the Wannier representation the 
Hamiltonian is 

H = y^t(m - n)6ss'clcj + ^^^gn^ [ui{q)d^ + H.c] 

+1 Yl - ^)^*^^- + Yl ^q(4^q + 1/2). (11) 

q 

One can transform it further using the site-representation also for phonons. The site repre- 
sentation of H(._ph is particularly convenient for the interaction with dispersionless modes, 
when cjq = ujq and the phonon polarization vector Gq = e are roughly g-independent. Intro- 
ducing the phonon site-operators = A^~^ e^'^'^dq one obtains in this case, 

He-ph = cuo ^ fl'(m - n)(e ■ e^^n)nms{d\, + d^), (12) 

n,m,s 

where 5'(m) is the dimensionless force acting between the electron on site m and the displace- 
ment of ion n, proportional to the Fourier transform of 7(q), and Cm^n = (m — n)/|m — n| 
is the unit vector in the direction from the electron on site m to the ion n. The real space 
representation is particularly convenient in parameterizing EPI in complex lattices. Atomic 
orbitals of an ion adiabatically follow its motion. Therefore the electron does not interact 
with the displacement of the ion, whose orbitals it occupies, that is (7(0) = 0. 



II. CONTINUUM POLARON 



If characteristic phonon frequencies are sufficiently low, the local deformation of ions. 



crystal lattice. This self-trapping phenomenon was predi cted by Landau 



more than 70 years ago . It was studied in gr eater detail (jPevreese 



Frohlich . 


1954; 


Pekar . 


1946; 


Rashba . 


1957) 



1996 



(Landau , 


1933) 


Fevnman . 


1955 



19571 ) in the effective mass approximation for the 



electron placed in a continuum polarizible (or deformable) medium, which leads to a so- 
called large or continuum polaron. Large polaron wave functions and corresponding lattice 
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distortions spread over many lattice sites. The self-trapping is never complete in the perfect 
lattice. Due to finite phonon frequencies ion polarisations can follow polaron motion if 
the motion is sufficiently slow. Hence, large polarons with a low kinetic energy propagate 
through the lattice as free electrons but with an enhanced effective mass. 

When the polaron binding energy Ep is larger than the halfbandwidth D of the elec- 
tron band, all states in the Bloch bands are "dressed" by phonons. In this strong-coupling 
regime, A = Ep/D > 1, the finite bandwidth becomes important, so the continuum ap- 
proximation cannot be applied. In this case the carriers are described as small or dis- 



ago ( 


'. ilaeles 




1963 




Holstein 




1959a 




1952; 


Yamashita and Kurosawa 


1 


1958[) 



Lang and Firsov 



1962 



Sewell 



1958 



Tyablikov 



19581 ). The first identifica tion of small po l arons in solids 



was made for non-stoichiometric uranium dioxide in Refs. (iDevreese 



1963 



Nagels et al. 



19631 ). Large and small polarons wer e discussed in a nurnber of review papers and textbooks 



for exami 


3le f 


Alexandrov and 


Uott . 


1994. 


1995; 


Add el . 


1968; 


Botteer and Brvksin 


Devreese 


19961: 


Firsov . 


I975I: 


Itoh and Stoneham . 


2001 


Mahan . 


1990; 


Mitra et al. 


Rashba . 


200^ 


5; 


Salje et al. . 


1995 


)• 



1985 



1987 : 



In many models of EPI the ground-state polaron ener gy is an analytical fu i iction 



of the coupling con s tant fo r any dimens i onality of space (IFehske and Trugman 



Gerlach and Lowen 



1991 



Hague et al. 



2006a 



Lowen 



198 



2007 



Peeters and Devreese 



19821 ). There is no abrupt (nonanalytical) phase transition of the ground state as the 
electron-phonon coupling increases. It is instead a crossover from Bloch states of band 
electrons or large polarons propagating with almost bare mass in a rigid lattice to heavily 
dressed Bloch states of small polarons propagating at low temperatures with an exponen- 
tially enhanced effective mass. The ground-state wave function of any polaron is delocalized 
for any coupling stren gth. This result holds for both finite-site models and infinite-site 
models (ILowen 1 . 119881 ). 



A. Pekar's polaron 

First let us briefiy discuss a single electron interacting with the lattice deformation in 
the continuum approximation, as studied by Pekar (1946). In his model a free electron 
interacts with the dielectric polarisable continuum, described by the static Eq and the optical 
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(high frequency) e dielectric constants. This is the case for carriers interacting with optical 
phonons in ionic crystals under the condition that the size of the self-trapped state is large 
compared with the lattice constant so that the lattice discreteness is irrelevant. 



1. Ground state 



Describing the ionic crystal as a polarisable dielectric continuum one should keep in 
mind that only the ionic part of the total polarisation contributes to the polaron state. 
The interaction of a carrier with valence electrons responsible for the optical properties is 
taken into account via the Hartree-Fock periodic potential and included in the band mass 
m. Therefore only ion displacements contribute to the self-trapping. Following Pekar we 
minimise the sum E{il)) of the electron kineti c energy and the potential eii ergy due to the 
self- induced polarisation field (here we follow (lAlexandrov and Mott I . Il995l )). 



dr 



2m 



ij{T) - P(r) ■ D(r) 



where 



D(r) = eV / dt 



(13) 



(14) 



is the electric field of the electron in the state with the wave function V'(i') and P is the 
ionic part of the lattice polarisation. Minimising E^ip) with respect to ilj*{r) at fixed P and 



Jdr 



1 one arrives at the equation of motion, 



r 1 

~--e rfr'P(r') ■ V'— - 
2m / ^ ^ r'- 



^(r) = Eotp{r), 



(15) 



where Eq is the polaron ground state energy. The ionic part of the total polarisation is 
obtained using the definition of the susceptibilities xo and P = (xo — x)D. The dielectric 
susceptibilities xo and x are related to the static and high frequency dielectric constants, 
respectively, (xo = (^o~ l)/47r£o5 and x = l)/47r£:), so that P = Then the equation 
of motion becomes 

^2 



2m 



dr' / dr' 



1 



■V, 



1 



Attk j j |r' — r"| |r' — r 

Differentiating by parts with the use of V^r~^ = — 47r5(r) one obtains 



^(r) = Eo^(r). (16) 



2m K j \r' 



J\\2 



^l){r) = Eoij{r) 



(17) 
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The solution of this nonhnear integro-differential equation can be found using a variational 
minimisation of the functional 

JW = J- f ,rlVm\' - I <i.<i.'Ml)lW, (18) 

2m J 2maB J |r — r| 

where = n/me^ is the effective Bohr radius. The simplest choice of the normalised trial 
function is ^^(r) = Ae~^/^^ with A = (irrp)"^^'^. Substituting the trial function into the 
functional yields JO') = T + ^U, where the kinetic energy is T = l/2mrp, and the potential 
energy U = —5/8maBrp. Minimising the functional with respect to yields the polaron 
radius, Vp = 160^/5, and the ground state energy Eq = T + U as Eq = — 0.146/ma^. This 
can be compared with the ground state energy of the hydrogen atom — O.Smge^, where me 
is the free electron mass. Their ratio is 0.3m/(meK^). In typical polar solids k > 4, so that 
the continuum polaron binding energy is about 0.25eV or less, if m ~ rrie- The potential 
energy in the ground state is U = —AT = 4Eq/3. 

The lowest photon energy u^in to excite the polaron into the bare electron band is 
^^min = \Eo\. The ion configuration does not change in the photoexcitation process of 
the polaron. However, a lower activation energy Wt is necessary, if the self-trapped state 
disappears together with the polarisation well due to thermal fluctuations, Wt = l-^ol — Ud, 
where Ud is the deformation energy. In ionic crystals 

Ud=l f drP(r)-D(r), (19) 



2 _ 

which for the ground state is Ud = 2|£'o|/3. Therefore the thermal activation energy is 
Wt = \Eq\/3. The ratio of four ch aracter i stic e nergies for the continuum Pekar's polaron is 
WT-.Ud-. torain : | f/| = 1 : 2 : 3 : 4 Jpekar I . llQSll 
Using Pekar's choice, 

ij{r) = A{1 + r/rp + pr^)e-'^'^^ (20) 

one obtains A = 0.12/rp''^, f] = 0.45/rp, the polaron radius Vp = LSla^ and a better 
estimate for the ground state energy, Eq = — 0.164/ma^ as compared to the result of the 
simplest exponential choice. 



2. Effective mass of Pekar 's polaron 



The lattice p ol arizat ion is responsible for the polaron mass enhancement 
( iLandau and Pekarl . Il948l ). Within the continuum approximation the evolution of 

11 



the lattice polarisation P(r, t) is described by a harmonic oscillator subjected to an external 
force ~ D/k: 

9^P(M) D(r,t) 
'^0 + (21) 

where ujq is the optical phonon frequency. If during the characteristic time of the lattice 
relaxation lJq^ the polaron moves a distance much less than the polaron radius, the 
polarisation practically follows the polaron motion. Hence for a slow motion with the velocity 
V -C uoaB the first term in Eq. (l2T]) is a small perturbation, so that 

P,M).i(D,M) -.„-«). (22, 

The total energy of the crystal with an extra electron. 



E = E{ip) + 2-KK j dr 



(23) 



is determined in such a way that it gives Eq. (l2Tl) when it is minimised with respect to P. 
We note that the first term of the lattice contribution to E is the deformation energy Ud, 
discussed in the previous section. The lattice part of the total energy depends on the polaron 
velocity and contributes to the effective mass. Replacing the static wave function iP{t) in all 
expressions for ?/;(r — \t) and neglecting a contribution to the total energy of higher orders 
than one obtains 

E = Eo + Ud + ^, (24) 

where 

m* = --r^-^ I rfrD(r) ■ V'D(r) (25) 

lITIUJqK J 

is the polaron mass. Using the equation V ■ D = — 4-7re|'?/'(r)p yields 

^* = ^ I dmv)\\ (26) 
Calculating the integral with the trial function Eq. fl2U|) one obtains 

m* ^ 0.02a^m, (27) 



where a is the dimensionless constant, defined as a = {e^ / K)y^m/2uJo. 

Concluding the discussion of the Pekar's polarons let us specify conditions of its existence. 
The polaron radius should be large compared with the lattice constant, rp ^ a to justify the 
effective mass approximation for the electron. Hence the value of the coupling constant a 
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should not be very large, a <^ {D/zluoY^'^, where D ~ z/2ma? is the bare half-bandwidth, 
and z is the lattice coordination number. On the other hand the classical approximation for 
the lattice polarisation is justified if the number of phonons taking part in the polaron cloud 
is large. This number is of the order of Ud/ojQ. The total energy of the immobile polaron 
and the deformed lattice is = — O.lOQa^cuo and Ud = 0.218a;^co'o, respectively. Then a 
is bounded from below by the condition Ud/ooo ^ 1 as ^ 5. The typical adiabatic 
ratio D/ujq is of the order of 10 to 100. In fact in many transition metal oxides with 
narrow bands and high optical phonon frequencies this ratio is about 10 or even less, which 
makes the continuum stron g-coupling polaron hard to be realised in oxides and related ionic 



compounds with light ions (lAlexandrov and Mott 



19951). 



B. Weak-coupling Frohlich polaron 



Frohhch (195 



(IFrohlich et al. 



) app lied the second quantisation form of the electron-lattice interaction 



1950l ) to describe the large polaron in the weak-coupling regime, a < 1, 



where the quantum nature of lattice polarisation becomes important, 

^ = + + h.c) + ^CUq(4rfq + 1/2). (28) 

q q 
The quantum states of the noninteracting electron and phonons are specified by the electron 

momentum k and the phonon occupation number ((iqrfq) = = 0, 1, 2, ...cxd. At zero 

temperature the unperturbed state is the vacuum, |0), of phonons and the electron plane 

wave 

|k,0) =e^^-'"|0). (29) 

While the coupling is weak one can apply the perturbation theory. The interaction couples 
the state Eq. fl29|) with the energy k"^ /2m and the state with the energy (k — q)^/2m + ujq of 
a single phonon with momentum q and the electron with momentum k — q, |k — q, Iq) = 
g«(k-q)-r|]^^\j^ The corresponding matrix element is (k — q, Iqlifg-p/ilk, 0) = V^. There is no 
diagonal matrix elements of Hf,-phi so the second order energy E]^ is 

P _ IT/ |2 

^ 2m ^ (k- q)V2m + o;o - P/2m' ^ ' 

There is no imaginary part of -Ek for a slow electron with k < qp, where 

Qp = min{muo/q + q/2) = {2muJoy^^, (31) 
13 



which means that the momentum is conserved. Evaluating the integrals one arrives at 

Ek = ; — arcsm — , (32) 



2m k \qp 
which for a very slow motion k <^ Qp yields 

^ -«^o + (33) 

Here the first term is minus the polaron binding energy, —Ep. The effective mass of the 
polaron is enhanced as 

TTl 

m* = - — ^~m(l + a/6) (34) 

due to a phonon "cloud" accompanying the slow polaron. The number of virtual phonons 
Nph in the cloud is given by taking the expectation value of the phonon number operator, 
Nph = (Xlq^q'^q)' where bra and ket refer to the perturbed state, 

l) = l°) + i: feV2rn-(k-'^q-)V2,n-c.J '^->- '^^^ 

For the polaron at rest (k = 0) one obtains Nph. = a/2. Hence the Frohlich coupling constant, 
a, measures the cloud 'thickness'. One can also calculate the lattice charge density induced 
by the electron. The electrostatic potential e0(r) is given by the average of the interaction 
term of the Hamiltonian 

e0(r) = (J]yqe*'iX + ^-C-), (36) 
q 

and the charge density p(r) is related to the electrostatic potential by Poisson's equation 
V0 = — 47rp. Using these equations one obtains 

/ X 1 ^ g^|Kipcos(q- r) 

which yields 

p(r) = -^^. (38) 

The mean extension of the phonon cloud, which can be taken as the radius of the weak 
coupling polaron, is = and the total induced charge is Q = / drp{r) = —e/n. 
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C. Lee-Low-Pines transformation 



One can put the Fr ohlich result on a variational basis by applying the Lee-Low-Pines 
(LLP) transformation (ILee et al. I . llQSSh. which removes the elec t ron coordinate, followe d 



by the displacement transformation ( 



Gurari. Il953 



Lee and Pines 



1952 



Tyablikov 



1952). 



The latter serves to account for that part of the lattice polarisation which follows the electron 
instantaneously. The remaining part of the polarisation field turns out to be small, if the 
coupling constant is not extremely large. In the opposite extreme limit, which is Pekar's 



expan sion in descending powers of a ( lAllcock 



„, ^^^^ — 

(Allcock. 


1963; 


Boeohubov . 


1950: 



Bogoliubov Jr. 



1994I ). Alternatively one can apply Feynman's path-integral formalism to remove the phonon 
f ield at the exp ense of a non-instantaneous interaction of electron with itself flll.Dt and also 



( ISchultz 



1962)). 



The transformation can be written as 

\N)=eMS)\N), 



(39) 



where S is an ant i- Hermit ian operator: = —S. In our case |A^) is a single-electron multi- 
phonon wave function. The transformed eigenstate, |A^) satisfies the Schrodinger equation, 
H\N >= E\N), with the transformed Hamiltonian 



H = exp(5)iJexp(-5). 



(40) 



If all operators are transformed according to Eg. (1401) the physical averages remain un- 
changed. LLP transformation eliminating the electron coordinate in the Hamiltonian is 
defined as 

W = ^5^(q-r)<c?q. (41) 



The transformed Hamiltonian is obtained as 

H = - J^qrfirfq)' + $^(Vqdq + H.c] 



^0$^(4rfq+l/2). 

q 



(42) 



The electron coordinate is absent in H. Hence the eigenstates |A^) are classified with the 
momentum K, which is the conserving total momentum of the system, |A^) = e^^'^^lNph), 
where \Nph) is an eigenstate of phonons. The number of virtual phonons is not small in the 
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intermediate coupling regime. Therefore one cannot apply the perturbation theory for H. 
However, one can remove the essential part of the interaction term in the Hamiltonian by 
the displacement canonical transformation, 

S = Y,fWc,-H.c., (43) 
q 

where c-number /(q) is determined by minimisation of the ground state energy. Assuming 
that the transformed ground state is the phonon vacuum e^\Nph) = |0), one obtains the 
energy as 

where 

Only the term independent of K needs to be retained in vj for a slow polaron with K <^ qp, 

Then the energy up to the second order in K is given by 

= -auo + (47) 

where the polaron mass is m* = m(l + a/6) as in Eq. flM|) . Lee, Low and Pines evaluated 
also the corrections due to off-diagonal parts of the transformed Hamiltonian and found that 
they are small. 

D. All-coupling continuum polaron 

1. Feynman theory 



Feynman devel oped a superior 



integral formalism (IFeynman 



all-coupling continuum-polaron theory using his path- 



19551 ). He got the idea to formulate the polaron problem into 
the Lagrangian form of quantum mechanics and then eliminate the field oscillators, ". . . in 
exact analogy to Q. E. D. . . . (resulting in) ... a sum over all trajectories . . . " . The resulting 
path integral (here limited to the ground-state properties) is of the form: 

(0,/5|0,0) = jvv{r)e^ (48) 
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S = exp 



-rdT + —n7: I I , , , —-ciTda 



2 2^/2 jQ Jo 



(49) 



where (5 = 1/T. Here, the Feynman units are used: h = 1, m = 1, uq = 1. Eq. (HH!) gives 
the amplitude that an electron found at a point in space at time zero will appear at the 
same point at the imaginary time (3. The interaction term in the action function S may be 
interpreted as indicating that at "time" r, the electron behaves as if it were in a potential 

-da (50) 



23/2 







which results from the electrostatic interaction of the electron with the mean charge den- 
sity of its "previous" positions, weighted with the function e"'"^^""'. This path integral (HHjl 
with (jinD has a great intuitive appeal: it shows the polaron problem as an equivalent one- 
particle problem in which the interaction, non-local in time or "retarded", occurs between 
the electron and itself. 

Subsequently Feynman introduced a variational principle for path integrals to study the 
polaron. He then simulated the interaction between the electron and the polarization modes 
by a harmonic interaction (with force constant k) between a hypothetical ("fictitious") 
particle with mass M and the electron. Within his model, the action function S (H9l) is 
imitated by a quadratic trial action (non-local in time): 



5*0 = exp 



rP f 2 n rP rP "i 



where the interaction potential (!50!) is replaced by a parabolic potential 

^ j\r{r)-v{a)fe~-\-'^\da (52) 

with the weight function e""^'"^"""'. The variational parameters C and w in Eq. flS21) are 
adjusted in order to partly compensate for the error of exploiting the trial potential (152]) 
instead of the true potential (l50l) . Following the Feynman approach, an upper bound for 
the polaron ground-state energy can be written down as 

E = E,- hm l(5-5o)o, (53) 

/3->oo p 

where S is the exact action functional of the polaron problem, while 5*0 is the trial action 
functional, which corresponds to the above model system, Eq is the ground-state energy of 
the model system, and 
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The parameters of the model system C and w are found from the variational condition that 
they provide a mimimum to the upper bound for the ground state energy of Eq. fl53|) . (For 
the details of the calculation, see subsection IIIl.JI ) At nonzero temperatures, the best values 
of the model par ameters ca n be determined from a variational princip le for the free energy 
( iFeynmaru . Il972l ). see Refs. (IKrivoglaz and Pekarl . 



1957 



Osaka 



19591). 




FIG. 1 Feynman-polaron energy as a function of a: the all-coupling theory. 

Eq. (l53|) constitutes an upper bound for the polaron self-energy at all a, which in Fig. [1] 
is compared with the results of weak- and strong-coupling expansions. The weak-coupling 
expansions of Feynman for the ground-state energy and the effective mass of the polaron 
are: 

-a-0.0123a2- 0.00064a3-... (a^O), (55) 



— = l + - + 0.025a2 + ... (a^O). (56) 
m 6 

In the strong-coupling limit Feynman found for the ground-state energy: 

-O.lOea^ - 2.83 - . . . (a^oo) (57) 



and for the polaron mass: 



3Z)l"; nnono„4 



m m 
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0.0202a* + ... (a^oo). (58) 



Becker et al. (1983), using a Monte Carlo calculation, derived the ground-state energy of 
a polaron as Eq = lim/3_»oo AF, where AF = Fp — Fp with Fp the free energy per polaron and 



FO = [3/ (2/5)] In (27r/?) thefreeener 



computation in Ref. (IBecker et al. 



g y per electron. The value Pujq = 25, used for the actual 
19831 ). corresponds to T/Tq = 0.04 (T^ = fkJi^o/kB] 



ti^Lo = is th e longitudina. 



authors of Ref. ( Becker et al. 



(LO) optical phonon frequency in conventional units). The 



19831 ) actually calculated the free energy AF, rather than the 



polaron ground-state energy. To investigate the importance of temperature effects on AF, 



(IPeeters and Devreese 



19851 ) considered the polaron energy as obtained by Osaka (1959), 



who generalized the Feynman polaron theory to nonzero temperatures: 



AF 

UJQ 



lln 



sinh 



^' sinh 



3 v^-w' 

4 V 



COth ^ - 

2 Pov 



L I V uyj Jo y/DiZ)' 

l) , and 



where [3q = Pljq, n {u 
D(u) = 



1/ (e^^ 



w"^ u 
^2 



u 
Jo 



w 



An (v) sinh' 



vu\ 



(59) 



(60) 



This result is variational, with variational parameters v and w, and gives an upper bound 
to the exact polaron free energy. — AF increases with increasing temperature and the effect 
of temperature on AF increases with increasing a. The values for the free energy obtained 
analytically from the Feynman polaron model are lower than the published MC results for 
a < 2 and a > 4 (but lie within the 1% error of the Monte Carlo results). Since the Feynman 
result for the polaron free energy is an upper bound to the exact result, we conclude that 
for a < 2 and a > 4 the results of th e Feynman model are closer to the exact result than 
the MC results of flBecker et a/.l . Il983l ). 



2. Diagrammatic Monte- Carlo algorithm 

Mishchenko et al. (2000) performed a study of the Fr ohlich polaron on the basis of th e 



Diagrammatic Quantum Monte Carlo (DQMC) method (iProkof'ev and Svistunov I . Il998l ). 
This method is based on the direct summation of Feynman diagrams for Green's functions in 
momentum space. The basic object of their investigation is the Matsubara Green's function 
of the polaron in the momentum (k)-imaginary time (r) representation 

G{k,r) = (vac|ck(r)4(0)|vac),r > 0, (61) 
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TABLE I Comparison between the free energy of the Feynman polaron theory, —{AF)p, and 



the Monte Carlo results of Ref. (Becker et al 



19M), -(AF)mc, for T/ Tn = 0.04. The relative 



differ ence is defined as A = 100 x [(AF)f — (AF)mc]/(A-F)mc- (From ([Peeters and Devreese 



19851 )) 



a 


-(AF)p 


-(AF)mc 


A(%) 


0.5 


0.50860 


0.505 


0.71 


1.0 


1.02429 


1.020 


0.42 


1.5 


1.54776 


1.545 


0.18 


2.0 


2.07979 


2.080 


-0.010 


2.5 


2.62137 


2.627 


-0.21 


3.0 


3.17365 


3.184 


-0.32 


3.5 


3.73814 


3.747 


-0.24 


4.0 


4.31670 


4.314 


0.063 



Ck(r) = exp(ifr)ckexp(— 



(62) 



In terms of a complete set {|z/)} of eigenstates of the Frohlich polaron Hamiltonian H, so 
that for a given k the relation Hluik.)) = E,y(k)\u{{k)) and if|vac) = Eo\v8ic), the expansion 
of the Green's function fl6Tl) 



G(k,r) = 5^ |(z/|4|vac)|2exp{-[E,(k) - Eo]r} 



(63) 



follows straightforwardly. For the calculations discussed in what follows, Eo = 0. 

The spectral function (Lehmann function) g\i{uj) is defined through the representation of 
the Green's function f l63|) in the form 



G(k, 



5k(fi) = 5^(5[f]-E,(k)]|(z/|4|vac) 



(64) 
(65) 



This spectral function is normalized, gk{Q)dQ = 1. It can be interpreted as the prob- 
ability that a polaron has momentum k and energy Q. The significance of the spectral 
function (165|) is determined by the fact that it has poles (sharp peaks) at frequencies, which 
correspond to stable (metastable) quasiparticle states. 
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If, for a given k, there is a stable state at energy E{h), the spectral function takes the 
form 

g^in) = zS'^S[n-Eik)]..., (66) 

where is the weight of the bare-electron state. The energy -E'gs(k) and the weight ^o^s 
for the polaron ground state can be extracted from the Green's function behaviour at long 
times: 

G(k, T > cuq') ZS"^ exp[-E(k)r]. (67) 
Similarly to Eq. fl6T]) . the A^-phonon Green's function is defined: 

GAr(k,r;qi,...,qAr) = (vac|c/q^(r)...rfq^(r)Cp(r)c{,(0)rfJj^(0)...c/J^^(0)|vac), r > 0, (68) 

N 

i-i 

From the asymptotic properties of the Green's functions fIBS]) at long times, the characteris- 
tics of the polaron ground state are found. In particular, the weight of the iV-phonon state 
for the polaron ground state is given by 

Giv(k,r > W(7^;qi,...,qjv) ^ zj^(qi, q^) exp[-E(k)r]. (69) 

A standard diagrammatic expansion of the above described Green's functions generates a 
series of Feynman diagrams. The following function is further introduced: 

oo „ 

P(k, r) = ^(k, r) + ^ / rfqi...(iqArG'Ar(k, r; qi, qjv), (70) 

N=l 

where are irreducible A^-phonon Green's functions (which do not contain disconnected 
phonon propagators). From Eqs. (!67l) . (!69l) and the completeness condition for the non- 
degenerate ground state 

oo „ 

+ E / c?qi-^^qA^4^^(qi,-,q^) = i (7i) 

N=l 

it follows that the polaron ground state energy is determined by the asymptotic behaviour 
of the function flTOjl : 

P(k, r > uj^^) ^ exp[-P(k)r]. (72) 

The function P(k, r) is an infinite series of integrals containing an ever increasing number 
of integration variables. The essence of the DQMC method is to construct a process, which 
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generates continuum random variables (k, r) with a distribution function that coincides 
exactly with P(k, r). Taking into account Eq. fITOl) and the diagrammatic rules, P(k, r) is 
identified with the distribution function 

oo „ 

Q{{y}) = / dXl...dx„^F.n,{im,{y} 

, Xi, . . . , X^j , (73) 

m=0 5„ 

where the external variables {y} include k, r, a and A^, while the internal variables de- 
scribe the topology of the diagram (labelled with ^m), times assigned to electron-phonon 
vertices and momenta of phonon propagators. The diagrammatic Monte Carlo process 
is a numeric procedure, which samples various diagrams in pa rameter space and colle cts 



19531) is 



statistics for Q{{y}) according to the Metropolis algorithm (IMetropolis et al. 
such a way that — when the process is repeated a large number of times — the result 
converges to the exact answer. The distribution function given by the convergent series 
fl73|) is simulated within the process of sequential stochastic generation of diagrams de- 
scribed by functions Fm{C,m, {y},xi, ...,Xm)- Further, using Eq. ( j64l) . the spectral function 
qkiu) is obtained applying a st ochastic optimization technique. We refer to the review 



( iMishchenko and Nagaosa 



2OO7I ) for further details on the DQMC and stochastic optimiza- 



tion, where information on the excited states of the polaron is also derived by the analytic 
continuation of the imaginary time Green's functions to real frequencies. 

DQMC confirms that for a > 1 the bare-electron state in the polaron wave function is 
no longer the dominant contribution and perturbation theory is not adequate. The bare- 
electron weight Zq for the polaron ground state, as a f unction of the po l aron coupling 



constant, rapidly vanishes for a ^ 3. It is suggested in (IMishchenko et all 120001 ) that in 



the interval 3 < a < 10 the polaron ground states smoothly transforms between weak- and 
strong-coupling limits. 

Below (see sections Illl.Ct IIII.Dt IIII.G|) earlier analytical studies and results on Frohlich 
polarons are compared with DQMC results. It would be beneficial, to have an indepen- 
dent numerical check of the DQMC results. The comparison of the DQMC results for the 
"low-energy" (Q < 0) pa rt of the spectral density ^ for the Frohlich polaron at a = 0.05 



( Mishchenko et al. 



2OOOI ) demonstrates perfect agreement with the perturbation theory re- 



^ It is worth of noticing, that the spectral density is not identical to the optical absorption coefRcient, which 
is discussed below in sections IIII.Bi IIII.CI and IIII.DI 
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suit: 



a 



go{n<Q) = -5[^ + a\ 



(74) 



(m and ujq are set equal to unity). The DQMC results for the "high-energy" (f2 > 0) 
part of the spectral density significantly differ from the perturbation-theory curve. This is 
attributed to the fact that for the Frohlich polaron the perturbation-theory expression for 

a e{VL - 1) 



go{^ > 0) 



(75) 



diverges as f2 ^ 1 and, consequently, the perturbation approach is no longer adequate. The 
main difference between the DQMC spectrum of the Frohlich polaron and the perturbation- 
theory result is the broad peak in the spectral density at Q ~ 3.5. This peak appears for 
significantly larger values of the coupling constant and its weight grows with a, see Fig. [2l 



As shown in the inset, near the threshold, = 1, the spectral density behaves as y/Q — 1. 

Next to Frohlich polarons, other polaron models have been extensively investigated using 
the Monte Carlo approach. In particular Alexandrov (1996) proposed a long-range discrete 
Frohlich interaction as a model of the interaction between a hole and the oxygen ions in 
cuprate superconductors. The essence of the model is that a charge carrier moves from site 
to site on a discrete lattice (or chain in ID) and interacts with all the ions, which reside at 
the lattice sites. Numerically rigorous polaron characteristics (ground-state energy, number 
of phonons in the polaron cloud, effective mass and isotope exponent) for such a lattice 
polaron, valid for arbitrary EPl, were o btained using a path-integral continuous-time quan- 



tum Monte Carlo techn i que ( CTQMC) (j Alexandrov and Kornilovitch 



1999 



2nn6a 



Hague et al. 



Spencer et al. I . 120051 ). This study leads to a "mobile small Frohlich" polaron 

Over the years the Feynman model for the Frohlich polaron has remained the most suc- 
cessful approac h to the prqblern . The analysis of an exactly solvab le ("symmetrical") ID- 



polaron model ( 



Devreese 



(IMishchenko et al. 



2000 



proach (jCataudella et al. 



1964 



Devreese and 



Titantah et al. 



2007 



^vrard 



1964 



19681 ). Monte Carlo schemes 



20011) and, re cently, a unifying variational ap- 



De Filippis et al. I . l2003l ) demonstrate the remarkable ac- 



wave- functions, which 



curacy of Feynman's path-integral method. Using the variationa 
combine both the Landau-Pekar and the Feynman approximations, (jCataudella et al 
found in the a — > oo limit: 

Eo = (-0.108507^2 - 2.67) cuo. 
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20071) 



(76) 



o. 
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FIG. 2 The spectral density of the FrohUch polaron for various values of the coupling constant: 
a = 0.5 (dotted line), a = 1 (solid line), and a = 2 (dashed line) with energy counted from the 



ground-state energ y of the polaron. " ' 



in the inset (from (jMishchenko et al 



le ini tial fragment of the spectral density for a = 1 is shown 



2000)) 



which is slightly lower than the variational Feynman's estimate, Eq. (I50l) . and at small 
values of a, 

Eq = -auo - 0.0123a2cJo, a -> , (77) 

i.e. the same result, as Eq. flS^ . In the latter limit the correct result for the electron 
self-energy is 

E = -auo - 0.0159a2u;o, (78) 

as found by Grosjean (1957) and confirmed by Hohler and Mullensiefen (1959), and Roseler 
(1968). 

Within the path-in t egral approach, Feynman et al . studied the mobility of polarons 



(IFeynman et al. 



1962 



Thornber and Feynman 



1970l ). Subsequently the path-integral 



approach to the polaron problem was generalized and developed to become a power- 
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ful tool to study optical a bsorption and cyclotron resonance (jPevreese et al. I . Il972a : 
Peeters and Devreese I . Il986l ). 



III. RESPONSE OF CONTINUUM POLARONS 



A. Mobility 



The mobility of larg e polarons was studied within various theoretical approaches. Frohlich 



( iFrohlich 



1937 



I954J ) pointed out the typical behavior of the large-polaron mobility 



/i oc exp(u;o/3), 



(79) 



which is characteristic for weak coupling. Within the weak-coupling regime, the mo- 
bility of the polaron was then derive d, e. g., using the Boltzmann equation in 



( iHowarth and Sondheirner 



in Ref. (ILow and Pines 
weak coupling: 



1953 



I955I ). In (jLangreth and Kadanoff I . Il964l ) it was shown that for 



Osaka 



19611) and starting from t he LLP-transformation 



[exp(^o/5) - 1] 

la 



o 



(80) 



A nonperturbative analysis was based on the Feynman polaron theory, where the mobility 
yU of the polaron using the path-integral formalism was derived by Feynman et al. (1962) 
(FHIP) as a static limit, starting fr om a frequency- dependent impedance function. Details 
of the FHIP theory are given in (jPlatzmanl . Il963l ). An approximate expression for the 



impedance function of a Frohlic h polaron at al 



strengths was obtained in Refs. (IFeynman et al. 



frequencies, teinperat ures and coupling 



1962 



Platzman 



19631 ) within the path- 



integral technique. Assuming the crystal to be isotropic, an alternating electric field E 
-EoGa; exp {iVLt) induces a current in the x-direction 



(81) 



The complex function z{Q) is called the impedance function. The electric field is considered 
sufficiently weak, so that linear-response theory can be applied. The frequency- depending 
mobility /i(fi) is defined as 

= Re[z{n)]~\ (82) 
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Taking the electric charge as unity, one arrives at j = (x), where (x) is the expectation value 
of the electron displacement in the direction of the applied field: (x) = E/iQz{Q). In terms 
of time dependent variables, 



(83) 



(x(r)) = - J iG{r - a)E{a)da, 

— oo 

where G{t) has the inverse Fourier transform 

oo 

G{n) = j G{T)exp{-inT)dT = [nz{n)]-^ 

— oo 

The expected displacement at time r is 



(84) 



(X 



T)) = TT[xUiT,a)pM'-\r, 



(85) 



Here, pa is the density matrix of the system at some time a, long before the field is turned 
on, and 



U (r, a) = T exp J Ws — XsE{s)] ds 



(86) 



is the unitary operator of the development of a state in time with the complete Hamiltonian 
H — xE, where H is the Frohlich polaron Hamiltonian, and T is the time ordering operator 



( iFeynman 



195ll ). Primed operators are ordered antichronologically: 



(87) 



G{t — a) can be represented as the second derivative 

1 / d'g 



G{T-a) 



2 \dride 



(88) 



e=ri=0 



of 



^ = Tr [U{b, a)paU'~\b, a)] , b ^ oo, a ^ -oo (89) 

with E{s) = e6{s — a) +77^(5 — r), E'{s) = e6{s — a) —ri6{s — T). The initial state is chosen at 
a definite temperature T, pa oc exp{—f3H). If the time a is sufficiently far in the past, FHIP 
assume that only the phonon subsystem was in thermal equilibrium at temperature T. The 
energy of the single electron and of the electron-phonon coupling are infinitesimal (of order 
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1/V) with respect to that of the phonons. With this choice of the initial distribution, the 
phonon coordinates can be ehminated from the expression fl89|) . and the entire expression is 
reduced to a double path integral over the electron coordinates only: 

g = j j exp{i(!>)Vr{t)Vr'{t), (90) 

where (taking m = 1) 



r^2 



%/2 



la 



23/2 




dt- j [E(t) ■r(t) -E'(t) ■r'(t)]cit (91) 

— oo 

exp(-2 \t - s\) + 2P(/5) cos(t - s) exp(+i \t - s\) + 2P{[3) cos(t - s) 



|r(t) — r(s) 



|r'(t) - r'(s) 



2 [exp(-i (t - s)) + 2P{(3) cos(t - s)] 



dtds. 



|r'(t) - r(s)| 

where -P(/3) = [^^ — l] ^ . The double path integral in ( l90ll is over paths which satisfy the 
boundary condition v{t)—v'{t ) = at times t appr oaching ±oo. The expression ( !90|) with 
f l9T]) is supposed to be exact (IFeynman et a/.l . Il962l ). Clearly to provide analytical solutions 
at all a presumably is impossible. 

Following Feynman's idea to describe the ground state energy properties of a polaron 
by introducing a para bolic "retarded" intera ction of the electron with itself (see subsection 
I1I.D|) . it is assumed in (IFeynman et a/.l . ll962l ) that the dynamical behavior of the polaron can 
be described approximately by replacing $ of Eq. fl90|) by a parabolic (retarded) expression 



dt- j [E(t) ■ r(t) - E'(t) ■ r'(t)] dt 

D — OO 

OO OO 

Y J J [e"™'*"'' +2P(/3w)cosw(t-s 



(92) 



+ [r'(t) - t'{s)Y [e+'"'l*-^l + 2P{pw) cos(t - s) 
-2 [r'(t) - r{s)f [e-™(*^'') + 2P{f]w) cos w{t - s)] ] dtds. 



The parameters C and w are to be determined so as to approximate $ as closely as possible. 
However, no variational principle is known for the mobility. At zero temperature, C and w 
are fixed at the values, which provide the best upper bound for the polaron ground state 
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energy (l53i) . At finite temperatures (T 7^ 0) the parameters C an d w are determined from the 
variational principle for the polaron free energy (IFeynmaru . Il972l ) . This way of selection of the 
model parameters C and w is based on the supposition, that "the comparison Lagrangian, 
which gives a good fit to the groun d-state energy a t zero temperature, will also give the 



dynamical behaviour of the system" (IFeynman et al. 



1962 



The analytical calculation of path integrals in (pOl) in (IFeynman et all Il962l ) was per- 
formed to the first term in an expansion of exp [i ($ — $0)]: 



g = j y exp (^$0) exp [i ($ - $0)] Vv{t)VY'{t) 



90 = J J exp{i^o)Dr{t)Vr\t), 

91 = j exp(i$o)($-$o)^^r(t)Pr'(t). 
Using dHl and dHHD, one finds fromd^ 

G{Q) ^Go{n) + Gi{Q). 
The first term in the rhs of Eq. (p6l) is 



with v"^ = w"^ + AC/w. The second term in the rhs of Eq. (!96l) is 

AC fi2 



S{u) 
D{u) 



W — 1/7^ 



00 

j {I- e^^") \mS{u)dh 


2a e™ + 2P(/3)cosM 



W I V 

—iu + 



w 



u 



1 - e™" + AP{f3v) sin^ 



(93) 
(94) 
(95) 

(96) 

(97) 

(98) 
(99) 

(100) 
(101) 



From (!84l) and (!96!l . the impedance results in the form 

nz(fi) ^ [Go(l^) + G'i(fi)]"^ . (102) 
Feynman et a/. (1962) suggested to use this expression expanded to first order in GiiO): 

VLz{n) ^ [Gq{Q)Y^ - [GQ{9)]-^Gi{n). (103) 
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Substitution o: 



(IFeynman et al 



(1971) and (p8|) into Eq. (ll03p leads to the final expression (Ref. 
19621 )) for the impedance function of the Frohlich polaron: 



(104) 



An alternative derivation of the impedance function of the Frohlich polar on, based on 



simple operator technique s, was presented in ( 



Peeters and Devreese 



1983ah . The FHIP 



result was worked out in (iPeeters and Devreesd . 119841 ) in detail to get a physical insight 
into the scattering processes incorporated in the FHIP approximation. For su fficiently low 
temperature the FHIP polaron mobility takes the form (IFeynman et alV 119621 ) 



w 



3e 



(105) 



, V ) 4mujQal3 

where v and w are variational functions of a obtained from the Feynman polaron model. 

Using the Boltzmann equation for the Feynman polaron model, Kadanoff (1963) found 
the mobility, which for low temperatures can be represented as follows: 



(106) 



. V / 2mco'o« 

The weak-coupling perturbation expans ion of the low-temperature pola ron mobility as found 
using the Green's function technique (ILangreth and Kadanoff I . Il964l ) has confirmed that 
the mobility derived from the Boltzmann equation is asymptotically e xact for weak cou pling 
(a ^ 1) and at low temperatures (T ^ ujq). As first noticed in (IKadanoff I . Il963l ). the 
mobility of Eq. (11051) differs by the factor of 3/{2j3ujQ) from that derived using the polaron 
Boltzmann equation a s given by Eq. (11061) . ^ As follows from this co mparison, the result o f 



(^Feynman et al 



19621 ) is not valid whe n T — » 0. As emphasized in (Fevnman et all Il962l ) 



1983b 



19841 ) the above dis- 



and later confirmed, in particular, in (IPeeters and Devreesd . 
crepancy can be attributed to an interchange of two limits in calculating the impedance. In 
FHIP, for weak electron-phonon coupling, one takes lim^^o hmQ,_»o, whereas limc^ohmn^o 
should be calculated. It turns out that for the asymptotically correct result the mobility 



In the asymptotic limit of weak electron-phonon coupling and low temperature , the FHIP polaron mobility 



of Eg. (| 105 |) differs by the same factor o f 3/(2/3wo) from the earlier result ( Howarth and Sondheimer 



1953 : 



Low and Pines 



1955; Osaka 



later publications ijKadanoff 
for /3 > 1. 



1963 



1961), which, as pointed out in Ref. (Fevnman et 



Langreth and Kadanoff 1 . 11964 : 



Peeters and Devreese 



1963 ) and in 



19841 ) , is correct 
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at low temperatures is mainly limited by the absorption of phonons, while in the theory 
of F HIP it is the emission of phon ons which gives the dominant contribution as T goes to 



zero ( Peeters and Devreesd . 



1983bh . 



The analysis based on the Boltzmann equation takes into account the phonon emis- 
sion processes whenever the energy of the polaron is above the emission threshold. The 
independent-collision model, which underlies the Boltzmann-equation approach, however, 
fails in the "strong coupling regime" of the Frohlich polaron, when the thermal mean free 
path becomes less than the de Bro glie wavelength: in this case, the Boltzmann equati on 



cannot be expected to be adequate (IFeynman et al. 



1962 



Hellwarth and Biaggio I . IiQQQI ). 



Experimental work on alkali halides and silver halides indicates that th e mobi 



tained from Eq. (11051) describes the experimental results accurately (see flBrown 



ityob 



1963 . 



19721 ) and references therein). Measurements of mobility as a function of temperature for 



photoexcited electrons in cubic n-type Bii^SiO^n are exp 



within the Feynman approach (IHellwarth and Biaggio 



ained well in terms of large polarons 



19991 ). The experimental findings 



on electron transport in crystalline Ti02 (rutile phase) probed by T Hz time-domai i i spec - 
troscopy were quantitatively interpreted within the Feynman model (IHendry et all l2004l ). 
One of the reasons for the agreement between theory based on Eq. (11051) and experiment 
is that in the path-integral approximation to the polaron mobility, a Maxwellian distribu- 
tion for the electron velocities is assumed, when applying the adiabatic switching on of the 
Frohlich interaction. Although such a distribution is not inherent in the Frohlich interac- 
tion, its incorporation tends to favor agreement with experiment because other mechanisms 
(interaction with acoustic phonons etc.) cause a Gaussian distribution. 



B. Optical absorption at weak coupling 



At zero temperature and in the weak-coupling limit, the optical absorption of a Frohlich 
polaron is due to the elementary polaron scattering process with the absorption of incom- 
ing photon and emission of an outgoing phonon. In the weak-coupling limit (a <^ 1) 
the polaron absorption coefficient for a many-polaron gas was first obtained by Gurevich 
e.t al. (1962). Their optical- absorption coefficient is equivalent to a particular case of 
(ITempere and Devreese I . l200ll ) with the dynamic structure factor .^ (q,^;) corresponding to 
the Hartree-Fock approximation. In (ITempere and Devreese I . l200ll ) the optical absorption 
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coefficient of a many-polaron gas was shown to be given, to order a, by 

2 1/"°° 

Re[ain)]=noe^-a^ dqq^ S (q, n - uo) , (107) 

where uq is the density of charge carriers. 

In the zero-temperature hmit, starting from the Kubo formula the optical conductivity 
of a single Frohlich polaron can be represented in the form 

+ I ^"^*(^^^*-i)E^^^^(^o|[i^.(t),i^_.(o)]|^o)^t, 

(108) 

where S = +0, Bq{t) = [V'qc/q(t) + V*^di^{t)]e''^'''^^\ and |^'o) is the ground-state wave 
function of the electron-phonon system. Within the weak coupling approximation, the 
following analytic expression for the real part of the polaron optical conductivity results 
from Eq. ffTOHj) : 

Rea (Q) = ^6 (Q) + ^^^/q^,Q (Q - uo) , (109) 

where G(x) = 1 if x > 0, and zero otherwise. The spectrum of the real part of the polaron 
optical conductivity fll09p is represented in Fig. [3l 

According to Eq. (11091) . the absorption coefficient of light with frequency > by free 
polarons for a ^ takes the form 

r(fi) = -^^|^(i^M-i)^/^ (no) 

nc 6m\r 

where c is the velocity of light in vacuum, n is th e refractive index of the medium, Up is the 



concentration of polarons. A simple derivation (iDevreese et al. \ . Il97ll ) using a canonical 
transformation method gives the absorption coefficient of free polarons, which coincides 
with the result (II 101) . The step function in (IllOp reflects the fact that at zero temperature 
the absorption of light accompanied by the emission of a phonon can occur only if the 
energy of the incident photon is larger than that of a phonon (f2 > c^o)- In the weak- 
coupling limit, according to Eq. (IllOp . the absorption spectrum consists of a "one-phonon 
line" . At nonzero temperature, the absorption of a photon can be accompanied not only by 
emission, but also by absorption of one or more phonons. Similarity between the temperature 
dependence of several features of the experimental infrared absorption spectra in high-Tc 
superconductors and the temperature dependence predicted for the optical absorption of a 
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FI G. 3 Polaron opt ical conductivity for a = 1 in the weak-coupling approximation, according 



to ( Devreese 



19721 ) A (5-like cen tral peak 



(Reprinted with permission from (jDevreesd . 



(at n 



0) is schematically shown by a vertical line. 



2006al ). ©2006, Societa Itahana di Fisica.) 



single Frohlich polaron flDevreese 



( iDevreese and Tempere 



19981). 



1972 



Devreese et al. 



1972al ) has been revealed in Ref. 



Experimentally, this one-phon on line has been observe d for free polarons in the infrared 



absorption spectra of CdO-films (IFinkenrath et al. 



19691 ). which is a weakly polar material 



with a ~ 0.74. The polaron absorption band is observed in the spectral region between 6 and 
20 /im (above the LO phonon frequency). The difference between the one-polaron theoretical 
absorption and experiment in the wavelength region where polar on absorption dominates 



the s pectrum has been explained as due to many-polaron effects (ITempere and Devreese 
200ll). 



C. Optical absorption at strong coupling 



The structure of the large polaron excitation spectrum constituted a central question at 
the earl y stages of t he de velopment of polaron theory. The exactly solvable polaron model 
of Ref. (IDevreese I . ll96J) was us ed to demonstrate the existe nce of the so-called "relaxe d 
excited states" of large polarons (IDevreese and Evrard 1 . 119641 ). In Ref. (IDevreese 1 . 1 19641 ). 
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an exactly solvable ("symmetric") ID-polaron model was introduced and a n alyse d . Th e 



further study of this model was performed in Refs. fjPevreese and Evrard 



1964 



19681). 



The model consists of an electron interacting with two oscillators possessing opposite wave 
vectors: q and — q. The parity operator, which changes dq and (i_q (and also d]^ and (ilq), 
commutes with the Hamiltonian of the system. Hence, the polaron states are classified into 
even and odd states with eigenvalues of the parity operator +1 and —1, respectively. For 
the lowest even and odd states, the phonon distribution functions are plotted in Fig. HI 
upper panel, for some values of the effective coupling constant A of this "symmetric" model. 
The value of the parameter x = [g^/mcjo]^^^ for these graphs was taken to be 1, while the 
total polaron momentum P = 0. In the weak-coupling case (A ~ 0.6) is a decaying 
function of A^. When increasing A, Wn acquires a maximum, e. g. at = 8 for the lowest 
even state at A = 5.0625. The phonon distribution function Wn has the same character for 
the lowest even and the lowest odd states at all values of the number of virtual phonons in 
the ground state (as distinct from the higher states). This led to the conclusion that the 
lowest odd state is an interna l excit ed state of the polaron. 



InRef. ( iMishchenko et 0,111200^ . the structure of the Frohlich polaron cloud was inves- 



tigated using the DQMC method. Contributions of A^-phonon states to the polaron ground 
state were calculated as a function of for a few values of the coupling constant a, see FigJH 
lower panel. As follows from FigHJ the evolution from the weak-coupling case {a = 1) into 
the strong-coupling regime {a = 17) was studied. The evolution of the shape and the scale of 



the distribution of the A ^-phonon states wi t 



within DQMC method (IMishchenko et al. 



1 incr easing a as derived for a Frohlich polaron 



20001) is in notabl e agreement with the results 



obtained w ithin the "symmetric" ID polaron model (IDevreese 



1964 



Devreese and Evrard 



1964 



19681). 



The insight gained as a result of those investigations concerning the structure of the 
excited polaron states, was subsequently used to develop a theory of the optical ab- 
sorption spectra o f pola rons. The first work was limited to the strong coupling limit 
( iKartheuser et al. Lll969l ). where the impact of the internal degrees of freedom of polarons on 



their optical properti es was revealed. T 



was treated in Ref. (IKartheuser et al. 



l e opti cal absorption of light by free Frohlich polarons 



19691 ) using the polaron states obtained within th e 



adiabatic strong-coupling approximation. It was argued in Ref. ( IKartheuser et al. 



19691), 



that for sufficiently large a (« ^ 3), the (first) relaxed excited state (RES) of a polaron is 
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FIG. 4 Upper panel: The phonon distribution functions Wn in the "symmet ric" polaron mod el 



for various values of the effective coupling constant Aatx=l,P = (from ([Devreese 



mi))- 



Lower panel: Distribution of multiphonon states i n the polaron cloud within DQMC r nethod for 



various values of a (from (jMishchenko et al. 



2000)). (Reprinted with permission from ([Devreesd . 



2006a|). ©2006, Societa Itahana di Fisica.) 
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a relatively stable state, which gives rise to a "resonance" in the polaron optical absorption 
spectrum. The following scenario of a transition, which leads to a "zero-phonon" peak in 
the absorption by a strong-coupling polaron, was then suggested. If the frequency of the 
incoming photon is equal to I^res = 0.065a^u;o, the electron jumps from the ground state 
(which, at large coupling, is well-characterised by "s" -symmetry for the electron) to an ex- 
cited state ("2p"), while the lattice pol arization in the filial stat e is adapted to the "2p" 



electronic state of the polaron. In Ref. (IKartheuser et al. I . Il969l ) considering the decay of 



the RES with emission of one real phonon, it is argued, that the "zero-phonon" peak can 
be described using the Wigner-Weisskopf formula which is valid when the linewidth of that 
peak is much smaller than ujq. 

For photon energies larger than I^res + ^lo? a transition of the polaron towards the first 
scattering state, belonging to the RES, becomes possible. The final state of the optical 
absorption process then consists of a polaron in its lowest RES plus a free phonon. A "one- 
phonon sideband" then appears in the polaron absorption spectrum. This process is called 
one-phonon sideband absorption. The one-, two-, ... K-, ... phonon sidebands of the zero- 
phonon peak give rise to a broad structure in the absorption spectrum. It turns out that the 
first moment of the phonon sidebands corresponds to the Franck-Condon (FC) frequency 



To summarise, following (IKartheuser et al 



1969I ). the polaron optical absorption spec- 



trum at strong coupling is characterised by the following features (at T = 0): 

a) An absorption peak ("zero-phonon line") appears, which corresponds to a transition 
from the ground state to the first RES at I^res- 

b) For Vt > Q-RES + ujq, a phonon sideband structure arises. This sideband structure 
peaks around Qpc- Even when the zero-phonon line becomes weak, and most oscillator 
strength is in the LO-phonon sidebands, the zero-phonon line continues to determine 
the onset of the phonon sideband structure. 



The basic qualitative strong coupling behaviour predicted in Ref. (IKartheuser et al. 



19691 ). namely, zero-phonon (RES) line with a broader sideband at the high-frequency side, 



was confirmed by later studies, as discussed below. 
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D. Optical absorption at arbitrary coupling 



The optical absorpti on of the Frohhch polaron was calculated in 1972 (jPevreese 



Devreese et al. 



1972 



1972a[ ) ("DSG") for the Feynman polaron model using path integrals. Un 



til recently DSG combined with (IKartheuser et al. 



196^ constituted the ba sic picture for 



the optical absorption of the Frohlich polaron. (iPeeters and Devreese I . Il983al ) rederived the 
DSG-result using the memory function formalism (MFF). The DSG-approach is successful 
at weak electron-phonon coupling and is able to identify the excitations at intermediate 
electron-phonon coupling (3 < a < 6). In the strong coupling limit DSG still gives an 
accurate first moment for the p olaron optical abso r ption but d oes not reproduce the b road 



phonon sideband structure (cf. (IKartheuser et al. 



19691 ) and fiGoovaerts et al. 



mm). A 



comparison of the DSG results with th e optical conductivity spe ctra given by recently devel- 



oped "approximation-free " num erical (IMishchenko et al. 



(De 


i^ilippis et al. 


2006 
2007 


), see also the 
)• 



2003 



20031 ) and app roximate analytical 



20061 ) a pproaches was carried o ut re cently ( 



De Filippis et al. 



2006), see also the review articles (ICataudella et all . 120071 ) and (IMishchenko and Nagaosa 



The polaron absorpti on coefficient ^ifi) of lig ht with frequency Q at arbitrary coupling 
was first derived in Ref . (IDevreese et al. I . Il972al ) . It was represented in the form 

ImS(fi) 



r(fi) 



Air 



;iii) 



ncm[n- ReS(fi)]^ + [ImS(l])]^ 

This general expression was the starting point for a derivation of the theoretical optical 
abs orption spectrura of a s ingle Frohlich polaron at all electron-phonon coupling strengths 
by (IDevreese et al. I . Il972al l. E(f2) is the so-called "memory function", which contains the 
dynamics of the polaron and depends on a, temperature and applied external fields. 
The key contribution of (IDevreese et al. \ . Il972al ) was to introduce in the form ( lllip 
and to calculate ReS(f2), which is essentially a (technically not trivial) Kramers-Kronig 



transform of the more sim ple function ImS 



for the Feynman polaron (IFeynman et al. 



Q). Only the function ImS(f2) had been derived 



19621 ) to study the polaron mobility fj, using the 



lim (ImS(fi)/fi). 



impedance function (152]) : fi 

The basic nature of the Frohlich polaron excitations was clearly reve aled through this 



polaro n optical absorption given by Eq. ( II 111) . It was demonstrated (IDevreese et al. 



1972al ) that the FC states for Frohlich polarons are nothing else but a superposition of 
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phonon sidebands and a relatively large value of the electron-phonon coupling strength 
{a > 5.9) is needed to stabilise the relaxed excited state of the polaron. It was, further, 
revealed that at weaker coupl ing only "scatter i ng states" o f the polaron play a significant 



role in the optical absorption (jPevreese et al 



1971 



1972ah . 



In the weak coupling limit, the absorption spectrum fill II) of the polaron is determined 
by the absorption of radiation energy, which is re-emitted in the form of LO phonons. As 
a increases between approximately 3 and 6, a resonance with increasing stabilit y appears 



in the optical absorption of the Frohlich polaron of Ref. (jPevreese et al. I . Il972al ) (see Fig. 



[5]). The RES peak in the optical absorption spectrum also has a phonon sideband-structure, 
whose average transition frequency can be related to a FC-type transition. Furthermore, 
at zero temperature, the optical absorption spectrum of one polaron exhibits also a zero- 
frequency "central peak" [oc For nonzero temperature, this "central peak" smears 
out and gives rise to an "anomalous" Drude-type low-frequency component of the optical 
absorption spectrum. 

For a > 6 .5 the polaron optical ab sorption gradually develops the structure qualitatively 
proposed in (IKartheuser et al. 1 . 119691 ): a broad LO-phonon sideband structure appears with 
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FIG. 5 Optical absorption sp ectrum of a Frohlich polaron for a = 4.5, a = 5.25 and a = 6 after 



Ref. (Devreese et al. 



1972al ) (DSG). The RES peak is very intense compared with the FC peak. 
The J-h ke central peaks (at fl = 0) are schematically shown by vertical lines. The DQMC results 



of Ref. ( Mishchenko et al. 



20031 ) are shown with open circles. 
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the zero-phonon ("RES") transition as onset. Ref. (IDevreese et al. I . ll972al ) does not predict 
the hroad LO-phonon sidebands at large couphng constant, although it still gives an accurate 
first Stieltjes moment of the optical ab sorption spectrum. 



Based on (IDevreese et al. I . Il972al ). it was argued that it is rather Holstein polarons 
that determine the optical proper ties of the charge carriers in oxides like SrTiOa, BaTiOa 



Muybrechts and Devreese 



19751). while Frohlich weak-coupling polarons could be identi- 



fied e. g. in CdO (IDevreese et al. I . Il97ll ). The Frohlich coupling constants of polar semi- 



conductors and ion ic crystals are generally too small to al 



( lEagles et all Il995l ) the RES-peaks of Ref. (IDevreese et al. 



ow for a static "RES". In Ref. 



1972al ) were invoked to explain 



the optical absorption spectrum of Pr2Ni04.22- The RES-like resonances in r(r2), Eq. (Illll) . 
due to the zero's of ^2 — ReS(r2), can effectively be displaced to smaller polaron coupling 
by applying an external magnetic field 5, in which case the contribution for what is for- 
mally a "RES-type resonance" arises at — — ReS(f2) = {ujc = eB/m is the cyclotron 



clearly observed and analysed in many sohds and structures, se e e. g. 



Hodby et al. 



1987 



Miura et al. 



1997 



Peeters and Devreese 



1986 



Devreese 




1996. 


2003; 


Shi et al. . 


199? 


) and 



references therein. 

Evidence for the polaron character of charge carriers in AgBr and AgCl was obtained 
through high-precision cyclotron resonance experiments in external magnetic fields up to 



16 T. The all-coupling magneto-absorption calculated in (jPeeters and Devreese 



19861 ) leads 



to the best quantitative agreement between theory and experiment for AgBr and AgCl. 
his quantitative in terpretation of t he cyclotron resonance exper iment in AgBr and AgCl 



( Hodby et al. 



19871 ) by the theory of (jPeeters and Devreese I . Il986l ) provided one of the most 



convincing and clearest demonstrations of Frohlich polaron features in solids. The energy 
spectra of polaronic systems such as shallow donors ("bound polarons"), e. g., the Dq and 
D^ centre s, constitut e the m ost complete and detailed polaron spectroscopy realised in the 
literature (jShi et al. I . Il993l ) . 

The numerical calculatio ns of the optical conductiv ity for the Frohlich polaron performed 
within t he DQMC method (iMishchenko et al. 1 . 120031 ) confirm the analytical results derived 



in Ref. ( Devreese et al. 



1972al ) for a ^ 3. In the intermediate coupling regime 3 < a < 6, 



the low-ener gy behaviour and the p osition of the RES-peak in the opt ical conductivity 



spectrum of (IMishchenko et al. I . l2003l ) follow closely the prediction of Ref. (IDevreese et al. 
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1972al ). There are some minor quantitative differences between the two approaches in the 



intermediate couphng regime: the dominant ( "RES" ) peak is less intense in the Monte-Carlo 
numerical simulations and the second ( "FC" ) peak develops l ess prominently. The follo wing 



qualitative differences exist between the two approaches: in flMishchenko et al. 



20031 ) ■ the 



dominant peak broadens for a > 6 and the second peak does not develop, but gives rise 
to a flat shoulder in the optical conductivity spectrum at a ~ 6. As a increases beyond 
a ~ 6, the DSG results for the OC do not produce the hroad phonon sideband spectru m of 
the RES-transition that was qualitatively predicted in Ref. (IKartheuser et al. 1 . 119691 ) and 
obtained with DQMC. 

An instructive comparison between the positions of the main peak in the optical absorp- 



tion spectra of Froh 



ich polarons obtained within the DSG and DQMC approaches has been 



performed recently (IDevreese and Klimin I . I2006I ). The main-peak positions, obtained within 



DSG, have been found in good agreement with the results of DQMC for all considered values 
of a. At large a the positions of the main peak in the DSG spectra are remarkably close 
to those given by DQMC. The difference between the DSG and DQMC results is relatively 
larger at a = 8 and for a = 9.5, but even for those values of the coupling constant the 
agreement is rather good. 

It is suggested that the RES-peak at a ~ 6 in the DSG-tr eatment, as a increases , 
gradually transforms onto a FC-peak. As stated above and in Ref. ( iDevreese et al. I . ll972al ). 
DSG predicts a much too narrow FC-peak in the strong coupling limit, but still at the 
"correct" frequency. The DSG spectrum also satisfies the zero and first moment sum rules 
at all a as discussed below. 



E. Main-peak line and strong-coupling expansion 



In order to describe the OC main pe ak line width at iri t ermed iate electron-phonon cou 
pling, the DSG approach was modified (IDe Filippis et al. I . l2006l ) to incl ude additional dis- 



sipat ion processes, the strength of which is fixed by an exact sum rule (jCataudella et al. 
200^ . 



Within the memory function formalism (MFF) (IGotze and P. Wolfie 



1972 



Mori 



19651 ) 



the interaction of the charge carriers with the free phonon oscillations is expressed in terms 
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of the electron density-density correlation function, 



X(q, t) = -iQ{t) (exp [iq ■ r{t)] exp [-iq ■ r(0)]) 



;il2) 



which is evaluated in a direct way (jPeeters and Devreese 



1983a] ) using the Feynman po- 



laron model, where the electron is coupled via a harmonic force to a fictitious particle that 
simulates the phonon degrees of freedom. Within this procedure the electron density-density 
correlation function takes the form: 



X™(q,t) = -tQ{t)exp [-iqH/2M] exp [-q^R{l - e-™*)/2M] 



(113) 



where R = (M — l)/v, M (the total mass of electron a nd fictitious par ticle), and v are 
determined variationally within the path integral approach (IFeynmaru . Il955l ) . The associated 
spectral function Am(q, uj) = — 2ImXm(q, u) is a series of 6 functions centered at /2M + nv 
{n is integer). Here q^ jIM represents the energy of the center of mass of electron and 
fictitious particle, and v is the energy gap between the levels of the relative motion. To 
include dissipation (jPe Filippis et al. 1 . 12006| ). a finite lifetime was introduced for the states 
of the relative motion, which can be considered as the result of the residual EPI not included 
into the Feynman variational model. To this end, in Xm{^.,t) the factor exp [— u^t] was 
replaced with {l+it/r)''"'^ which leads to the replacement of 6 functions by Gamma functions 
with mean value and variance given respectively by q^ /2M + nv and nv/r. The parameter of 
dissipation r is determined by the third sum rule for A(q, cu), which is additional to the first 
two sum rules that are already satisfied in the DSG model without damping. As expected, 
r turns out t o be of the order of OJ^^ . If broadening of the qscillato r levels is neglected, the 



DSG results ( jPevreese et al. 



1972a 



Peeters and Devreese 



1983al ) are recovered. 



Starting from the Kubo formula, the strong-coupling polaron optical conductivi t y can 



be evaluated usin g the strong-coupl ing expansion (S C E), Refs. ( 



Cataudella et al. 



Devreese and Klimin 



2007 



20061) SCE 



De Filippis et al. 1 . 120061 ). In Refs. flDevreesd . l2007al Jbl: 

has been extended. In order to apply the extended SCE for the polaron OC, a scaling trans- 
formatio n of the coordinates and moments of the electron-phonon system is made following 



AUcock f AUcock 



19631 ). Y = a ^x, p = —iad/d^, and q = aq. This transformation allows 



us to see explicitly the order of magnitude of different terms in the Hamiltonian. Expressed 
in terms of the new variables, the Hamiltonian can be written as a sum of two terms, which 
are of different orders in powers of a, H = Hi + H2, where Hi ~ is the leading term, and 
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Ho 



a 



is the kinetic energy of the vibrating ions. The next step is the Born-Oppenheimer 



approximation (jAllcockl . Il963l ). which neglects the non-adiabatic transitions between difFer- 



ent polaron levels in calculating the dipo 



(jDevreesd . 



2007b 



Devreese and Klimin 



e-dip ole correlation function of the Kubo formula 



20061). 



Figure E] shows the polaron OC spectra for different values of a calculated numerically 
using the extended SCE with different approximations. The OC spectra calculated within 
the extended SCE approach, taking into account both the Jahn- Teller effect and the cor- 
rections of order are shown by the solid curves. The OC obtained with the leading- 
term strong-coupling approximation taking into account the Jahn- Teller effect and with the 
leading term of the Landau- Pekar adiabatic approximation are plo tted as dashed and dash- 



dotted curves, respectivel y. The full circles show the DQMC data (jPe Filippis et al. 



Mishchenko et al. 



2006 : 



20031). 



The polaron OC band of Fig. [6] obtained within the extended S CE generalises the 



Gaus sian-like polaron OC band (as given e.g. by Eq. (3) of Ref. ljDe Filippis et al. 



20061)') thanks to (i) the use of the numerically exact strong-coupling polaron wave functions 



(iMiyakd . Il975l ) and (ii) the incorporation of both static and dynamic Jahn- Teller effects. 



The polaron OC broad structure calculated within the extended SCE consists of a series 
of LO-phonon sidebands and provides a realisation — with all LO-phonons in volved for a 



19691 ). 



given a — of the strong-coupling scheme proposed in (iKartheuser et al. 

As seen from Fig. [6l the polaron OC spectra calculated within the asymptotically exact 
strong-coupling approach are shifted towards lower frequencies as compared with the OC 
spectra calculated within the LP approximation. This shift is due to the use of the numer- 
ically exact (in the strong-coupling limit) energy levels and wave functions of the internal 
excited polaron states, as well as the numerically exact self-consistent adiabatic polaron po- 
tential. Furthermore, the inclusion of the corrections of order a° leads to a shift of the OC 
spectra to lower frequencies with respect to the OC spectra calculated within the leading- 
term approximation. The value of this shift Af2„^o/^LO ~ —1.8 obtained in the present 



calculation, is c lose to the LP value AVL^^P /ujq = 41n2 — 1 1.7726 (cf. f Allcock 



1956; 



Feynmaru . 



19721 )). The distinction between the OC spectra calculated with and without the 



Jahn- Teller effect is very small. 

Starting from a ^ 9 towards larger values of a, the agreement between the extended SCE 
polaron OC spectra and the numerical DQMC data becomes gradually better, consistent 
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FIG. 6 The polaron OC calculated within the extended SCE taking into account corrections of 
order (solid curve), the OC calculated within the leading-term strong-coupling approximation 
(dashed curve), with the leading term of the Landau-Pekar (LP) adiabatic approximation (dash- 
dotted curve), and the nume rical DQMC data (open circles) for a = 7, 9, 13 and 15. (From Ref. 



(Devreese and Klimin 



20061 1. 



with the fact that the extended SCE for the polaron OC is asymptotically correct in the 



strong-coupling limit. The results of the extended SCE are qua^ 



the interpretation advanced in Ref. (IKartheuser et al. I . Il969l ). In (IKartheuser et al. 



itatively consistent with 



only the 1-LO-phonon sideband was taken into account, while in Ref. (IGoovaerts et al. 



196^) 



I973I ) 2-LO-pho non emission was i nclude d. The extended SCE carries on the program 



started in Ref. (IKartheuser et al. 



19691 ). The spectra in Fig. [HI in the strong coupling 



approximation, consist of LO-phonon sidebands to the RES (which itself has negligible 
oscillator strength in this limit, similar to the optical absorption for some colour centres in 
alkali halides). These LO-phonon sidebands form a broad FC-structure. 
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F. Comparison between optical conductivity spectra obtained by different methods 



A comparison between the optical conductivity spectra obtained with the DQMC method, 
exte nded MFF, SCE and DS G for different values of a is shown in Figs. [3, [HI taken from 
Ref. (jPe Filippis et al. 1 . 12006| ). The key results of the comparison are the following. 

First, as expected, in the weak-coupling r egime, both the extended MFF with phonon 



broad ening and DSG ( 



Devreese et al. 



data ( Mishchenko et al. 



1972al ) are in very good agreement with the DQMC 



20031 ). showing significant improvement with re s pect to the weak- 



coupling perturbation approach (IDevreese et al. 



19711 : 



Gurevich et al. 



19621 ) which pro- 



vides a good description of the OC spectra only for very small values of a. For 3 < a < 6, 
DSG predicts the essential structure of the optical absorption, with a RES-transition grad- 
ually building up for increasing a, but underestimates the peak width. The damping, intro- 
duced in the extended MFF approach, becomes crucial in this coupling regime. 

Second, comparing the peak and shoulder energies, obtained by DQMC, with the peak 



energies, given by 



VIFF , and the FC transition energies from the SCE, it is concluded 



(IDe Filippis et al. I . I2006I ) that as a increases from 6 to 10 the spectral weights rapidly switch 
from the dynamic regime, where the lattice follows the electron motion, to the adiabatic 
regime dominated by FC transitions. In the intermediate electron-phonon coupling regime, 
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tended MFF (solid line) and DSG (jPevreese et al. 



c ulated with the DQMC method (cir 



1972a 



Peeters and Devreese 



cles), ex- 



1983a|) (dotted 



line), for four different values of a. The arrow indicates the lowe r -frequ ency feature in the DQMC 



data. (Reprinted with permission from Ref. (|De Filippis et al. 
Physical Society.) 



20061 ^. ©2006 by the American 
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6 < a < 10, both adiabatic FC and non-adiabatic dynamical excitations coexist. For still 
larger coupling a > 10, the absorption spectrum consists of a broad FC-structure, built of 
LO-phonon sidebands. 

In summary, the accura te numerical results obtained from DQMC and from the recent 



analytical approximations (jPe Filippis et al. 



2003 



20061 ) confirm the ess ence of the mech- 



anism for the optical ab sorption of Frohl ich polarons, proposed in Refs. (jPevreese 



Devreese et al. 



1972 



1972a] ) combined with (IKartheuser et al. I . Il969l ) and do add important 



new insights. 



G. Sum rules for the optical conductivity spectra of Frohlich polarons 



In this section several sum rules for th e optical conduct i vity sp ectra of Frohlich polarons 
are applied to test the D SG approach (IDevreese et al. I . Il972al ) and the DQMC results 
(IMishchenko et al. 1 . 120031). The values o f the p olaron effective mass for the DQMC approach 
are taken from Ref. (IMishchenko et a/.l . l2000l ). In Tables HTl and UTTl we show the polaron 
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FIG. 8 a), b) and c) Comparison of the optical conductivity calculated with the DQMC method 
(circles), the extended MFF(solid line) and SCE (dashed line) for three different values of a. d) 
The energy of the lower- and higher-frequency features obtained by DQMC (circles and triangles, 
respectively) compared (i) with the FC transition energy calculated from the SCE (dashed line) 
and (ii) with the energy of the peak obtained from the extended MFF (solid line). In the inset the 
weights of Franck-Condon and adiabatically connected transitions are shown as a function of a . 



We have used for r] the value 1.3. (Reprinted with permission from Ref. (jDe Filippis et al. 
(c)2006 by the American Physical Society.) 



200a). 
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ground-state Eq and the zero, Mq and first, Mi, frequency moments calculated using the 
optical conductivity spectra: 



Mo = /f ™ Rea (fi) 
Ml = /f """" fiRea 



(114) 

(115) 



wher e flmax is the upper value of the frequency available from Ref. (IMishchenko et al. 
2mm . and 



Mn = 



TT 

2m* .11 



Rea (fi) 



;il6) 



Here m* is the polaron mass, the optical conductivity is calculated in units n^e^ j (mujQ), m* 
is measured in units of the band mass m, and the frequency is measured in units of cjq- The 
values of fimax are: f2max = 10 for a = 0.01, 1 and 3, f2max = 12 for a = 4.5, 5.25 and 6, 



^max = 18 for a = 6.5, 7 and 8. 



TABLE II Polaron parameters Mn^^M^^M^ 



(Reprinted with permission from ( 



Devreesd . 



obtain ed from the diagrammatic Monte Carlo results 



2006al ). ©2006, Societa Italiana di Fisica.) 



a 


^(DQMC) 


^*(DQMC) 


^(DQMC) 


^DQMC)/^ 


j^(DQMC) 


0.01 


0.00249 


1.0017 


1.5706 


0.634 


-0.010 


1 


0.24179 


1.1865 


1.5657 


0.65789 


-1.013 


3 


0.67743 


1.8467 


1.5280 


0.73123 


-3.18 


4.5 


0.97540 


2.8742 


1.5219 


0.862 


-4.97 


5.25 


1.0904 


3.8148 


1.5022 


0.90181 


-5.68 


6 


1.1994 


5.3708 


1.4919 


0.98248 


-6.79 


6.5 


1.30 


6.4989 


1.5417 


1.1356 


-7.44 


7 


1.3558 


9.7158 


1.5175 


1.2163 


-8.31 


8 


1.4195 


19.991 


1.4981 


1.3774 


-9.85 



The optical conductivity derived by DSG fjPevreese et al. 



1972al ) exactly satisfies the 
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TABLE III Polaron p arameters Mn^ M- i , Mn obtained within the path-integral approach (Reprinted 
with permission from (jPevreesd . l2006al ). ©2006, Societa Itahana di Fisica.) 



a 


^(DSG) 


j^*(Feynman) 


^^(DSG) 




T-i(Feynman) 


0.01 


0.00248 


1.0017 


1.5706 


0.633 


-0.010 


1 


0.24318 


1.1957 


1.5569 


0.65468 


-1.0130 


3 


0.69696 


1.8912 


1.5275 


0.71572 


-3.1333 


4.5 


1.0162 


3.1202 


1.5196 


0.83184 


-4.8394 


5.25 


1.1504 


4.3969 


1.5077 


0.88595 


-5.7482 


6 


1.2608 


6.8367 


1.4906 


0.95384 


-6.7108 


6.5 


1.3657 


9.7449 


1.5269 


1.1192 


-7.3920 


7 


1.4278 


14.395 


1.5369 


1.2170 


-8.1127 


8 


1.4741 


31.569 


1.5239 


1.4340 


-9.6953 



sum rule ( jPevreese et al. 1 . 119771 ) 



/-oo 

2m* ^ ^ 2 



Since the optical conductivity obtained from the DQMC results (IMishchenko et al. I . l2003l ) 



(117) 



is known only within a limit ed interval of frequenci es 1 < f2 < ilmax; the integral in Eq. 



(11161) for the DSG-approach (IDevreese et al 



interval as for the Monte Carlo results (IMishchenko et al. 



1972a) is ca: 



culated over the same frequency 



20031). 



The comparison of the resulting zero frequency moments M^ 



(DQMC) 




and Mn^^^'* with each 



other and with the value 7r/2 = 1.5707963... corresponding to the right-hand-side of the sum 
rule (HIT]) shows that the difference M^°^^^^ - M^ 



(DSG) 




on the interval a < 8 is smaller 



than the absolute value of the contribution of the "tail" of the optical conductivity for 
Q > f2max to the integral in the sum rule (11171) : 



/ Rea(°s^) (Q) dQ ^ - - M^^ 
I o 2 



(118) 



Within the accuracy determined by the neglect of the "tail" of the part of the spectrum for 
Vl > fimax) the contribution to the inte gral in the sum rule (1117 1 ) for the optical conductiv- 



ity obtained from the DQMC results (IMishchenko et al 



20031 ) agrees well with that for 
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the op tical conductivity found within the path-integral approach in Ref. (jPevreese et al. I . 
1972al ). Hence, the conclusi o n foll ows that the optical conductivity obtained from the DQMC 
results ( Mishchenko et al. . 200^ satisfies the sum rule flllTp within the aforementioned 
accuracy. 



H. Optical absorption spectra of continuum-polaron gas 



For the weak-coupling regime, which is realized in most polar semiconduc tors, the ground- 



state energy of a gas of interacting continuum polarons has been derived in ( iLemmens et al. 
19771 ) by introducing a variational wave function: 



I^ldb) = U\(p) Iv^ci) , 



(119) 



where \ipci) represents the ground-state many-body wave function for the electron (or hole) 
system, |0) is the phonon vacuum and U is a. many-body unitary operator. U defines the 
LDB-canonical transformation for a fermion gas interacting with a boson field: 



I i=l q J 



(120) 



where represent the p osition of the N const ituent electrons (or holes). The /q were 

' [3)- It may be emphasized that Eq. 



determined variationally (iLemmens et al. 



1977 



alt hough it appears like a straightforward generalization of the one-particle transformation 



in (ITomonaga . 



19471 ). constitutes — especially in its implementation — a nontrivial extension 



of a one-particle approximation to a many-body system. An advantage of the LDB-many- 
polaron canonical transformations introduced in (ILemmens et all 119771 ) for the calculation 
of the ground state energy of a polaron gas is that the many-body effects are contained in 
the static structure factor of the electron (or hole) system, which appears in the analytical 
expression for the energy. Within the approach the minimum of the total ground-state 
energy per particle for a polaron gas lies at lower density than that for the electron gas. 

The non-degenerate system of interacting polarons in polar doped insulators was analysed 
by Fratini and Quemerais (2000) using a simplified Feynman-type polaron model. In the 
low-density limit, the ground state of the many-polaron system is the Wigner lattice of 
polarons. With increasing density depending on the value of a, one of the following two 
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scenarios is possible: (i) the melting of the polaron Wigner lattice fo r a < a* and (ii) the 
dissociation of the polarons for a > a* ( iFratini and Quemeraisl . |2000| ). 



The LDB-canonical transformation has been fr uitfully applied in the theory of optical 



absorption spectra of many-polaron systems. In (ITempere and Devreese 



200 If ), starting 



from the LDB-many-polaron canonical trans 



wave function introduced in (ILemmens et al. 



or mat ion and the variational many-polaron 



19771 ). the optical absorption coefficient of 



a many-polaron gas has been derived. The real part of the optical conductivity of the 
many-polaron system is obtained in an intuitively appealing form, given by Eq. (11071) . 

This approach to the many-polaron optical absorption allows one to include the many- 
body effects to order a in terms of the dynamical structure factor S'(k, Vt — culo) of the 



electron (or hole) sy stem. The exper imental peaks in t he mid-infrared optica. 



spectra of cuprates (ILupi et al. 



19991 ). and manganites flHartinger et al 



absorption 



20041 ) have been 

adequately interpreted within this theory. The many-polaron approach descri bes the exper 



imental optical c ondu ctivity better t han the single-polaron approxi mations (Emin 



Gurevich et al. 



1993 



19621 ) . Note that in (ITempere and Devreese 1 . 1200 ll ) , like in (lEagles et al. 



19951 ). coexistence of small and Frohlich polarons in the same solid is inv olved. 



The optical conductivity of a many-polaron gas was investigated in (jCataudella et al. 



19991 ) in a different way by calculating the correction to the dielectric function of the electron 



gas, due to the electron-phonon interaction with variational parameters of a single-polaron 
Feynman model. A suppression of the optical absorption of a many-polaron gas as coni pared 



to the one-polaron optical absorption of Refs. (jPevreese 



1972 



Devreese et al. 



1972ah with 



increasing density has been found. Such a suppression is expected because of the screening 
of the Frohlich interaction with increasing polaron density. 



I. Ripplopolarons 



An interesting 20 sys t em consists of ele c trons on films of liquid He 



(IJackson and Platzma: 



3 



1981 



Shikin and Monarkha 



19731 ). In this system the elec- 



trons couple to the ripplons of the liquid He, forming "ripplopolarons". The effective 
coupling can be relatively large and self-trapping can result. The acoustic nature of the 
ripplon dispersion at long wavelengths induces the self-trapping. Spherical shells of charged 
particles appear in a variety of physical systems, such as fuUerenes, metallic nanoshells. 
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charged droplets and neutron stars. A particularly interesting physical realization of 
the spherical electron gas is found in multielectron bubbles (MEBs) in liquid helium-4. 
These MEBs are 0.1 fim - 100 /im sized cavities inside liquid helium, that contain helium 
vapour at vap our pressure and a nanomet e r-thick electron laye r, anchored to the surface 



of the bubble (Albrecht and Leiderer 



1987 



Volodin et al. 



19771 ). They exist result of 



equilibriu m between the surface tension of liqu i d he 



electrons (jSalomaa and Williamsl . 



1981 



Shikin 



19781). 



ium and the Coulomb repulsion of the 



Recently proposed experimental schemes to stabilize MEBs ( 



Silvera, 



200 



lated theoretical investigation of their properties (see e. g. (ITempere et al. 



have stimu- 



200ll)). The 



dynamical modes of MEB were described by considering the motion of the helium surface 
("ripplons") and the vibrational modes of the electrons together. In particular, the case 
when the ripplopolarons form a Wigner lattice was analyzed. The interaction energy be- 
tween the ripplons and the electrons in the multielectron bubble is derived from the following 
considerations: (i) the distance between the layer electrons and the helium surface is fixed 
(the electrons find themselves confined to an effectively 2D surface anchored to the helium 
surface) and (ii) the electrons are subjected to a force field, arising from the electric field 
of the other electrons. To study the ripplopolaron Wigner lattice at nonzero temperature 
and for any valu e of the electron-ripplon coupling, the variational path-integral approach 
(IFeynmaru . Il955l ) has been used. The destruction of the ripplopolaron Wigner lattice in a 
MEB occurs through the dissociation of ripplopolarons. Below a critical pressure (on the 
order of 10^ Pa) the ripplopolaron solid will melt into an electron liquid. This critical pres- 
sure is nearly independent of the number of electrons (except for the smallest bubbles) and 
is weakly temperature dependent, up to the helium critical temperature 5.2 K. This can 
be understood since the typical lattice potential well in which the ripplopolaron resides has 
frequencies of the order of THz or larger, which correspond to ^ 10 K. 

The new phase that was predicted in (ITempere et al\ . l2003al ). the ripplopolaron Wigner 
lattice, will not be present for electrons on a fiat helium surface. At the values of the pressing 
field necessary to obtain a strong enough electron-ripplon coupling, the fiat helium surfac e 
is no longer stable against long- wavelength deformations (IGor'kov and Chernikoval . Il973l ). 



Multi-electron bubbles, with their different ripplon d ispersion and the prese nce of stabilizing 



factors such as the energy barrier against fissioning (ITempere et al. 



2003bl ). allow for much 



larger electric fields pressing the electrons against the helium surface. The regime of A^, p, T 
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parameters suitable for the creation of a ripplopolaron Wigner lattice lies within the regime 



that wou 



d be achievab le in recently proposed experiments, aimed at stabilizing multielectron 



bubbles ( 1Silveral . l200ll ). The ripplopol aron Wigner lattice and its melting traiisition might be 



detected by spectroscopic techniques (IFisher et al 



11979 



for example the transverse phonon modes of the lattice ( 



Grimes and Adams 



Deville et al. 



198J). 



19791 ) probing 



J. Polaron scaling relations 



The form of the Frohlich Hamiltonian, Eq. (l28|) . in n dimensions is the same as in 3D, 
except that now all vectors are n- dimensional. In this subsection we take m = cuq = 1. In 

= 2\/2na/Viq^ . The interaction coefficient 



3D the EPI matrix element is well know n, \Vc . 



in n dimensions becomes ( jPeeters et all Il986l ) 



a 



;i2i) 



with Vn the volume of the n-dimensional crystal. 

The only difference between the model system in n dimensions and the model system 



in 3D is t 



lat no w one deals with an n-dimensional harmonic oscillato r. Directly fo. 



(IFeynman 



19551 ). the variational polaron energy was calculated in Ref. (jPeeters et al. 



lowing 



19861 ) 



_ n{v-wf _ r(^)" TOO e-* j±. 

2V2V{^) Jo ^/B^)' 



where 



w 



— w"^ 



2v^ 



(1 - e--) . 



(122) 



(123) 



In or der to facilitate a comparison of E for n dimensions with the Feynman result (IFeynmanl . 
19551) for 3D, 

(a 



3 (f — w)' 



roc 



it is convenient to rewrite Eq. (11221) in the form 



dt. 



n 



EnD (a) = T7 



a 



-.dt 



(124) 



(125) 



^ 2nr(f) /dTW' 

The parameters w and v must be determined by minimizing E. In the case of Eq. fll25p one 
should minimise the expression in the square brackets. The only difference of this expression 
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from the rhs of Eq. fll24p is that a is multiphed by the factor 



(126) 



2nr(f) 

This means that the minimizing parameters w and v in riD at a given a will be exactly the 
same as those calculated in 3D with the Frohlich constant chosen as a„Q!: 



Vnv> (a) = ^^3D (a„a) , w^d (a) = wsd (fln") 



(127) 



Comparing Eg. (11251) to Eg. fll24p. the followin g scaling relation (jPeeters and Devreesd . 



1987 



Peeters et al 



1986 



Xiaoguang et all Il985l ) is obtained 



EnD (a) 



ft 

-EsD (ana) 



(128) 



where a„ is given by Eq. fll26p . As discussed in Ref . (iPeeters et all . Il986l ) , the above scaling 
relation is not an exact relation. It is valid for the Feynman polaron energy and also for 
the ground-state energy to order a. The next-order term (i.e., a^) e. g. no longer satisfies 
Eq. (11281) . The reason is that in the exact calculation (to order a^) the electron motions 
in different space directions are coupled by EPI. No such coupling appears in the Feynman 
polaron model; this is the underlying reason for the validity of the scaling relation for the 
Feynman app roximation 
In Refs. 



1987 



Devreese et al. 



Peeters et al 



1977 



Peeters and Devreese 



1983a 



Peeters and Devreesd . 



19861 ). scaling relations were obtained also for the impedance function, 



Znj) {a; Q) = (a„a;; Q), the effective mass and the mobility of a polaron. In the impor- 



tant particular case of 2D, the scalin^ 



Peeters et al. 



1986 



g relations take the form (IPeeters and Devreese 



Xiaoguang et al. 



1987 



19851): 



E2D (a) 
Z2B («; ^) 



2^ f3n 



m. 



2D 



[a 



m 



3D 



/i2D (a) 



(m 



/^3D 



The scaling relati ons ( 



Monte Carlo results (ITitantah et al. 



Peeters and Devreese 



'3D 

Sir 



-a 



(129) 
(130) 
(131) 
(132) 



19871 ) can be checked for the path integral 



200 ll ) for the polaron free energy given in 3D and in 2D 



for a few values of temperature and for some sele cted values of a. They follow very closely 



the scaling relation of the form given by Eq. fll29p ((Devreese 



2006a). 
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IV. DISCRETE HOLSTEIN AND FROHLICH POLARON 



When the couphng with phonons increases the polaron radius decreases and becomes 
of the order of the lattice constant. Then all momenta of the Brillouin zone contribute to 
the polaron wave function and the effective mass approximation cannot be applied. This 
regime occurs if the characteristic potential energy Ep due to the local lattice deformation 
is compared or larger than the half-bandwidth D. The strong coupling regime with the 
dimensionless coupling constant 

A = § > 1 (133) 
is called the small or discrete lattice polaron. In general, Ep is expressed as 



(134) 



for any type of phonons involved in the polaron cloud. For the Frohlich int eraction with 



optic al phonons one obtains Ep ^ g^e^/vrK, where qd is the Debye momentum (lAlexandrov 



19961 ). For example, with parame ters appropriate for copper oxides gg ^ g ^ 5 and 



0.7 A ^ one obtains Ep ^ 



0.6 eV (^Alexandrov and Bratkovsky 



2000 



Eagles I . Il966l ). The 



exact value of Ac when the continuum (large) polaron transforms into the small one, depends 

the electron-phono n 



on the lattice structure, phonon frequency dispersions and the radius o 
interaction, but in most cases the transformation occurs around Ac — 1 (lAlexandrov 



2001ah 



A. Holstein model 



For comparison we briefly introduce a small polaron created by a short-range EPI, which 
i s known as the Ho lstein polaron. Its main features are revealed in the simple Holstein model 



(IHolstein 



1959bl ) of two vibrating molecules and the electron hopping between them. A 



simplified version of the model is defined by a two-site Hamiltonian describing the electron 
tunnelling between sites 1 ( "left" ) and 2 ( "right" ) with the amplitude t and interacting with 
a vibrational mode of an ion, placed at some distance in between, Figj9l 



H = t{c\c2 + 4ci) + Hph + H, 



e—phi 



(135) 



Here we take the position of an atomic level in the rigid lattice as zero, and q annihilates 
the electron on the left, « = 1, or on the right, i = 2, site. 
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a-polarised 



c-polarised | 




o 



"left" "right" 

FIG. 9 Electron tunnels between sites 1 ("left") and 2 ("right") with the amplitude t and interacts 
with c-axis or a-axis polarised vibrational modes of the ion, placed in between. 

The vibration part of the Hamiltonian in this toy model is 

1 d"^ kx"^ 
^'' = -2Md^^^—^ 

where M is the ion mass, k = MuJq is the spring constant, and x is the ion displacement. 

The electron-phonon interaction, H^-ph, depends on the polarization of vibrations. If the 

ion vibrates along the perpendicular direction to the hopping (in "c" -direction, FiglHD we 

have 

He-ph = fcx{c\ci + 4c2), (137) 

and 

He-ph = faX{c\ci - 4C2), (138) 

if the ion vibrates along the hopping ("a" direction). 

The wave-function of the electron and the ion is a linear superposition of two terms 
describing the electron on the "left" and on the "right" site, respectively, 

ip = [u{x)c\+v{x)cl]\0), (139) 

where |0) is the vacuum state describing a rigid lattice without the extra electron. Substi- 
tuting tp into the Schrodinger equation, Htp = Etp, we obtain two coupled equations for the 
amplitudes, 

{E - fa,cx - Hph)u{x) = tv{x), (140) 
{E ± fa,cx - Hph)v{x) = tu{x), (141) 
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where + and — in the second equation correspond to "a" and "c" polarised vibrations, 
respectively. There is the exact solution for the "c" -axis polarization, when a change in the 
ion position leads to the same shift of the electron energy on the left and on the right sites, 

u{x) = uxn{x), (142) 

v{x) = vxnix), (143) 

where u and v are constants and 

Xn{x) = l^-^^^^y' H^[ix-fJk)iM^^ (144) 

is the harmonic oscillator wave-function. There are two ladders of levels given by 

= -Ep±t + uJo{n + l/2) (145) 

with Ep = f^/2k. Here i^n(0 ^^e Hermite polynomials, and n = 0, 1, 2, 3, ... . Hence the 
c-axis single-ion deformation leads to the polaron level shift but without any renormalisation 
of the hopping integral t. In contrast, a-polarized vibrations with the opposite shift of the 
electron energy on the left and on the right sites, strongly renormalise the hopping integral. 
There is no simple general solution of the Holstein model in this case, but one can find it in 
two limiting cases, non-adiabatic, when t <^ uq and adiabatic, when t ^ ujq. 



1. Non-adiabatic Holstein polaron 

In the non-adiabatic regime the ion vibrations are fast and the electron hopping is slow. 
Hence one can apply a perturbation theory in powers of t to solve 



0, 



(146) 



E - faX - Hph -t \ u{x) 

-t E + faX - Hph J [ v{x) _ 

One takes t = in zero order, and obtains a two- fold degenerate ground state [m''''(x), f ''^(x)], 
corresponding to the polaron localised on the left (/) or on the right (r) sites. 



u (x) = exp 



{x + fa/ky 



v^{x) = 



(147) 



and 



u''{x) = 0, 
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v^{x) = exp 



MujQ 



{x - fa/kf 



(148) 



u{x) 




u\x) 









= a 








v{x) 









if{x) 



with the energy Eq = —Ep + ujq/2, where Ep = fH^k. The eigenstates are found as hnear 
superpositions of two unperturbed states, 

'1,'■(T^^ r n 

(149) 

Here the coefficients a and (3 are independent of x. The conventional secular equation for E 
is obtained, multiplying the first row by u\x) and the second row by v^{x), and integrating 
over the vibration coordinate, x, each of two equations of the system. The result is 

E-Eo i 
i E-Eo 

with the renormalised hopping integral 

i J^^dxu^{x)v''{x) 
t= JZdx\u^{xW 

The corresponding eigenvalues, E± are 

E± = uJo/2 - Ep± i. 



det 







(150) 



(151) 



(152) 



The hopping integral splits the degenerate level, as in the rigid lattice, but an effective 
'bandwidth' 2t is significantly reduced compared with the bare one. 



i = texp{—2Ep/ujQ). 



(153) 



This polaron band narrowing originates in a small overlap integral of two displaced oscillator 
wave functions u^(x) and v^(x). 



2. Adiabatic Holstein polaron 



In the adiabatic regime, when t ^ ujq, the electron tunnelling is fast compared with 
the ion motion. Hence one c an apply the Born-Oppenheimer adiabatic approximation 



(IBorn and Oppenheimer 1 . 119271 ) taking the wave function in the form 



u 






Ua 


[xY 




[x) _ 








y 






'x)_ 



(154) 
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Here Ua{x) and Va{x) are the electron wave functions obeying the Schrodinger equation with 
the frozen ion deformation x, i.e. 



Eaix) - faX -t 

-t Ea{x) + faX 



Uaix) 
Va{x) 



0. 



(155) 



The lowest energy level is found as 



Ea{x) = -^/{faXy+t^. 



(156) 



Ea{x) together with kx^ jl plays the role of a potential energy term in the equation for the 
'vibration' wave function, xi?^\ 

1 9^ kx^ 



t2 



Xix) = Ex{x) 



(157) 



2M dx^ 2 

Terms with the first and second derivatives of the electron wave-functions Ua{x) and Va{x) are 
small compared with the corresponding derivatives of x{^) the adiabatic approximation, 
so they are neglected in Eq.( ll57p . As a result we arrive at the familiar double- well potential 



problem, where the potential energy U{x) = kx'^/2 — (fax)'^ + has two symmetric 
minima, separated by a barrier. Minima are located approximately at 



Xm ^fa/k 



(158) 



in the strong-coupling limit, Ep 3> t, and the potential energy near the bottom of each 
potential well is about 

f/(x) = -i?, + ^fc^. (159) 



If the barrier were impenetrable, there would be the ground state energy level Eq = —Ep + 
ujq/2, the same for both wells. The underbarrier tunnelling results in a splitting of this level 
2t, which corresponds to a polaron bandwidth in the lattice. It can be estimated using the 
quasi-classical approximation as 



t oc exp 



-2 / p{x)dx 

'0 



(160) 



where p{x) = ^y2M[U{x) — Eq] ^ {Mky^'^\x — fa/k\ is the classical momentum 
Estimating the integral one finds the exponential reduction of the "bandwidth" , 



i oc exp(— 2£'p/a;o), 



:i6ii 



56 



which is the same as in the non-adiabatic regime. Holstein found corrections to this ex- 
pression up to terms of the order of which allowed him to improve the exponent and 
estimate the pre-exponential factor as 



EpUjQ 



IT 



Here = 2Ep/ujQ and 



A more accurate expression for t was obtained in (lAlexandrov et al. I . Il994al ) 



-9^ 



(162) 



(163) 



(164) 



where now g"^ = g'^{P — [ln(2A(l + /3))]/4A^}, and A = Ep/t. This expression takes into 
account the phonon frequency renormalisation, jS = ujq/ujq = 



'1 — 1/4A^, and the anhar- 
monic corrections of the order of 1/A^ to the turning point Xm. in Eg. (11601) . The term in 
front of the exponent in Eqs. (I162II164P differs from t of the non-adiabatic case, Eq. (11531) . 
It is thus apparent that the perturbation approach covers only a part of the entire lattice 
polaron region, A > 1. The upper limit of applicability of the perturbation theory is given by 



t < EpOJo- For the remainder of the region the adiabatic approximation is more appropri- 
ate. A much lower effective mass of the adiabatic small polaron in the intermediate coupling 
region compared with that estimated fro r n the p erturbation theory expression, Eq. (11531) . is 



1994"] ). The double well potential disappears at 



revealed in Eq. (11641) (lAlexandrov et al. 
A = Ac = 0.5, where the renormalised phonon frequency cJq is zero. 

The Holstein polaron model can be readily generalized for i nfinite lattices (section 



Similar models were used, f o r inst ance, in studies of dissipation (jCaldeira and Leggett 



1981 



Castro Neto and Caldeira 
systems. 



I991I ) and effects of decoherence in open quantum-mechanical 



B. Lang-Firsov canonical transformation 



The kinetic energy is smaller than the interaction energy as long as A > 1. Hence a 



self-consistent approach to the discrete (or lattic e) polaron problem with 



is possible with the "1/A" expansion technique (ILang and Firsov 



1962 



PI o f any range 



19631 ) on infinite 



lattices with any type of EPI conserving the electron site-occupation numbers, and any 
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phonon spectrum. The technique treats the electro n kinetic ener g y as a perturbation, and 
can be apphed for multi-polaron systems as well flAlexandrov I . Il983l ) (see IV1.B|) . It is 
based on the fact, known for a long time, that there is an analytical exact solution of the 
any — number polaron problem in the extreme strong-coupling limit, A ^ oo. Following 
Lang and Firsov one applies the canonical transformation to diagonalise the Hamiltonian, 
Eq. (fTT!) . The diagonahsation is exact, if t(m) = (or A = oo): 



H = e^He 



where 



5* = - ^ [ni(q)(iq - H.c. 



(165) 



(166) 



is such that = —S. The electron and phonon operators are transformed as 



C7 G Co G 



and (iq = e^d^e ^ . The result is 



and 



CiXi 



rfq = rfq - ^niM*(q) 



(167) 



:i68) 



where Xi = exp 



Eq^i(q)c^q - H.c. 



The Lang-Firsov transformation shifts the ions to 
new equilibrium positions. In a more general sense it changes the boson vacuum. As a 
result, the transformed Hamiltonian takes the following form 



where 



aij = t{m - n)5ss'XlXj 



is the renormalised hopping integral depending on the phonon operators, and 
Vij = Vc{m - n) - — ^ |7(q)|2^qCos[q ■ (m - n)] 



(169) 



(170) 



(171) 



is the interaction of polarons, Vij = v{m — n), comprising their Coulomb repulsion and the 
interaction via the lattice deformation. In the extreme infinite-coupling limit, A ^ cxd, we 
can neglect the hopping term of the transformed Hamiltonian. The rest has analytically 
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determined eigenstates and eigenvalues. The eigenstates |A^) = |nj,nq) are sorted by the 
polaron n^s and phonon occupation numbers. The energy levels are 

1 X 

E = + Ep)^ni + - ^VijUiUj + ^u^{n^ + 1/2), (172) 

i i^j q 

where rij = 0, 1 and = 0, 1, 2, 3, ....oo. 



1. "1/X" expansion and polaron band 



The Hamiltonian H in zero order with respect to the hopping describes localised polarons 
and independent phonons, which are vibrations of ions relative to new equilibrium positions 
depending on the polaron occupation numbers. The middle of the electron band is shifted 
down by the polaron level-shift Ep due to the potential well created by lattice deformation. 



Importantly the phonon frequencies remain unchanged in this limit at any po 



n. At finite A and n there is a softenin g 



1992a 



Alexandrov and Capellmann 



of p 



l onons 6uJo of the order of uJnnlX^ (lAlexandrov 



1991 



Alexandrov et al. 



(the initial paper on the phonon renormalisat ion (lAlexandrov and 



dieting 5uJn oc 1/A was subsequently corrected (lAlexandrov 



1992 



aron density , 



Gobel et al 



1992a 



Capellmann 



Alexandrov et al. 



1994h 



1991 ) pre 



1992 



Gobel et al. 



19941 )). Interestingly the optical phonon can be mixed with a low- 



laronic plasmon forming a new excitation, "plasphon" , which was proposed in (lAlexandrov 



requency po- 



1992a) as an exp 



(IRietschel et al. 



anation of the anomalous phonon mode splitting observed in cuprates 



19891). 



Now let us discuss the 1/A expansion. First we restrict the discussion to a single-polaron 
problem with no polaron-polaron interaction. The finite hopping term leads to the polaron 
tunnelling because of degeneracy of the zero order Hamiltonian with respect to the site 
position of the polaron. To see how the tunnelling occurs we apply the perturbation theory 
using 1/A as a small parameter. The proper Bloch set of A^-fold degenerate zero order 
eigenstates with the lowest energy {—Ep) of the unperturbed Hamiltonian is 



'^'0) = ;;^E^-exp(^k-m)|0), 



;i73) 



where |0) is the vacuum, and is the number of sites. By applying the textbook pertur- 
bation theory one readily calculates the perturbed energy levels. Up to the second order in 
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the hopping integral they are given by 



-Ep + ek-J2 



l(k,0|E^,j^»i4cj|k^r^q) 



(174) 

k',n, -q^q'^q 

where |k',nq) are the exited states of the unperturbed Hamiltonian with one electron and 
at least one real phonon. The second term in this equation, w hich is linear with r espect to 
the bare hopping t(m), describes the polaron-band dispersion flAlexandrov 1 . 12000| ). 

ek = ^t(m)e-^''('")exp(-zk-m), (175) 

m 

where 

q 

is the band-narrowing factor at zero temperature. 

The third term, quadratic in t(m) , yields a negative k.- independent correction to the 
polaron le yel-shift of the order of 1/A^, and a small correction to the pqlaron band dispersion. 



Eq. fll75p (lAlexandrov 



1992a 



Gogohn 



1982 



Kudinov and Firsov 



19971 ). The correction 



to the level shift is due to polaronic hops onto a neighbouring site with no deformation 
around it. As any second order correction this transition shifts the energy down by an 
amount of about —t'^{m)/Ep. It has little to do with the polaron effective mass and the 
polaron tunneling mobility because the lattice deformation does not follow the electron. 
The polaron hops back and forth many times (about e^^) "waiting" for a sufficient lattice 
deformation to appear around neighbouring site n. Only after the deformation around the 
neighbouring site is created does the polaron tunnel onto the next site together with the 
deformation. 



2. Temperature effect on the polaron hand 

Let us now analyse the temperature dependence of the polaron bandwidth, which is 
determined by the average of the multiphonon operator, Eq. fll70p . 

((X/X,)) = \[{{eM<W\ - H.c] exp[n,(q)rfq. - H.c.])). (177) 
q 

Here the double angular brackets correspond to quantum as well as statistical averages of 
any operator A with the Gibbs distribution. 



{n~E„)/Ti 



,\A\v)=Tr{e^''-^^'^A}, 



:i78) 
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where Q is the thermodynamic potential and |z/) are the eigenstates of H with the eigenvalues 
Ey. An operator identity exp(A + B) = exp{Aexp{B)exp{—[A,B]/2) is instrumental. It is 
applied when the commutator [A, B] is a number. The identity allows us to write 



,K(q)4-H.c.] [«,{q)dq-H.c.] ^ a*4 -adq -|q|V2 



Uii<ih*i (q)-"i (q)t':i (q)]/2 



(179) 



Quantum and statistical averages are calculated by expanding the exponents in the trace as 
((e"*''e-"'^)) = il-p)J2 J^P^'^-^rPf^NiN - 1) x ... x {N - n + 1), (180) 

Ar=On=0 ^ 

where we dropped the phonon and site quantum numbers for transparency. Here p = 
exp(— cUqj^/T), so that a sin gle-mode phonon partitio n function is Zp^ = 1/(1 —p). Eq. fllSOp 
can be written in the form (ILang and Firsov I . Il962l ) 

N I |2n 



n=0 ^ ^ N=0 

Taking the sum over A^, ^^=0^*^ = V(l differentiating it n times yield n\ in the 

numerator, after which the series over n turns out equal to 



((e"*'^^e-"'^)) = e" 



(182) 



where n^^ = [exp(a;q/T) — 1] ^ is the Bose-Einstein phonon distribution function. Collecting 
all multipliers one finally obtains 



((a,,)) = T(m - n)S,,, exp (^-^ ^ |7(q)|2[l - cos(q ■ m)] coth ^ j , 



(183) 



with the zero-temperature limit given by Eq.f ll75l) . 

The small-polaron band i s exponentially narrow. H e nce, on e can raise a concern about 



its existence in real solids (jPe Mello and Ranninger 



1998al ). At zero temperature the 



perturbation term of the transformed Hamiltonian conserves the momentum because all 
off-diagonal matrix elements vanish. 



(k,0|^a„clc,|k',0) = 



;i84) 
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if k 7^ k'. The emission of a single high-frequency phonon is impossible for any k because 
of the energy conservation. The polaron half-bandwidth is exponentially reduced, 



w 



De 



;i85) 



and it is usually less than the optical phonon energy cuq {g'^ is about DX/ujq). Hence, there 
is no damping of the polaron band at T = caused by optical phonons, no matter how 
strong the interaction is. The phonons "dress" the electron and coherently follow its mo- 
tion. However, at finite temperatures the simultaneous emission and absorption of phonons 
is possible. Moreover the polaron bandwidth shrinks with increasing temperature because 
the phonon-averaged hopping integrals depend on temperature , Eg. (11831) . For high tem- 
peratures, T ^ uJo/2, the band narrows exponentially as ~ De~^^^°, where 

q 

On the other hand, the two-phonon scattering of polarons becomes more important with 
increasing temperature. One can estimate this scattering rate by applying the Fermi-Dirac 
golden rule. 



\q,q' 



(186) 



where the corresponding matrix element is 



M. 



qq 



^ (k + q - q', - 1, riq/ 1| aijcjcj |k, nq, nq/) . 



For simplicity we consider the momentum independent 7(q) = 7o and cUq = cjq. Expanding 
o"jj-operators in powers of the phonon creation and annihilation operators one estimates the 
matrix element of the two-phonon scattering as Mqq/ ^ A^~^w7Qy^?2^Jr^r+Ty. Using this 

l/2w one 

(187) 



^W7^V^raq(raq/ + 1 

estimate and the polaron densi t y of s tates (DOS), Pp(0 = X]k^(^ " ^k) 



obtains (lAlexandrov and Mott 



19951 ) 



T 



where = [exp{ujQ/T) — 1]"^ is the phonon distribution function. 

The polaron band is well defined, if 1/r < w, which is satisfied for a temperature 
range T < Tmin ~ ^o/^'^lo about half of the characteristic phonon frequency for rele- 



vant values of 7q. At higher tem peratures t 



ne mco 



inates in the polaron dynamics (iHolstein 



1959bl : 



lerent thermal activated hoppin 



Lang and Firsov 



1962 : 



g dom- 



Sewell 



1958 : 
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Yamashita and Kurosawa 



When the optical phonon freq uencies are excep t ional 



19581 ). and the polaron states are no longer the Bloch states. 

y high (i.e. about lOOOK as in high- 



temperature superconductors (IRietschel et al. I . Il989l )) lattice polarons are in the Bloch 



states in the relevant range of temperatures, where the Boltzmann kinetic theory with renor- 
malised energy spectrum is applied. 



C. Discrete Frohlich polaron at strong coupling 



T he narrowing of the band and the p olaron effective mass strongly depend on the range of 



EPI ( lAlexandrov 



1996 



Eagles I . Il969l ). Let us compare the small Holstein polaron (SHP) 



formed by the zero-range EPI and a small polaron formed by the long-range (Frohlich) 



interaction, which we refer to as t he small or discrete Frohlich po 



aron (SFP). We use the 



19991) 



real-space representation of H^-ph ( lAlexandrov and Kornilovitch 

He-ph = ^ /(m - n){nfii. 

with the normal coordinate at site n 

= ^(2A^MWq)^l/2giq-n^^ ^ ^ 

q 

and the force between the electron at site m and the normal coordinate .^n, 

/(m) = iV-i5^7(q)(M^2)'/V^-'". 



(188) 



;i89) 



(190) 



In general, there is no simple relation between the polaron level-shift Ep and the exponent 
of the mass enhancement. This relation depends on the form of EPI. Indeed for EPI with 
a single dispersionless phonon mode, tUq = uq, one obtains 



E„ 



2Muj, 



m 



and 



;i9i) 



(192) 



where a is the primitive lattice vector. In the nearest-neighbour approximation the effective 
mass renormalisation is given by m*/m = e^^ where 1/m* = d'^e^^/dk^ at ^ is the 
inverse polaron mass. If the interaction is short-ranged, /(m) = K(5m,o (the Holstein model), 
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FIG. 10 A one-dimensional model of the lattice polaron on chain m interacting with displacements 
of all ions of another chain n (t is the polaron hopping integral along chain m). 

then g'^ = Ep/ujQ. Here k is a constant. In general, one has g'^ = '-fEp/uo with the numerical 



coefficient 



l-E„./(m)/(m + a) 



^ En/^(n) ' ^ ^ 

which might be less than 1. To estimate 7 let us consider a one-dimensional chain model 

with the long-range Coulomb interaction between the electron on one chain (m) and ion 

vibrations of another chain (n), po larised in the direction perpendicular to the chains 

(lAlexandrov and Kornilovitch I . Il999l ). FigJTOl The corresponding force is given by 



/(m - n) 



(194) 



(|m -n|2 + 1)3/2- 

Here the distance along the chains |m — n| is measured in units of the lattice constant a, the 
inter-chain distance is also a, and we take a = 1. For this long-range interaction we obtain 
Ep = 1.27KV(2Mcjg), g^ = 0.49kV(2Mcu[|) and g^ = 0.39Ep/uJo . Thus the effective mass 
renormalisation in the non-adiabatic regime is much smaller than in the Holstein model, 
roughly as m*gpp oc (^^^//p)^^^ in units of m. An analytical solution of a two-site single 
electron system interacting with many vibrating ions of a lattice via a long-range Frohlich 



EPI found that SFP is also several 



regime (lAlexandrov and Yavidov 



orders of magnitude lighter than SHP in the adiabatic 



20M). 



Another interesting point is that the size of SFP and the length, over which the distortion 
spreads, are different. In the strong-coupling limit the polaron is almost localised on one 
site m. Hence, the size of its wave function is the atomic size. On the other hand, the 
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ion displacements, proportional to the displacement force /(m — n), spread over a large 
distance. Their amphtude at a site n falls with the distance as |m — n|~'^ in the one- 
dimensional model, FigJTUl The polaron cloud (i.e. lattice distortion) is more extended 
than the polaron itself. Such polaron tunnels with a larger probability than the Holstein 
polaron due to a smaller relative lattice distortion around two neighbouring sites. For a 
short-range EPl the entire lattice deformation disappears at one site and then forms at its 
neighbour, when the polaron tunnels from site to site. Therefore 7 = 1 and the polaron 
is very heavy already at A ~ 1. On the contrary, if the interaction is long-ranged, only a 
fraction of the total deformation changes every time the polaron tunnels from one site to its 
neighbour, and 7 is smaller than 1. 



D. Effect of dispersive phonons 

A lighter mass of SEP compared with the nondispersive SHP is a generic feature of 
any dispersive electron-phonon interaction. For example, a short-range interaction with 



dispersive acoustic phonons (7(q) oc l/q-*^^^, oc q) also 



strong-coupling regime compared with SHP flAlexandrov 



eads to a lighter polaron in the 



2OO3I ). Even within the Holstein 



model with the local (intramolecular) EPl the dispersion of phonon frequencies is a vital 
ingredient since nondispersive phonons mi ght lead to a divergent site jump probability of 
polarons (lYamashita and Kurosawa I . Il958l ). Importantly the comprehensive studies of the 
molecular Holstein Hamiltonian, in which the dispersive features of the phonon spectrum 
are taken into account , found mu c h low e r values of the polar on mass compared with the 



non-dispersive model (jZoli 



1998 



2000 



Zoli and Das 



200J). In those studies the 1/A 



perturbation theory based on the standard Lang-Eirsov (LE) and the variational (modified) 
MLE transformation of the molecular Holstein Hamiltonian with dispersive phonons has 
been applied for ID, 2D and 3D lattices including the second-order corrections as in Eg. (11 741) . 



H 



-t > clcj + gujo > Hi 



{di + d\) + ^uj^{d\d^ + l/2), 
q 



(195) 



where (iq is the Eourier transform of di. The phonon dispersion has been modeled using 
the intermolecular first neighbours force constant, Mcuf, which yields, for example, in ID 
case = uJq/2 + cj^ + {uJq/^ + uj^ujlcosq + oj'lY^'^. MLE improves the convergence of 
the 1/A perturbation series by introducing a suitable variational parameter Aq in the LE 
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transformation, Eq. fll66p . as (IChatterjee and Das 1 . 120001 : IChatterjee et al. 1 . 12003| ) 



S =^ Smlf 



q,k 



Cu_L„Ck(dq - dl 



(196) 



As a result the polaron mass converges to much lower values when the phonon dispersion 
is introduced, in particular in the adiabatic regime. There is a continuous mass enhance- 
ment whose abruptness is significantly smoothed for the largest values of the phonon dis- 
persion, similar to the SFP discussed above. While a phase transition is ruled out in the 



single-ele ctron Holstein Hami 



strength (IGerlach and Lowen 



t onian, where the gr ound state energy is analytic in the EPI 



1991 



Lowen 



19881 ). a crossover from more itinerant to more 



self-trapped be h aviour may be ident i fied a s a rather sudden e yent i n the adiabatic regime 



(lEagles 



20031 ). 



1969 



Emin and Holstein 



1976 



Kalosakas et al. 



1998 



Zheng and Avignon 



E. All-coupling discrete polaron 



1. Holstein model at any coupling 



During past twenty years significant efforts were directed towards the extension of the 
weak and strong-coupling perturbation lattice polaron theories to the intern i ediate region of 



the re 



1975 



evant pa r amet 



Gogolin 



ers, A ~ 1, and cuo/t ~ 1. It was argued (lAlexandrov 



1992a 



Firsov I . 



19821 ) that the expansion parameter is actually l/2zA^, so the analytical 



strong-coupling expansion in powers of 1/A might have a wider region of applicability than 
one can expect using simple physical arguments (i.e. A > 1). However, it has not been clear 
how fast the expansion converges. 

Kudinov and Firsov (1997) developed the analytical approach to the two-site Holstein 
model by the use of the expansion technique, which provides the electronic and vibronic 
terms as well as the wave functions and all correlation functions in any order of powers of 
t. They have found the exponential reduction factor in all orders of the 1/A perturbation 
expansion. On the other hand, the corrections to the atomic level were found as small 
as 1/A^ rather than exponential in agreement with the conventional second-order result, 
Eq.f ll75p . Chatterjee and Das (2000) studied the same problem for any coupling within the 
perturbative expansion combined with MLF, Eq. (11961) . and MLF with a squeezing canonical 
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transformation (iHang I . Il988l ). exp(S'), where S = a{didi — d\d\). Using two variational 
parameters introduced by MLF and squeezing transformations allows for very good conver- 
gence of the 1/A perturbation series even in the near-adiabatic regime, uj^/t > 0.5, where 
the conventional 1/A expansion shows poor convergence. These studies also showed that the 
region of the parameters of the Holstein model, where neither week nor strong-coupling per- 
turbation analytical methods are applicable, is rather narrow. A semi-analytical approach 
to the solution of two coupled differential equatio ns. Ea.(|140l 1411). o f the Holstein model in 



the whole parameter space has been proposed in (IHan et al. 



2OO2I ) based on the coherent- 



state expansion of u{x),v{x). These authors obtained the recursive relations among the 
expansion coefficients, allowing for highly accurate numerical solutions, which agree well 
with those by MLF method in the weak- and strong-coupling regimes. The deviation from 
the MLF solution in the intermediate-coupling regime implies that MLF misses some higher 
order correlation term s. A continued fraction an alytical solution of the two-site Holstein 



model was der ived by (jCapone and Ciuchi 



2OO2I ) as for a related model in quantum optics 



(jSwain I . Il973l ). In practice it also requires some truncation of the infinite phonon Hilbert 
space. Finally all Green's functions for the two - site H olstein-Hubbard model (HHM) were 
derived in terms of continued fractions (IBerciu 1 . 120071 ). 

Numerical results obtained by different methods actually show that the ground 
state energy (about —Ep) is not very sensitive to the parameters, while the effec- 
tive mass and the bandwidth strongly depend on the polaron size, EPI range, and 
the adiabatic ratio, oJo/t. Several methods exist for nur aerical simulations of lat 



tice po la rons. 



Barisic 



2004 



They include exact 



Capone et al 



Kongeter and Wagner 



Stephan 



1996 



1990 



1997 



diagonalisatio n 



Fehs. 



^e et al. 



Marsiglio 



Wellein et al. 



1993 



(ED ) (lAlexandrov et al. 



1995 



19951: 



1994a 



Fehske and Trugman 



Ranninger and Thib 



19961) , the glob a 



other advanced variational met hods ( 



Monte-Car l o (QMC) algorithms 



1984 



2004 



1985 



Fradkin and Hirsch 



2005b 



Ku et al 



f Berger et al. 



1983 



-local (GL) ( 



I995I: 



2002 



Romero et al. 



Trugman et al. 



m 



2007 



1992 



200 



De Raedt an 



Hirsch and Fradkin 



Hohenadler and von der Linden 



group (DMRG) (jJeckelmann and White 



1998 



1982 



d Lagendiik 



1983 



19981) and 



, quantuna 



1982 



1983 . 



Hohenadler et al. 



2007cf). density rn a trix renormalization 



Jeckelmann et al. 



1999f). contiriuous-time QMC f|Alexandrov and Kornilovitch 



Kornilovitch 



1998 



1999 



Spencer et al. 



1999 



1999 



Zhang et al. 



Hague et al. 



2006a 



2OO5I ). and diagrammatic QMC (IMacridin et al. 
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2004 



Mishchenko et al. 



2000 



Prokof 'ev and Svistunov 



19981 ). The methods vary in 



accuracy and versatihty, and, combined together, can provide all the polaron proper- 
ties of interest in the entire space of model parameters. On the other hand, ED suf- 
fers from the necessary truncation of the phonon Hilbert space, especially at strong cou- 
plings and low phonon frequencies (even then, the total Hilbert space is huge, reducing 
the number of sites and leading to poor momentum resolution), DMRG cannot easily 
handle long-range interactions, diagrammatic QMC and ED are inconvenient in calcu- 



(Fehske and Trueman . 




2007; 


H 


ohenadler and von der Linden 


Mishchenko and Nagaosa. 


2007 


). In numerical analysis of pol 



2007c 



Kornilovitch 



2007 



proach is needed where each method is employed to calculate what it does best. 

Until recently most numerical studies were performed on the Holstein model 
(i.e with zero-range EPI). Reliable results for the in t ermed i ate region we r e ob 



tained using ED of vibr a ting c 



De Mello and Ranninger 
199,4 



Marsiglio 



19971. 



usters fjAlexandrov et al. 



1998a 



Fehske et al 



Ranninger and Thibblin 



1992 



1995 



Stephan 



994i 



Capone et al. 



Fehske and Trugmar 



19961 : 



Wellein et al. 



1997 



2007 



19961). 



Taking as a measure of the polaron kinetic energy the correlation function tg// = (~^(cjc2 + 
c\ci)) (here ci,2 are annihilation operators on t he 'left' and 'right' sites of the Holstein m odel) 



one might doubt of the Lang-Firsov approach (IDe Mello and Ranninger 



1997 



1998al ). since 



this correlation function is much larger than the small polaron bandwidth. However, ap- 
plying the 1/A expansion up to the second order in t one obtains the numerical tp ff very 



close to the pertur bation tip in the strong-coupling regime, A > 1, (lAlexandrov 



Firsov et al. 



2000; 



19991) 



At 

UJ 



t 

A' 



(197) 



teff - tlF = -texp 

with A = 2Ep/t. Here only the first exponential term describes the true coherent tunnelling, 
while the second term describes the correction to the middle of the polaron band owing to 
the virtual "back-forth" transitions to the neighbouring site, as discussed above fllV.B.ll) . 
The m ain contribut i on to tpjf conies from the second-order term lowering the middle of the 



band (lAlexandrov 



1992a 



2000 



Gogolin 



1982 



Kudinov and Firsov 



from the polaron-transport related first term (see also (IDe Mello and Ranninger 



1997), rather than 



1998b 



19991 )). Comparing the analytical expression, Eq. (ll97p with the numerically calculated tg// 



one confirms that the Holstein-Lang-Firsov approach is asymptotically exact both in the 
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non-adiabatic and adiabatic regimes (lAlexandrov 1 . 12000| ). if the second-order correction is 
taken into account. 



The numerica 



(lAlexandrov et al. 



diago nahsation of the two-site-one-electron Holstein model 



1994al ) shows that the first-order term of the 1/A perturbation 



theory describes well the polaron bandwidth in the non-adiabatic regime for all values 
of the coupling constant. There is no agreement in the adiabatic region, where the first 
order perturbation expression, Eg. (11531) . overestimates the polaron mass by a few orders 
of magnitude. A poor convergenc e of the p e rturba tion expansion is due to appearance of 
the familiar double- well potential (IHolstein I . Il959bl ) in the adiabatic limit. The tunnelling 
probability is extremely sensitive to the shape of this potential. The splitting of levels for 
the two-site cluster is well described by the Holstein quasi-classical formula generalised 
for the intermediate coupling in Eq. fll64p . While the small Holstein polaron is only a few 
times heavier than the bare (unrenormalised) electron in a wide range of the coupling for a 
moderate adiabatic ratio un/t > 1, it becomes very heavy in the adiabatic regime and for 
the strong coupling (lAlexandrov et al. I . Il994al ). 



2. Holstein polaron in infinite lattices 



A number of other independent numerical results proved that "by the use of the Holstein 
approximation and the canonical Lang-Firsov approach with appropriate corrections, one 
obtains an excellent estimate of the coherent bandwidth in both adiabatic and non-adiabatic 
(strong- coupling) regi r nes'' 



(VED) (IBonca et al. 



1999 



ehske et al. 



1997) 



Fehske and Trugman 



I I part icular the elaborate variational ED 



20071 ) provides an exact numerical solu- 



tion of the Holstein crystal model, Eq. (ll95p . in any dimension. In contrast to finite-lattice 
ED, it yields the ground-state energy which is a variational bound for the exac t energy in the 



thermodynamic limit. Figure fTTl shows the effective mass computed by VED (IBonca et al. 



19991 ) in comparison with global-local (GL) and DMRG methods, m* is obtained from 
m/m* = {l/2t)d'^Ek/dk'^, k ^ 0, where mo = l/(2t) is the rigid band mass (here the lattice 
constant is a = 1). In the intermediate coupling regime where VED gives an energy accuracy 
of 12 decimal places, one ca n calculate the effective mass accurately (6-8 decimal places) 
(IFehske and Trugman 1 . 120071 ). 

In Fig. [TT] the parameters span different physical regimes including weak and strong 
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IG. 11 Logarithm of the p olaron effective mass in ID Holstein model as a function of g (after 



(jFehske and Trugman 



20071 )). VED results (full lines) were obtained operating repeatedly L = 20 



times with the off-diagonal pieces o f the Holstein Hamilt onian. For comparison GL (global-local) 



results (dashed lines) are included (IRomero et al 



(jJo/t, are DMRG results (jJeckelmann and White 



1993) 



19981 ') ■ Open symbols, indicating the value of 



coupling, and low and high phonon frequency. There is good agreement between VED and 
GL away from strong coupling and excellent agreement in all regimes with DMRG results. 
DMRG calculations are not based on finite-/^ calculations due to a lack of periodic boundary 
conditions, so they extrapolate the effective mass from the ground state data using chains 
of different sizes, which leads to larger error bars and demands more computational effort. 
There is no phase transition in the ground state of the model, but the polaron becomes 
extremely heavy in the strong coupling regime. The crossover to a regime of large polaron 
mass is more rapid in adiabatic regime, i.e. at small ujo/t- 

Unifying systematization of the crossovers between the different polaron behaviors in 
one dimension was obtained in terms of quantum to classical, weak coupling to strong 
coupling, adiabatic to nonadiabatic , itiner ant to self-trapped polarons and large to small 
polarons (10. Barisic and S. Barisic I . |2008| ). It was argued that the relationship between 
various aspects of polaron states can be specified by five regimes: the weak-coupling regime, 
the regime of large adiabatic polarons, the regime of small adiabatic polarons, the regime of 
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small nonadiabatic (Lang-Firsov) polarons, and the transitory regime of small polarons for 
which the adiabatic and nonadiabatic contributions are inextricably mixed in the polaron 
dispersion properties. Holstein polarons in three-dimensional anisotropic lattices have been 
recently studied using a variational ED technique, which provide s highly accurate results 
for the polaron mass and polaron radius (lAlvermann et al. I . l2008l ). Varying the anisotropy 
( lAlvermann et al. 1 . 12008| ) demonstrated how a polaron evolves from a one-dimensional to a 
three-dimensional quasiparticle and showed that even the local Holstein interaction leads to 
an enhancement of anisotropy in charge- carrier motion. 

For analyzing all-coupling polarons in more complex lattices the continuous-time path- 
integral quantum Monte -Carlo algorithm (CTQM C) is ideally suited. The algorithm is 



formulated in real space (IKornilovitch 



1998 



20071 ) and based on the analytical integration 



over phonon degrees of freedom introduced by Feynman (1955) and on an earlier numer- 
ical implementation in discrete time by De Raedt and Lagendijk (1982,1983,1984,1985). 
CTQMC introduced two critical improvements. Firstly, formulation in continuous imagi- 
nary time eliminated errors caused by the Trotter slicing and made the method numerically 



exact for any st rength of EP 



imaginary time (IKornilovitch 



. Secondly, introduction of tw i sted boundary conditions in 



1998 



Kornilovitch and Pike 



19971 ) enabled calculation of 



polaron effective masses, spectra and the densities of states (DOS) in any dimensions in 
infinite lattices. 

The polaron action, obtained by analytical integration over phonon degrees of freedom, is 
a functional of the polaron path in imaginary time r(r). It is given by the following double 
integral 



AW)] = / drdr'-^/^(e-(^/-l-'l) + e-(^/-l-'l))$o[r(r),r(r')] 

7\r) 

+ I rfrdr'"-^e-^(^-^')(<|.Ar[r(r),r(r')]-<fo[r(r),r(r')]):, (198) 

$Ar[r(r), r(rO] = /m[r(r)]Wr[r(r')] :, (199) 

m 

where the vector Ar = r(/3) — r(0) is the difference between the end points of the polaron 
path, P = t/T, u) = ujo/t, and /m(n) = /in(n)/K (see also Eq. fll94p ). From this starting 
point, the polaron is simulated using the Metropolis Monte-Carlo method. The electron 
path is continuous in time with hopping events (or kinks) introduced or removed from the 
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path with each Monte-Carlo step. From this ensemble, various physical properties may be 
computed, in particular, the ground state energy, the number of phonons in the polaron 
"cloud" and the polaron band energy spectrum. 



ek — Co = — ]im ^ ln(cos(k ■ Ar)) :, 



(200) 



where k is the quasi momentum. By expanding this expression in small k, the z-th component 
of the inverse effective mass is obtained as 



^ = lim i((Ar,)2) : 



(201) 



Thus the inverse effective mass is the diffusion coefficient of the polaron path in the limit of 
the infinitely long "diffusion time" /3. C TQMC has enabled accurate analysis of models with 



long-range electron-phonon i nteractions (lAlexandrov and Kornilovitch 



1999 



2006a 



19991 



Spencer et al. 1 . 120051 ) and a model with anisotropic electron hopping (IKornilovitch 



Hague et al. 



All numerical results confirm gross polaronic features well understood analytically by 



( iHolstein I . Il959bl ) and others both in the non-adiabatic and adiabatic regimes. A great 
power of numerical methods is the ability to calculate an entire polaron spectrum, and 
the polaron DOS p{E) = ^^6{E — E]^ + Eq), in the whole parameter space. The 
coherent part of the spectrum, ek, possesses an interesting property of flattening at large 



lattice momenta in th e adiabatic limit, t ^ Un (iFehske and Trugman 



1998. 


1999; 


Stephan . 


1996: 



Wellein et al. 



1996 



2007 



Wellein and Fehske 



Romero et al. 



19971 ). In the weak- 



coupling limit, the flattening can be r eadily understood as hybridiz ation between the bare 



electron spectrum and a phonon mode (ILevinson and Rashba 



19731 ). The resulting polaron 



dispersion is cosine-like at small k and fiat at large k. As a result the polaron DOS shou ld be 



peaked cl ose to the top o 



CTQMC flKornilovitch 



the polaron band. Exact VED (IFehske and Trugman 1 . 120071 ) and 



20071 ) calculations have confirmed that this is indeed the case for 
the short-range EPI. At small ujq, DOS develops a massive peak at the top of the band. With 
increasing uq, the polaron spectrum approaches the cosine-like shape in full accordance with 
the Lang-Firsov non-adiabatic formula. The respective DOS gradually assumes the familiar 
shape of the tight-binding band with renormalised hopping integrals. However, at small-to- 
moderate uo/t in two and three dimensions, the bottom half of the polaron band contains a 
tiny minority of the total number of states, so that the system's responses will be dominated 
by the states in the peak. 
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In the long-range model (11940 DOS is much closer to the tight-binding shape than the 
Holstein DOS at the same parameters (IHague et al. I . l2006al ). The polaron spectrum and 
DOS is another manifestation of the extremity of the Holstein model. Long-range EPI 
removes those peculiarities and makes the shape of polaron bands closer to the 1/A expansion 
results. 



3. Discrete Frohlich polaron at any coupling 

As discussed above (lIV.Cp the lattice polaron mass strongly depends on the radius of EPI. 
So does the range of the applicability of the analytical 1/A expansion theory. The theory 
appears almost exact in a wide region of the Frohlich EPI. Eg. (11941). for which the exact 



polar on mass was calculated with CTQMC algorithm in (lAlexandrov and Kornilovitch 



19991 ). 

At large A (> 1.5) SEP was found to be much lighter than SHP in agreement with the 
analytical results ( 1IV.CI) . while the large Frohlich polaron (i .e. at A < 1) was heavier than 



the 



arge Holstein polaron with the same binding energy (lAlexandrov and Kornilovitch 
19991 ). The mass ratio m*pp/m*fjp is a non-monotonic function of A. The effective mass of 
the small Frohlich polaron, m^p(A) is well fitted by a single exponent, which is e^'^^'^ for 
= t and e^'^^ for ujq = 0.5t. The exponents are remarkably close to those obtained with 
the Lang-Firsov transformation, e°'^^^ and e^'^^^, respectively. Hence, in the case of the 
Frohlich interaction the transformation is perfectly accurate even in the moderate adiabatic 
regime, ujo/t < 1 for any coupling strength. It is not the case for the Holstein polaron. If the 
interaction is short-ranged, the same analytical technique is applied only in the non-adiabatic 
regime Uo/t > 1. 



An important question about polaron properties also involves the effects of 



the e lectron-phonon interaction. Unscreened EPI makes polarons very mobile (lAlexandrov 



screening on 



19961 ). which leads to strong effects even on the qualitative physical properties of the polaron 



gas. Modeling the screening ef f ects, the interaction force between electrons and phonons was 



introduced in (ISpencer et al. 



20051 ) in the form of a screened discrete Frohlich interaction. 



It describes EPI of holes with c-axis polarised lattice distortions, which has been suggested 
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as the relevant electron-phonon interaction in the cuprates (lAlexandrov I . Il996l ). 

The CTQMC polaron mass for the one-dimensional lattice with the screened Frohlich 
EPI, Eg. (12021) is shown in Fig{T2] at four different values of the screening length, Rg^ — > 
(the short-range Holstein interaction), Rsc = 1, -Rsc = 3, and Rsc oo (the non-screened 
Frohlich interaction). 

An important observation is evident from the plot of ln(m*/mo) against A shown in 
figure f|T2l) . At intermediate and large couplings (that is, in the transition and small polaron 
regions), altering the value of Rsc has a dramatic effect on the effective mass. For example, 
the non-screened Frohlich polaron is over 10^ times "lighter" than the Holstein polaron at 
A = 4, and over 10^ times lighter at A = 5. It is apparent from the above results for the 
screened Frohlich model that as -Rsc increases (from Holstein to Frohlich), the QMC results 
move, in general, closer to the "1/A" expansion over the entire range of A. That is, the 
"1 / A" expansion becomes generally more applicable as the range of EPI increases (as well 
as with increasing ui). 

The polaron features for the long-range EPI were also investigated by extending a varia- 
tional approach prev iously proposed for the study of systems with local (Holstein) coupling 



(jPerroni et all |200J). The ground state spectral weight, the average kinetic energy, the 
mean number of phonons, and the electron-lattice correlation function were calculated for 
a wide range of model parameters focusing on the adiabatic regime. A strong mixing of 
electronic and phononic degrees of freedom even for small values of EPI was found in the 
adiabatic case due to the long-range interaction. 

F. Isotope effect on the polaron mass and the polaron band-structure 

There is a qualitative difference between ordinary metals and polaronic (semi) conductors. 
The renormalised effective mass of electrons is independent of the ion mass M in ordinary 
metals (where the Migdal adiabatic approximation is believed to be valid, see below), because 
A does not depend on the isotope mass. However, when electrons form polarons dressed by 
lattice distortions, their effective mass m* depends on M through m* = mexp{'~fEp/ujQ), 
in the strong-coupling limit. Here the phonon frequency depends on the ion mass, so that 
there is a large polaronic isotope effect (PIE) on the ca rrier mass w i th the carrier mass 



isotope exponent am* = (1/2) ln(m*/m) as predicted in (lAlexandrov I . Il992bl ). in contrast 
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FIG. 12 The logarithm of the effective ma s s for s creening lengths Rsc = 0, 1, 3, oo (increasing size of 



circles) versus A at w = 1 (jSpencer et al. 



20051 ). At intermediate and strong couphng, decreasing 
the value of R^c dramatically increases the effective mass. The curves tend to the strong-coupling 
analytical result (dashed lines). 



to the zero isotope effect in ordinary metals (see (lAlexandrov 



effect was found ex perimentally in cu prates fjKhasanov et 



Zhao et al. 



19971 ) and manganites (lAlexandrov et al. I . l200ll ) . More recent high resolution 



al. 



2004 



2003) for more details). The 



Zhao and Morris 



1995 



angle resolved photoemission spectroscopy (ARPES) also revealed a complicated isotope 
effect on the whole ban d-structure in cupra te superconductors depending on the electron 
energy and momentum flGweon et al. 1 . 120041 ) . 



PIE in the intermediate re gion of parameters was ca 



field appro ximation (DM 
algorithm ( Hague et al. I . 



T) (jpeDDeler and Millis 



2006 



3 



2002 



c ulated using the dynamic mean 



Pad et al. 



Kornilovitch and Alexandrov 



2004 



2mm. and C 



Spencer et al. 



QMC 



20051) 



Importantly and the effective mass averaged over dimensions are related to the critical 
temperature isotope exponent, a = —dlnTc/d\nM, of a (bi)polaronic superconductor as 



a = ar 



m/m* 
X- He 



(203) 



where /Xc is the Coulomb pseudo-potential (lAlexandrov 



1992b 



20031 ). 
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FIG. 13 Isotope effect on polaron spectrum and density of states. Top row: the 2d Holstein 
polaron at A = 1.2. Left: polaron spectrum at a; = 0.80 and 0.75 (in units of t). Middle: the 
isotope exponent for each k-point. Right: the density of states for t he two frequencies. Bottom row . 



the s ame for the discreate 2d Frohlich polaron at A = 2.4 (after (jKornilovitch and Alexandrov 



200J)). 



In the adiabatic regime the isotope effect on the polaron mass does not fully represent the 
isotope effect on the vast majority of polaron states, in particular for the short-range EPI, so 
that additional insight can be gained from the isotope effect on the entire polaron spectrum 
(jKornilovitch and Alexandrov 1 . 120041 ). The isotope effect on polaron spectrum and DOS in 
d = 2 is illustrated in Fig [131 The ratio of the two phonon frequencies, u = 0.80 1 and 0.75 1, 
has been chosen to roughly correspond to the substitution of ^^O for ^^O in complex oxides. 
One can see that the polaron band shrinks significantly, by 20-30%, for both polaron types. 
The middle panels show the isotope exponents on spectrum points calculated as 

1 {uj) Aek 



(204) 



2 (ek) Aa; ' 

where the angular brackets denote the mean value of either the two frequencies or of the 
two energy values. An interesting observation is that Ok of the Frohlich polaron is roughly 
independent of k (±10%). In the Holstein case ak dips in the vicinity of the F point. 

Unconventiona l isotop e effects as observed in high-temperature superconducting cuprates 



( iKhasanov et al. I . |200J), were also explained by polaron formation stemming from the 
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coup ling; to the particular quadrupolar Q(2)-type phonon mode (IBussmann-Holder et al. 



20051 ). Polaronic isotop e effects in the spectral fu nction of strongly correlated systems were 



numerically stu died by ( 



Fratini and Ciuchi 



model and by (IMishchenko and Nagaosa 
EPI. 



2005f) in the framework of the Holstein- Hubbard 



20061 ) using the extended t — J model including 



G. Jahn- Teller polaron 



The density-displacement EPIs discussed above are not only possible types of t he electron- 



lattic e coupling. Examples of other types are the Su-Schrieffer- Heege r (SSH ) EPI (ISu et al. 
19791 ) ( for recent path-integral results fo r SSH polarons see jzoh], |2Qo3)), and the Jahn- 



Teller (JT) EPI (jjahn and Teller I . 119371 ). In the former, the lattice deformation is cou- 
pled to electron kinetic energy, while JT involves a multidimensional electron basis and 
a multidimensional representation of the deformation group. The JT interaction is active 
in some molecules and crystals of high point symr aetry, and it has served as a guide in 
the search for high-temperature superconductivity (IBednorz and Miiller I . Il988l ). Later on 
JT-type EPIs were widely discussed in connection with cuprate and other high-T^ supercon- 



ductors (see, for example, ( Alexandrov et 



Mihailovic and Kabanov 



see, for example, R efs . 



1995 



Roder et al. 



2001 



Miiller 



al. 



1996 



Gunnarsson 



1997 



Mertelj et al. 



2005 : 



2000)) and co 



Bussmann-Holde 



199 



: and Bishop 



Yunoki et al. 



1998 



3)- 



ossal-magnetoresistance manganites 



1999 



Edwards 



2002 



Millis et al. 



The simplest model of the JT interaction is the E ® e interaction (iKanamori I . Il960l ) 
that describes a short-range coupling between twice-degenerate electronic levels (ci,C2) 
and a local double-d egenerate vibron mode {Crj). The Hamiltonian reads ( here we follow 



(IKornilovitch 



20071 )) 



HjT= - tY^ (cn'lCnl + cl,2Cn2^ 
(nn'> 

1 / 



K 



+ —^{Cl + r,l) 



(205) 



2M V^Cn dr]lj ' 2 
The symmetry of the interaction ensures the same coupling parameter n for the two phonons. 
Because the ionic coordinates of different cells are not coupled, the model describes a coUec- 
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tion of separate clusters that are linked only by electron hopping. To relate the Hamiltonian 



to niore realistic situations, p 



1999 



lonon dispersion should be added ( lAllen and Perebeinos 



Mihailovic and Kabanov 



200l|). 



An important property of the E ® e interaction is the absence of an exact analytical 
solution in the atomic limit t = 0. Here, in contrast with the Holstein and Frohlich EPI, 
the atomic limit is described by two coupled partial differential equations for the electron 
doublet '4)ifi{C,il)- At large couplings, however, the elastic energy assumes the Mexican 
hat shape and the phonon dynamics separates into radial oscillatory motion and azimuthal 
rotary motion. This results in an additional pre-exponential factor oc k in the ion overlap 
i ntegral, l eadin g to the effective mass mjrp/m = (2/71(7)^/^ expgf^, where g"^ = k? /2Muj'^ 



f lTakadal . l2000h . 



A path integral appr oach to Hami 



Its details are found in ( Kornilovitch 



tonian fl205p was developed by Kornilovitch (2000). 



20071 ). Because there are two electron orbitals, the 



electron path must be assigned an additional orbital index (or colour) a = 1,2. Colour 1 
(or 2) of a given path segment means that it resides in the first (second) atomic orbital of 
the electron doublet. There is a difference between the two phonons. Phonon ( is coupled 
to electron density, like in the Holstein case. The difference from the Holstein is that 
the direction of the force changes to the opposite when the electron changes orbitals. In 
contrast, phonon rj is coupled to orbital changes themselves: the more often the electron 
changes orbitals, the more "active" is 77. Discrete orbital changes are analogous to electron 
hops between discrete lattice sites, and as such are associated with "kinetic orbital energy" . 
Phonon couplin g to orbita l chan ges is analogous to phonon coupling to electron hopping in 



the SSH model (ISu et al. 



19791). 



The JT polaron properties can be calculated with no approximations using QMC algo- 
rithm. The effective mass, spectrum, and density of states are obtained as for the conven- 
tional lattice polarons. Most propertie s behave similarly t o those of the d = 3 Holstein 



polaron at the same phonon frequency (iKornilovitch 1 . 120001 ). For example, the kinetic en- 



ergy sharply decreases by absolute value between A = 1.2 and 1.4. The JT polaron mass is 
slightly larger at the small to intermediate coupling, but several times smaller at the strong 
coupling. This non-monotonic behaviour of the ratio of the JT and Holste in masses was later 



confirmed by accurate variational calculations (lEl Shawish et al. 



20031 ). ahhough in that 



work the JT polaron (and bipolaron) was investigated in one spatial dimension. The relative 
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lightn ess of the JT polaron is consistent with Takada's result mentioned above (ITakada I . 
2000|). Finally, the density of JT polaron states features the same peak at the top of the 
band, caused by the spectrum flattening at large polaron momenta, see Fig. (IT5]) . 

In summary, the local character of the JT interaction and the independence of vibrating 
clusters result in the same extremity of polaron properties as in the 3d Holstein model. One 
could expect that either a long-range JT interaction or ph onon dispersion will soften the 
sharp polaron features and make JT polarons more mobile (IKornilovitch 1 . 120071 ). 



V. RESPONSE OF DISCRETE POLARONS 



A. Hopping mobility 

Studies o : 



f lloffe 



lattice polarons allowed for a theoretical understanding of low mobility problem 
19571 ) of many "poor" conductors, where an estimate of the mean-free path yields 
values much lower than the lattice constant. Transport properties of lattice polarons depend 
strongly on temperature. For T lower than the characteristic phonon ener gy polaroii kinetic s 
is the Boltzmann kinetics of heavy particles tunnelling in the narrow band (IHolstein I . ll959bl ). 
However at higher temperatures the polaron band collapses fllV.B.2p, and the transport is 



1975; 


Holstein . 


1959bl: 


Khr 


leer . 


Yamashita and Kurosawa . 


1958 



iumps of polarons from site to site 



1979 



Lang and Firsov 



1963 



1968 



flAppel 



1968 



Firsov 



Mott and Davis 



1982 



)• 



The 1/ A perturbation expansion is particularly instrumental in calculating hopping mo- 
bilities and optical absorption coefficients. Applying the canonical transformation fllV.Bp 
and singling out the diagonal coherent tunnelling in the polaron band one can write the 
Hamiltonian as 



H = Hr, + H, 



ph 



H., 



inti 



(206) 



where Hp = ^k^kcj^Ck is the "free" polaron contribution (here we drop the spin), Hph = 
^qCi;q((i]j(iq + l/2) is the phonon part, and = Z'Ei^ — fi is the renormalised (polaron) band 
dispersion. The chemical potential fi includes the polaron level shift —Ep, and it could also 
include all higher orders in 1/A corrections to the polaron spectrum, which are independent 
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of k, Eq. fll75p . The band-narrowing factor Z' is defined as 

'Ylim ^(m)e^^^*-"'^ exp(— zk ■ m) 



Z' 



(207) 



Em^(m)exp(-«k- m) 

which is Z' = expl—'-fEp/uj)) where 7 < 1 depends on the range of EPI and phonon frequency 
dispersions ( IIV.CI) . The interaction term Hint comprises the polaron-polaron interaction, 
Eq. fll7ip . and the residual polaron-phonon interaction 



p—ph 



El 



(208) 



where (aij) means averaging with respect to the bare phonon distribution. In the framework 
of the single-polaron problem one neglects the polaron-polaron interaction and treats Hp_ph 
as a perturbation. 

The motion of the s mall polaron at h igh temperatures is a random walk consisting of 



steps from site to site (IHolstein 



1959bl ). Holstein calculated the probability W for the 



hop of a small polaron to a neighbouring site. He suggested that the random walk was 
a Markovian process. For such processes the diffusion coefficient is given by D ~ a^W, 
where W is the hopping probability. The only term in the polaronic Hamiltonian, which 
changes the phonon occupation numbers, is the polaron-phonon interaction Hp^ph- The 
nearest neighbour hopping probability in the second order is 



w = 2n{j2\mp-pkm's (y1 

j \ q 



^qK -<) ), 



(209) 



where \i) and |j) are the eigenstates of Hq corresponding to the polaron on site i with 
phonons in each phonon mode and the polaron on the neighbouring site j with phonons, 
respectively. Replacing the 6 function in Eq. fl209p by the integral yields 



£*[exp(li:Wq)P|l--«l-a)l^)-l 



-iO+|T| 



(210) 



The integration over r is performed using a saddle-point approximation by allowing a finite 
phonon frequency dispersion Sou <ti u^. Expanding cos(ci;qr) in powers of r, and using the 
Einstein relation fi = eD/T one obtains the hopping mobility. 



l^h = ea 



2T(E„T)i/2 



-Ea/T 



(211) 
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where 



T 



= iV l7(q)r[l - cos(q • a)] tanhK/4T) 



(212) 



is the activation energy, which is half of the polaron level shift Ep for the Holstein EPI. This 
expression is applicable at T ^ 0^0/2 and /ujQ^EaT) <^ 1. Importantly, increasing the EPI 
range diminishes the value of Ea further. The hopping mobility = crh/ne ~ exp{—Ea/T) 
can be below ea^ ~ Icm^/Vs, which is the lowest limit for the Boltzmann theory to be 
applied. Within the Boltzmann theory such a low mobility corresponds to the mean free 
path / < a, which is not a reasonable result. 

Holstein suggested that at low temperatures T < 0^0/2 In 70, there is the ordinary Boltz- 
mann transport in momentum space, but in the narrow polaron band. According to 



(IHolstein 



1959bl ). the transition from the band regime to the hopping regime occurs when 



the uncertainty in the polaron energy becomes comparable to the width of the polaron 
band. If phonons dominate in the scattering, the polaron mobility, fi = fit + fJ'h decreases 
when the temperature increases from zero to some Tmin, where it is at minimum, because 
its Boltzmann part, fit, falls down due to an increasing number of phonons, while the hop- 
ping part, fih, remains small. However fi increases above Tmin, due to the thermal activated 
hopping. There is a wide temperature range around Tmin where the thermal activated hop- 
ping still makes a small cont ribution to the cond uctivity, but the uncertainty in the polaron 



band is already significant fiFirsov 



special diagrammatic technique (iFirsov 



19751 . 



2007 



1975 



Th e polaron transport theory req uires a 



2007 



Lang and Firsov 



1965 



1963 



19681 ) and a 



19661 ) in this re- 



conditional-probability function description (IKudinov and Firsov 
gion. 

The transverse conductivity a^y and the Hall coef ficient Rh = (J xy/ Hal^ of lattice po- 



larons can be calculated with the Peierls substitution (jPeierls 



19331 ) 



tim — n) 



t(m — n)e' 



-jeA(m)-(m— n) 



(213) 



which is a fair approximation if the magnetic field, B, is weak compared with the atomic field 
eBa? « 1. Here A(r) is the vector potential, which can be also time dependent. Within the 
Boltzmann theory the sign of the Hall coefficient Rh — ±l/en depends on the type of carriers 
(holes or electrons) and the Hall mobility fin = Rh<^xx is the same as the drift mobility fit up 
to a numerical factor of the order of 1. The calculat i ons o f the hopping Hall cur r ent similar to 



those of the hopping conductivity (IBottger et al. 



2007 



Bryksin and Firsov 



1970 



Emin 



81 



1971 



Friedman and Holstein 



19631 ) show that the Hall mobility depends on the symmetry 



of the crystal lattice and has nothing in common with the hopping mobility, neither with 
respect to the temperature dependence and even nor with respect to the s ign. In particular 
for hexagonal lattices three-site hops yield (IFriedman and Holstein I . Il963l ) 



TT 



ea 



-Ea/3T 



(214) 



:i2EaT)V2 

with the same sign for electrons and holes. The activation energy of the Hall mobility is three 
times less than that of the hopping mobility. In cubic crystals, the hopping Hall effect is 



gover ned by four-site hops. The four-site calculations in (IBryksin and Firsov 



1970 



Emin 



197ll ) gave the Hall mobility with the "normal" sign depending on the type of carriers. 



The Hall conductivity and resistivity of strongly localized electrons at low temper- 
atures and small magnetic f ields strongly depend on frequency, the size of a sample 



( lEntin-Wohlman et al. 



19951) ■ and on a magnetic order like, for example in ferromagnetic 



(Ga,Mn)As (IBurkov and Balents I . l2003l ). In the presence of the spin-orbit interaction each 
hopping path acquires a spin-dependent phase factor of the same form as that in a perpen- 
dicular (to the 2D system) magnetic field, which leads to spin accumulation and spin-Hall 



effects ( iBottger et al. 



2007 



Entin-Wohlman et al. 



20051). 



B. Optical conductivity 



One of the fingerprints of lattice polarons is the frequency {uj) and temperature de- 
pendence of their mid-infrared (MIR) conductivity a{uj). In the low frequency a nd low 
temperature region, where the tunneling band transport operates (IHolstein I . Il959al jbl). the 
conductivity acquires the conventional Drude form. 



(215) 



l + (c^r)2' 

where the transport relaxation time r may be frequency dependent because of the narrow 
band. For high (mid-infrared) frequencies, well above the polaron bandwidth {uj ^ w) but 
below the interband gap, the Drude law is not applied. In this frequency region one can 
apply the generalised Einstein relation o-{ijj) = eD{uj)/uj, where D{uj) = a '^Wju), and W(uj ) 
is the hopping probability of the absorption of the energy quantum u (IHolstein I . Il959bl ). 
The number of nearest neighbour transitions per second with the absorption of a photon of 
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the energy u is given by the Fermi-Dirac golden rule, 



(216) 



and with the emission, 



As a result one obtains 



2irC£\m,^vK\3)?5[Y, 

3 



(217) 



W{uj) = W- -W^ = 2t^e-^3' sinh(cj/2T) 



dre 



cos cur 



• (218) 



As in the case of the hopping mo bility discussed above t 



the saddle-point approximation (IBottger and Bryksin 



sinh u;q/2T 
le int egral over r is calculated using 



19851 ) 



a{uj) = ah 



Tsinh(cj/2T) 



where 



1 a;^ 

<^ = y l7(q)r[l - cos(q ■ a)]— —-3- -. 

2ArZ^i'^^^i^ ^^sinh(u;q/2T) 



(219) 



(220) 



For high temperatures, T >> ujq/2, this expression simplifies as 



2u{EaTyn 



exp 



4^a) 



(221) 



Here an = nefih and n is the atomic density of polarons. The frequency dependence of 



the MIR conductivity has a form o: 



with the half-width Ay/EaT (lEagles 



an asymmetric Gaussian p e ak ce ntered at u; = 4E, 



1963 



Klinger 



1963 



Reik 



19631 ). According to the 



Franck-Condon principle, the position of the ions is not changed during an optical transition. 
Therefore the frequency dependence of the MIR conductivity can be understood in terms 
of transitions between the adiabatic levels of the two-site Holstein model flIV.A.2p . The 
polaron, say, in the left potential well absorbs a photon through the vertical transition 
to the right well, where the deformation is lacking, without any change in the molecular 
configuration. The photon energy required to excite the polaron from the bottom of the well 
is ^ 2Ep, which corresponds to the maximum of o"(ci;) in Eq. fl22ip . The main contribution 
to the absorption comes from the states near the bottom with the energy of the order of 
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T. The corresponding photon energies are found in the interval 2Ep ± ^^^>EpT in agreement 
with Eg. (12211) . For low temperatures T < a;n/2 the h alf-width of the MIR maximum is 
about ~ ^J EpUjQ rather than ~ \/ EpT flEmin I . Il993l ). The optical absorption of small 
polarons is distinguished from t hat of large polarons fjllLFD by the shape and the temperature 



dependence. Their comparison (jCalvani 



2001 



Emin 



1973 



19931 ) shows a more asymmetric 



and less temperature dependent mid-infrared (MIR) absorption of large polarons compared 
with that of small polarons. The high-frequ ency behavi or of the optical absorption of small 



polarons is described by a Gaussia n decay ( 



Reik 



19721 ). Eg. (12211) . while for large polarons 



it is much slower power law dPevreese I . ll972l ). 

Many materials are characterised by intermediate values of EPI, A ~ 1, which requires 
an extension of the theory of optical absorption to the crossover region from continuous 
to lattice polarons. The intermediate coupling and frequency regime has been inaccessible 
for an analytical or semi-analytical analysis, with an exception of i n finite spatial dimen- 



sions , where DMFT yie 



1997 



Fratini and Ciuchi 



as re. 



2003 



iable results jBenedetti and Zeyher 



1998 



Ciuchi et al. 



Sumi 



19741 ). DMFT treats the local dynamics exactly, but 
it cannot account for the spatial correlations being important in real finite-dimensional sys- 
tems. Nevertheless DMFT allows one to address the intermediate coupling and adiabaticity 
regimes in the absorption not covered by the standard small polaron theory, where qualita- 
tively new features arise. In particular, the optical absorption exhibits a reentrant behavior, 
switching from weak-coupling-like to polaronic-like upon increasing the temperature, and 
sharp peaks with a nonmonoton i c tern perature dependence emerge at characteristic phonon 
frequencies (IFratini and Ciuchi I . l2006l ) . 

Earlier ED studies of the lattice polaro n absorption were 



ID and 2D clusters in the Holstein mode 



Fehske et al. 



1997 



Wellein and Fehske 



flAlexandrov et al. 



imited to small 2 to 10-site 



1994b 



Capone et al. 



1997 



19981 ). The optical absorption occurs as energy is 



transferred between the electromagnetic field and the phonons via the charge carriers. The 
vibration energy must be capable of being dissipated. Hence, using ED one has to introduce 
some contin uous density of phonon state s, or a phonon lifetime, which makes MIR absorption 
to be finite (lAlexandrov et al. I . Il994bl ). As a result one obtains a fair agreement between 
ED absorption spectra and the analytical results, Eq. (l22ip . in the strong coupling limit as 
far as a smooth part of uo dependence is concerned. The MIR conductivity occurs much 
more asymmetric in the intermediate coupling region than in the strong coupling regime, 
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and it shows an additional oscillating superstructure corresponding to a different spectral 
weight of the states with a different number of virtual phonons in the polaron cloud. 

More recent ED, VE D and a kernel polynomial method (KPM) (for a review of KPM 
see ( iWeifie et al. 1 . 12006| )) allowed for numerical calculations of lattice polaro n properties ir i 



the Holstein mode^ 



El Shawish et al. 



in the whole parameter ra nge on fairly large systems (IBarisic 



2004 



2003 



Schubert et al. 



20051 ) . The main sig nature of optical conductivity 



for intermediate-to-strong EP coupling in the adiabatic regime (IFehske and Trugman 



20071) 



is that the spectrum is strongly asymmetric, which is a lso characteristic for rath er large po- 
larons ( IIII.FI) . as observed in cuprate su perconductors (jMihailovic et al. . 1990 ). perovskite 
tungsten bronzes (IRuscher et al. 1 . 12008| ) and many other doped insulators. Importantly, the 
weaker decay at the high-energy side mee ts the experimental findings for many polaronic 



materials like Ti02 (IKudinov et al. I . Il969l ) even better than standard small-polaron theory. 



At larger EPI a more pronounced and symmetric maximum appears in the low-temperature 
optical response. When the phonon frequency becomes comparable to the electron trans- 
fer amplitude different absorption bands appear, whi ch can be classified accordin g to the 
number of phonons involved in the optical transition (IFehske and Trugman 1 . 120071 ). 



C. Spectral function of discrete strong-coupling polarons 



The polaron problem has the exact solution. Eg. (11721) . in the extreme infinite-coupling 
limit, A = oo, for any EPI conserving the on-site occupation numbers of electrons. For a 
finite but strong coupling 1/A perturbation expansion is applied. Importantly, the analytical 
perturbation theory becomes practically exact in a wider range of the adiabatic parameter 
and of the coupling constant for the long-range Frohlich interaction fllV.E.Sp . 

K eeping this in mind, let us calcula t e the polaron spectral function in the first order in 



1/A (lAlexandrov and Ranninger 



1992 



Alexandrov and Sricheewin 



2000l ) . We can neglect 



Hp_ph in the first-order if 1/A ^ 1. To understand spectral properties of a single polaron 



we also neglect the polaron-polaron interaction. Then the energy levels are 



Era = X^^k'^k + ^LU^[n^ + 1/2], 



(222) 



where .^k is the small polaron band dispersion with respect to the chemical potential, and 
the transformed eigenstates \rh) are sorted by the polaron Bloch-state occupation numbers, 
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rik = 0, 1, and the phonon occupation numbers, riq = 0, 1, 2, oo. 

The spectral function of any system described by quantum numbers m, n with the eigen- 



values Em is defined as 



A(k, uj) = 7r(l + e-"/^)e^/^ ^ e"^"/^ | {n\ Ck \m)f + cu). (223) 

It is real and positive, yl(k, cj) > 0, and obeys the sum rule 

1 

- / dujA(k,uj) = 1. (224) 
The matrix elements of the electron operators can be written as 

{n\ Ck |m) = ^ ^ e"*''-'" (n| qX^ |m) (225) 



111 

using the Wannier representation and the Lang-Firsov transformation. Now, applying the 
Fourier transform of the 5-function, the spectral function is expressed as 

1 f°° 1 

A{k,uj) = - rfte*-*^ e'k.(n-m) ^ ^226) 



■((c,(t)X,(t)cJx;)) + ((cjX]c,(t)X,(t))) 

Here the quantum and statistical averages are performed for independent polarons and 
phonons, therefore {{ci{t)Xi{t)Xjcl)) = {{ci{t)Cj)) {{Xi{t)Xj)) . The Heisenberg free-polaron 
operator evolves with time as Ck(t) = Ckexp(— i^k^)? so that 

((Q(t)cl)) = ^ E = (227) 

k',k" 

— ^[1 - n(k')]e'^'-('"-")-'«'''*, 
k' 

{{4c^m = ^5^n(k')e^^'-(-'^)-^«-'* (228) 
k' 

where n(k) = [1 + exp,^k/^]^^ is the Fermi-Dirac distribution function of polarons. The 
Heisenberg free-phonon operator evolves in a similar way, (iq(t) = dqexp{—iujqt), so that 

((X,(t)Xj)) = J]((exp[Mi(q,t)c/q - H.c] exp[-uj{q)d^ - H.c.])), (229) 
q 

where Uij{q,t) = Mjj(q) exp(— icjqt). This average is calculated using the operator identity, 
as in Eq. ffTTTj) . 

{{Mt)X])) = exp 1-^ E l7(q)l'A(m - n,t)| , (230) 
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where 



/q(m,t) = [1 — cos(q ■ m) cos(a;qt)] coth ^ + i cos(q ■ m) sin(ci;qt) 



(231) 



Here we used the symmetry of 7(— q) = 7(q), so that terms containing sin(q ■ m) have 
disappeared. 

To proceed with the analytical results we consider low temperatures, T ^ cjq, when 
coth(c<jq/2T) ^ 1. Expanding the exponent in Ea.fl230l) and performing summat ion over 
m, n, k' and integration over time we arrive at (lAlexandrov and Sricheewin 1 . 12000| ) 



A(k,c.) = 5^U-)(k,a;) + 4+)(k, 



1=0 



where 



qi,...q; 



n!-=i l7(qr 
{2Nyi\ 



(232) 



(233) 



1 - k- ^ q. 



6 CJ-^CJq^-e 



r=l 



(234) 



qiv-q; 



n 



k+ ^ qj 5 L+ 5^ ^q.-^k+EUi q. ' 

^ r=l / \ r=l / 



andZ = exp(-(2iV)-iEql7(q)P). 

Clearly Eg. (12321) is in the form of the perturbative multi-phonon expansion. Each con- 
tribution ylp''(k, cj) to the spectral function describes the transition from the initial state 
k of the polaron band to the final state k± Yll=i with the emission (or absorption) of / 
phonons. 



The 1/A expansion result, Eq. 



is different from the conventiona 



of met a' 



Migdal 



l ie elec trons coupled to phonons in the Migdal-Eliashberg theory (lEliashberg 



spectral fu r iction 



1960 



19581 ). There is no imaginary part of the self-energy since the exponentially small 



(at low temperatures) polaronic damping, Eq. fll87p . is neglected. Instead EPI leads to the 
coherent dressing of electrons by phonons, and phonon "side-bands" . The spectral function 
of the polaronic carriers comprises two different parts. The first (/ = 0) k-dependent coherent 
term arises from the polaron band tunnelling. 



Alr\k,u;) + Al^'(k,uj 



(235) 
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The spectral weight of the coherent part is suppressed as Z ^ 1. However in the case of the 
Frohhch interaction the effective mass is less enhanced , since the band narrowing factor Z ^ 
Eq. fl207p . in = Z'Ei^ — /i is large compared with Z (lAlexandrov and Sricheewin 1 . 12000| ). 

The second incoherent part AincohO^,^) comprises all the terms with / > 1. It describes 
excitations accompanied by emission and absorption of phonons. We notice that its spectral 
density spreads over a wide energy range of about twice the polaron level shift Ep, which 
might be larger than the unrenormalised bandwidth 2D in the rigid lattice without phonons. 
On the contrary, the coherent part shows a dispersion only in the energy window of the order 
of the polaron bandwidth, w = Z'D. It is interesting that there is some k dependence of 
the incoherent background as well, if the EPI matrix element and/or phonon frequencies 



depend on q. Only in the Holstein model with the short-range dispersio nless e-p 



1 interaction 



7(q) =7o and tUq = ujq the incoherent part is momentum independent (IMahan 



1990) 



Anco/^(k,^) = vr-J^^x (236) 
1=1 

5^ {[1 - n (k')] S (u-luo-^^,) + n (k') 6 (u+luo-^^,)} . 



As soon as we know the spectral function, different Green's functions (GF) are readily 
obtained using their analytical properties. For example, the temperature GF is given by the 
integral 

,A{k,u;') 



where Uk = TiT(2k + 1), = 0, ±1 



model (lAlexandrov and Ranninger 



dJ ^ ' ' . (237) 
2, ... . Calculating the integral we find in the Holstein 



1992|) 



lU)r. 



N 2^ 2H\ ^ 

1=1 k' 



n(k') 



n(k') 



(238) 



Here the first term describes the coherent tunnelling in the narrow polaron band while the 
second k-independent sum describes the phonon-side bands. 



D. Spectral function of discrete all-coupling polarons 



The spectral function, Eq. 



satisfies the major sum rule, Eq. 



However the 



higher- momentum int egrals, /"^^ dujuj^ACk , u) with p > 0, calculated using Eq. (12381) . differ 



from the exact values (IKornilovitch 



2002) by an amount proportional to 1/A. The difference 
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is due to a partial "undressing" of high-energy excitations in the side-bands, which is beyond 
the first order 1/A expansion. A rather accurate Green's function of the Holstein polaron has 



been recently obtained by sum ming all t 



its free propagators' momenta flBerciu 



le diagr ams, but with each diagram averaged over 



2006a 



b|). The resulting Green's function satisfies 



exactly the first six spectral weight sum rules. All higher sum rules are satisfied with great 
accuracy, becoming asymptotically exact for coupling both much larger and much smaller 
than the free particle bandwidth. 

The spectral properties of the Hols tein model in a wider parameter range have 



been studied numerical 



El Shawish et al. 



2003 



Fehske et al. 



y using ED (see (jAlexandrov et al. 



1997 



1994b 



Fehske and Trugman 



Capone et al. 



2007 



Loosed al. 



1997 



2006 



Ranninger and Thibblin 



1992 



Schubert et al. 



2005 



Wellein and Fehske 



19981), and ref- 



ID Holstein model is near 



erences therein). At the weak EPI, the electronic spectrum o: 
unaffected for energies below the phonon emission threshold (IFehske and Trugman I . l2007l ). 
Hence, the renormalised dispersion nearly coincides with the tight-binding cosine band 
(shifted oc Ep) up to some kx-, where the phonon energy intersects the bare electron band 
with the familiar flattening effect ( 1IV.F|) . Reaching the intermediate EPI (polaron crossover) 
regime a coherent band separates from the rest of the spectrum [kx ^ vr] . At the same time 
its spectral weight becomes smaller and will be transferred to the incoherent part, where 
several sub-bands emerge. In the strong-coupling case the coherent quasi-particle absorption 
band becomes extremely narrow and its bandwidth approaches the strong-coupling result. 
The incoherent part of the spectrum carries most of the spectral weight consisting of a 
sequence of sub-bands separated in energy by cjq; in agreement with the analytical results. 
Effects of the finite-r ange EPI on the spe ctral properties of lattice polarons have been 



studied numerically by flFehske et al 



in the framework of the ID model of (Alexandrov and Kornilovitch 



20021) using exact Lanczos di a gona 



isation method 



19991 ). The polaron 



band-structure has been calculated in agreement with the analytical and CTQMC results 
IIV.E.3[ The optical absorption of lattice polarons with a finite-range (Frohlich-type) EPI 
has been found similar to the continuous-polaron absorption ( lIILFp for all EPI strengths. 
The polaron features due to EPIs with dif ferent coupling rang es were also investigated in 
the framework of the variational approach (jPerroni et alV |2005| ). The ground-state energy, 
the spectral weight, the average kinetic energy, the mean number of phonons, and the 
electron-lattice correlation function were calculated for the system with coupling to local 
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and nearest-neighbor lattice displa cements and compared with the long-range case. As in 
(lAlexandrov and Sricheewin 1 . 12000| ) a substantially different mass renormalisation compared 
with the coherent weight reduction, Z <^ Z\ was found for the finite-range EPI. 



VI. BIPOLARON 

A. Polaron-polaron interaction 

Polarons interact with each other, cf. Eg. (11711) for small polarons. The range of the 
deformation surrounding the Frohlich polarons is quite large, and their deformation fields are 
overlapped at finite density. Taking into account both the long-range attraction of polarons 
owing to the lattice deformations and their direct Coulomb repulsion , the residual long-range 



inter action turns out rather weak and repulsive in ionic crystals (lAlexandrov and Mott 



19951 ). In the long- wave limit (g ^ vr/a) the Frohlich EPI dominates in the attractive part. 



but polarons repel each other at large distances, |m — n| ^ a, 

t;(m-n) = — > 0. (239) 

eolm — n| 

The Frohlich EPI nearly nullifies the bare Coulomb repulsion, if eo ^ l? but cannot over- 
screen it at large distances. 

Considering the polaron-phonon interaction in the multi-polaron sys tem we have to tak e 
into account the dynamic properties of the polaron response function (lAlexandrov 1 . 120001 ). 
One may erroneously believe that the long-range Frohlich EPI becomes a short-range (Hol- 
stein) one due to screening of ionic potentials by heavy polaronic carriers. In fact, as distinct 
from large polarons, small polarons cannot screen high-frequency optical vibrations because 
their renormalised plasma frequency i s comparable with o r even less than the phonon fre- 



quency in the strong-coupling regime (lAlexandrov I . Il992al ) . The small-polaron plasma fre- 



quency is rather low due to the large static dielectric constant, e^'^ X^'vn ionic lattices and 
the enhanced polaron mass vn* ^ m. If cjq > the singular behaviour of the Frohlich 
EPI, 7(q) ~ 1/g, is unaffected by screening. Polarons are too slow to screen high-frequency 
crystal field oscillations. As a result, EPI with high-frequency optical phonons in ionic solids 
remains unscreened at any density of polarons. 

Another important issue is a possibility of the Wigner crystallization of polarons. Because 
the net long-range interaction is relatively weak, a relevant dimensionless parameter = 
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m*ae /eo(47rn/3)i/3 is not very large in ionic semiconductors. The Wigner crystallization 
appears around ^ 100 or larger, which corresponds to the atomic density of polarons 
n < 10"^ with = 30 and m* = 5m. This es t imate tells us that sufficiently mobile small 



polarons can be in a liquid state (lAlexandrov 



20001) at substantial doping levels, however 



they can be crystallised at low doping (see IIILH|) . 

When the short-range deformation and molecular-type (i.e. Holstein) EPIs are added 
to the Frohlich interaction, two polarons attract each other at a short distance of about 



the lattice constant. Then, owing to a nar row band, two lattice polarons easi 



tightly bounc 


state. 


i.e. a small t 


ipolaron ( 


Alexandrov and Rannineer 


1981a 


b; 


Anderson , 


1975; 


Aubrv . 


1993; 


Beni et al. . 


1974; 


Street and Mott . 


1975 


). One can estimate the cou- 



y form a local 



pling c onstant A and the a diabatic ratio Wo/t, at which the small bipolaron "instability" 
occurs (lAlexandrov I . l2000l ). The characteristic attractive potential is |f| = D/[\ — fic), 
where /ic is the dimensionless Coulomb repulsion, and A includes the interaction with all 
phonon modes. The radius of the potential is about a. In three dimensions a bound state 
of two attractive particles appears, if |f | > 7r^/8m*a^. Substituting the polaron mass, 



m 



[2aH] ^ exp{'y\D/ujQ), we find 



— < ilzX)-' In 



IT 



(240) 



[42;(A-/iJJ 

As a result, small bipolarons form at A > /Xc + tt^/4:Z, which is almost independent of the 
adiabatic ratio. 

The formation of sm all bipolarons is closely related to a negative U -c enter model 



Street and Mott 



dLegget 



3) 



19751 ). Starting from initial works by Eagles f lEagles 1 . 1 19691 ) and Legget 



19801) this model receivec 



superconductivity (IMicnas et al. 



part icular attention in connection with high-temperature 



1990). However, in using the negative Hubbard U model, 



one has to realize t hat this model, which predi c ts a s mooth crossover from Cooper pairs 



to real-space pairs (iNozieres and Schmitt-Rinkl . Il985l ). cannot account for the polaron- 
bipolaron crossover with the increasing EPI. The essential effect of the polaron band- 
narrowing, which is responsib l e for high critical temperatures in the model of polaronic 
superconductors (lAlexandrov I . Il983l ). is missing in the negative Hubbard U model. 
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B. Holstein bipolaron 



The attractive energy of two small polarons is generally much large r than the polaron 



1983 



bandwidth, which allows for a consistent treatment of small bipolarons ( Alexandrov 

irnn 

Alexandrov and Ranninger I . ll981al Jbl). Under this condition the hopping term in the trans 



formed Hamiltonian if is a small perturbation of the ground state of immobile bipolarons 
and free phonons, 

H = Ho + Hpert, (241) 

where 

Ho = IY. ^-^cJc.c, + ^ + 1/2] (242) 



and 



pert 



ijclcj 



Anderson 



1975 



Aubry 



1993 



1995 



Beni et al. 



(243) 



arons, which are the ground state 



Let us first discuss the dynamics of onsite hipo 
of the system w i th the Holste i n non- dispersive EPI flAlexandrov and Raiininger 



1974 



Street and Mott 



3- 



1981a 



19751). The on- 



site bipolaron is formed if 2Ep > U, where U is the onsite Coulomb correlation energy (the 
Hubbard U). The intersite polaron-polaron interaction, Eg. (11711) . is just the Coulomb repul- 
sion since the phonon mediated attraction between two polarons on different sites is zero in 
the Holstein model. Two or more onsite bipolarons as well as three or more polarons cannot 
occupy the same site because of the Pauli exclusion principle. Hence, bipolarons repel single 
polarons and each other. Their binding energy, A = 2Ep — U, is larger than the polaron 
half-bandwidth, A ^ w, so that there are no unbound polarons in the ground state. Hp^rt, 
Eq. fl243p . destroys bipolarons in the first order. Hence it has no diagonal matrix elements. 
Then th e bipolaron dynamics is described b y the use of another canonical transformation 
exp(5'2) flAlexandrov and Ranninger I . Il981al ). which eliminates the first order of Hp^rt, 

{f\^tjcUj\p) 



Ef - Ep 



(244) 



Here Ef^p and |/), \p) are the energy levels and the eigenstates of Hq. Neglecting the terms 
of the order higher than {w/A)'^ one obtains 



ff 



e^^ He- 



Si 



ff 



(245) 
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{Hh)ff' ^ (i/o)//'-^Xl Yl (f\^^i'4ci'\p){p\^jj'c]cj,\f') X 



2 



Ep-Ef Ep~Ef 

The bipolaronic Hamiltonian Hf, is defined in the subspace |/), |/') with no single (unbound) 
polarons. On the other hand, the intermediate bra {p\ and ket \p) refer to configurations 
involving two unpaired polarons and any number of phonons. Hence we have 

Ep-Ej = A + J2^^^ fc - ^i) , (246) 

ci,u 

where ra^'^ are the phonon occupation numbers (0, 1, 2, 3...oo). 

The lowest eigenstates of Hb are in the subspace, which has only doubly occupied 
^ms'^ms'lO) empty |0) sites. On-site bipolaron tunnelling is a two-step transition. It 
takes place via a single polaron tunneling to a neighbouring site. The subsequent tunnelling 
of its "partner" to the same site restores the initial energy state of the system. There are no 
real phonons emitted or absorbed because the (bi)polaron band is narrow. Hence we can 
average Hh with respect to phonons, 

Hh = Ho-t J2 5Z ~ m')t(n - n') x (247) 

m^m',s n^n',s' 

POO 

Jo 

Here $J^JJ^'(i) is a multiphonon correlator, 

^^Ut) = {{x}{t)Mt)x]Xr)), (248) 



which can be readily calculated as (lAlexandrov I . Il992al ) 



*mm'(*) = e-^'^'^-'^'^e-^?'^"""') X (249) 

i^.(m,m,n,n) --^^ 

where -Fq(m, m', n, n') = cos[q ■ (n' — m)] +cos[q ■ (n — m')] — cos[q ■ (n' — m')] — cos[q ■ (n — 
m)]. 

Taking into account that there are only bipolarons in the subspace, where Hh operates, 
one can rewrite the Hamiltonian in terms of the creation b]^ = c^j^jcj^j^ and annihilation 
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CmiCmt operators of singlet pairs (lAlexandrov and Ranninger I . Il981al ): 



m 

E 



1 



nm + 



(250) 



t,,(m - m')&Jn&m' + -jvim - va!)n^n^, 



Here TT-m = fej^fem is the bipolaron site-occupation operator, y(m — m') = 4t;(m — m') + 
f *^^)(m — m') is the bipolaron-bipolaron interaction including the direct polaron-polaron 
interaction w (m — m') and a repulsive correlation energy, 



2z / rfre--^$-^(r). 

'0 



(251) 



This additional repulsion appears since a virtual hop of one of two polarons of the pair 
is forbidden, if the neighbouring site is occupied by another pair. The bipolaron transfer 
integral, is of the second order in t(m) 



^^(m — m') = — 2zt^(m — m') 



dre-^^-<l>--;(r) 



(252) 



The multiphonon correlator is simplified for dispersionless phonons at T ^ cuq as 



-2g2(m-m') 



exp [2g 



- m')e-^^"* 



which yields (lAlexandrov and Kabanov I . Il986l ) 



t(ui) 



2t\m) 



1=0 



and 



(253) 



(254) 



A " ^ mi + liUo/A) 

When A ^ C(Jo, we can keep the first term only with / = in the bipolaron hopping inte- 
gral, Eq. fl252p . In this case the bipolaron half-bandwidth zt{a.) is of the order of 2w'^ / [zA). 



However, if the bipolaron binding energy is large, A ^ cuq, the bipolaron bandwidth dra- 
matically decreases proportionally to e~^^^ in the limit A ^ oo. This limit is not realistic 
since A = 2Ep — Vr < 2q'^ujn. In a more realistic regime, ujq < A < 2g'^LjQ, Eg. (12521) yields 
(lAlexandrov and Kabanov I . Il986l ) 

2^/2^1"^ (m) 



tb{m.) 



V^oA 



exp 



uo \ A 



(255) 
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On the contrary, the bipolaron-bipolaron repulsion, Eg. (12540 . has no small exponent in the 



limit A oo, ' oc D / A. Together with the direct Coulomb repulsion the second order 
v^"^^ ensures stability of bipolarons against clustering. 

Interestingly the high-temperature behavior of the bipolaron bandwidth is just opposite 
to that of the small polaron bandwidth. While the polaro n band collapses with i ncrea sing 
temperature (lIV.B.2p . the bipolaron band becomes wider (IBryksin and GoFtsev I . Il988l ). 



tb[m) oc —= exp 



T 



Ep + A' 
2T 



(256) 



for T > ujq. 

The hopping conductivity of strong-coupling on-site bipolarons in the Holstein-Hubbard 
model (HHM) was f ound small compared with t he hopping conductivity of thermally ex- 
cited single polarons (IBryksin and Gol'tsev I . ll988l ). However, as the frequency of the electric 
field increases, the dominant role in the optical absorption is gradually transferred to bipo- 
larons at low temperatures. Like in the single-polaron case the bipolaron optical absorption 
can be estimated using the Franck-Condon principle which states that optical transitions 
take place without any chang e in the nuclear configuration. The corresponding analysis 
(IBryksin and Voloshin I . Il984l ) shows that the absorption coefficient of light by the on-site 
bipolaron has three Gaussian peaks located at frequencies Q = AEa, 8Ea — U and IQEa- The 
lowest peak corresponds to the absorption by single thermally excited polarons. The highest 
peak is due to the shakeoff of phonons without dissociation of the bipolaron while the main 
central peak is the absorption involving dissociation. A generalisation of the Franck-Condon 
pri nciple for the optical absorption of int ersite bipolarons with a finite-range EPl was given 
by (lAlexandrov and Bratkovsky I . Il999bl ) . 

The optical absorption and single-particle spectral functions of the bipolaron in ID HHM 
in the wh ole range of pararneters w ere calculated using E D on a 2-site cluster wit h two 



electrons (Alexandrov et al. 



1994bl ). and more recently in (IHohenadler et al. 1 . 12005| ) using 



the cluster perturbation theory (CPT). The latter allowed one to calculate the spectrum 
at continuous wave vectors and to find pronounced deviations (e.g. noncosine dispersions) 
of the bipolaron band structure from a simple tight-binding band due to an important 
contribution from the next-nearest-neighbour hoppings. 



Treating phonons classi c ally in the extreme adiabatic limit (lAubry 



1993 



1995 



2007 



Proville and Aubry 



1998 



19991 ) found along with the onsite bipolaron (SO) also an 



95 



anisotropic pair of polarons lying on two neighboring sites, i.e. the inte rsite bipolaron 
{SI). Such bipolarons were originally hypothesized in flAlexandrovl . Il99ll ) to explain the 
anomalous nuclear magnetic relaxation (NMR) in cuprate superconductors. The intersite 
bipolaron could take a form of a "quadrisinglet" {QS) in 2D HHM, where the electron den- 
sity at the central site is 1 and "1/4" on the four nearest neigbouring sites. In a certain 
region of U, where QS is the ground state, the double-well potential barrier, which usually 
pins polarons and bipolarons to the lat tice, depresses to almost zero, so that adiabatic lattice 



bipolar ons ca n be rather mobile (see (lAubry 



20071 ) for more details). 



(jSill . Il999l ) investigated the stability of the bipolaronic phase in HHM using a modified 
Lang-Firsov variational transformation with on-site and nearest-neighbour lattice distor- 
tions. There is a critical on-site Hubbard repulsion Uc below which the bipolaronic phases 
are stable for a fixed electron-phonon coupling. In the absence of on-site repulsion, bipola- 
ronic phases are stable over the entire range of electron-phonon coupling for one dimension, 
whereas there is a critical electron-phonon coupling for formation of a stable bipolaron in 
two and three dimensions. Mobile 51 bipolarons were found in ID HHM using variational 
methods also i n the non- and near- adiabatic re gimes with dynamical quantum phonons 



( iBonca et al. 



2000 



La Magna and Pucci 



19971 ). The intersite bipolaron with a relatively 



small effective mass is stable in a wide region of the parameter s of HHM due to both ex- 



change and nonadiabaticity effects (ILa Magna and Pucci 



19971 ). Near the strong coupling 



limit the mobile 51 bipolaron has an effective mass of the order of a single Holstein po- 
laron mass, so that one should not rule out the pos sibility of a superco nducting state of 



5*1 bipolarons with s or d-wave symmetry in HH M 



Bonca et al. 



2000). The recent di- 



agrammatic Monte Carlo study (IMacridin et al. I . |2004| ) found SI bipolarons for large U 
at intermediate and large EPI and established the phase diagr am of 2D HHM, compr ising 
large and small unbound polarons, 5*0 and 5*1 domains. Ref. (IMacridin et al. I . |200J) em- 
phasised that the transition to the bound state and the properties of the bipolaron in HHM 
are v ery different from bound sta tes in the attractive (negative U) Hubbard model without 
EPI ( iRobaszkiewicz et al. I . ll98ll ). 



The two-dimensiona 



many-body HHM was examin ed within a fluctuation-based effective 



cumulant approach by (IHakioglu and Zhuravlevl . Il998l ) confirming that the numerical results 



on systems with finite degrees of freedom fllV.Ep can be qualitatively extended to the systems 
with large degrees of freedom. When the electron-electron repulsion U is dominant, the 
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transition is to a Mott insulator; when EPI dominates, the transition is to a bipolaronic 
state. In the former case, the transition was found to be second order in contrast to the 



transition to the bipolaronic state, which is first order for larger values of U (IKoUer et al. 
2004bl ). 



For a very strong electron-phonon coupling polarons become self-trapped on a single 
lattice site and bipolarons are on-site singlets. A finite bipolaron mass appears only in 
the second order of polaron hopping, Eq. fl252p . so that on-site bipolarons might be very 
heavy in the Holstein model with the short-range EPI. As a result the model led some 
authors to the conclusion that the formation of itine r ant sniall po larons and bipolarons in 



real materials is unlikely (jPe Mello and Ranninger 



bipolaronic superconductivity is impossible (lAnderson 



1997, 



1993 



1998af). and high-ternperat ure 



Chakraverty et al. 



19981). 



C. Continuum Frohlich bipolaron 



While polarons repel each other at large distances, two large polarons can be also bound 
into a large bip olaron by an exchange interaction even without any additional EP I but 



the Frohlich one (Adamowski 



Kashirina et al. 



2005 



1989 



Bassani et al. 



Kochetov et al. 



mi 



1991 



Hiramoto and Tovozawa 



Mukhomorov 



1982 



1985 



Suprun and Moizhes 



1982 



Verbist et al. 



the reviews ( iDevreese 



1990 



1991 



19981: 



Vinetskii 



1961 



Vinetskii and Giterman 



Smondyrev and Fomin 



19571 ) (see also 



19941 )). Large bipolarons in the con- 



tinuum limit are referred to as Frohlich bipolarons. Besides the Frohlich coupling constant, 
a, the Frohlich bipolaron energy depends also on the dimensionless param eter U , a mea- 



sure f or th e strength of the Coulomb repulsion between the two electrons (IVerbist et al. 



1990 



199ll ). U = (e^/ecjo) v^™Jo- In the discussion of bipolarons the ratio r] = e / Eq of the 



high-frequency (electronic) and static dielectric constants is important (0 < ?7 < 1). The fol- 
lowing relation exists between U and a: U = Only values of U satisfying the inequality 
U > V2a have a physical meaning. It was shown that bipolaron formation is favoured by 
larger values of a and by smaller values of rj. 

Intuitive arguments suggesting that the Frohlich bipolaron is stabilised in going from 
3D to 2D had been given, but the first quant i tative analy sis based on the Feynman path 



integral was presented in Refs. fIVerbist et al 



1990 . 



199ll ). The conventional condition for 



bipolaron stability is E^p < 2Epo\, where Epoi and E^p denote the ground-state energies 



97 



15 



10 - 



c 
o 

JO 

u 

<D 5 - 



1 1 7Z 

Bipolaron Regiori^ 



■■ YBa^CUjO, 
U= 1.526 a 
LajCuOc 




a 



3D PHASE DIAGRAM . 



15 



10 

.2 

°- Ik 




2D PHASE DIAGRAM 

YBa,Cu30,- 

i7=1.526a 
La2Cu04- 

(7= 1.537a^!^'Bipolaron Region 



■ 1^ j^-.<*critical point a^saiZ.S 



X 



2 4 6 8 

Coupling constant a 



10 



FIG. 14 The stability region for bipolaron formation in 3D (a) and in 2D (b). The dotted line 
U = \pia separates the physical region ([/ > \/2a) from the non-physical [IJ < \f^a). The shaded 
area is the stability region in physical space. The dashed (dotted) "characteristic line" U = 1.537a 
{U = 1.526a) is determined hy U = \/ 2a/(l — e /en) wher e we took the experimental values 



e = 4 and eq = 50 for La^CuO^Je = 4.7 



experimental data of (|Genzel et al 



19891 : 



Genzel et al 



Kress et al 



19891 ) and Eq = 64.7 calculated using the 



19881 ) for YBa2Cu307). The critical points 



ttc = 6.8 for 3D and a,^ = 2.9 for 2D are represented as full dots. (Reprinted with permission after 



Ref. ( Verbist et al. 



1990). ©1990 by Elsevier.) 



of the polaron and bipolaron at rest, respectively. From this condition it follows that the 
Frohlich bipolaron with zero spin is stable (given the effective Coulomb repulsion between 
electrons) if the electron-phonon coupling constant is larger than a certain critical value: 
a > ac- 

A "phase-diagram" for the two continuum polarons — bipolaron system is shown in 
Fig. [14] for 3D and for 2D. The Frohlich coupling constant as high as 6.8 is needed to 
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allow for bipolaron formation in 3D. The confinement of the bipolaron in two dimensions 
facilitates bipolaron formation at smaller a. This can be shown using the PD-type scaling 
relation betw een the free energies F in tw o dimensions F^oip-i (3) and in three dimensions 



F3o(a, [/, /5) (IVerbist et al. 



199Q 



199l|) 



i^2D(«,t/,/5) 



2 ^ .Stt Stt^^ ^, 



(257) 



According to Eq. f l257p . the critical value of the coupling constant for bipolaron formation 



Uc turns out to scale with a factor 37r/4 ^ 2.36 or ac 



(2D) 



a 



(3D) 



/2.36. From Fig. [Mb it is 



seen that bipolarons in 2D can be stable for a > 2.9, a domain of coupling constants which 
is definitely realised in several solids. The "characteristic line" U = 1.526a for the material 
parameters of YBa2Cu307 enters the region of bipolaron stability in 2D at a value of a 
which is appreciably smaller than in the case of La2Cu04. This fact suggests YBa2Cu307 
as a good candidate for the occurrence of stable Frohlich bipolarons. 

An analytical strong-coupling asymptotic expansion in inverse powers of the electron- 
phonon coupling constan t for the large bipolaron energy at T = was derived in 



(ISmondyrev et al. 



19951) 



i?3D(a, u) 



2a' 
"Stt 



A{u) - B{u) + 0{a-'^), 



(258) 



where the coefficients are closed analytical functions of the ratio u = U/a: 



A{u) 



2V2u f 1 + — 
' 128 



3/2 



U 



512 



(259) 



and for B{u) see the above-cited paper. The scaling relation ( 1257P allows one to find the 
bipolaron energy in two dimensions as well. 

The stability of bipolaron s has also beeii exam ined with the use of operator techniques 



with a variational approach (IBassani et al 



19911 ). The bipolaron is bound if the electron- 
phonon coupling constant a is larger than ~ 6 in three dimensions and larger than ~ 2 
in two dimensions, provided the ratio rj = e/eq is smaller than a critical value rjc which 
depends on a. The critical value rjc is larger in the two-dimensional case than in the three- 



dimensional on e. The bipolaron rad ius is shown to be o: 



The results of (IBassani et al. I . Il99ll ) and (IVerbist et al. 



the order o f a few polaron radii. 



1990 



confirm each other. Furthermore, bipolaron states obtained in (IBassani et al. I . Il99ll ) under 



19911) tend to qualit atively 



the assumption that the total linear momentum is conserved, have intrinsically high mobility. 
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In the framework of the renewed interest in bipolaron theory afte r the discovery o f high 



supercond uctivity, an analysis of the optical absorption by large (jPevreese et al. 



Emin 



19931 1 and small (lAlexandrov and Bratkovsky 



1995 



1999bl ) bipolarons was given. 



D. Discrete strong-coupling Frohlich bipolaron 



The Holstein model is an extreme polaron model, and typically yields the highest possible 
values of the (bi)polaron mass in the strong coupling regi me (except the case w hen the 
lattice vibrations are polarised along the hopping direction (ITrugman et al. I . l200ll ) ) . Many 
doped ionic lattices are characterized by poor screening of high-frequency op tical phonons 



and they are more appropriately described by the finite-range Frohlich EPI (lAlexandrov 



19961 ). The unscreened Frohlich EPI provides relatively light lattice polarons fllV.Cp and also 



"superlig ht" small bipo 



in HHM ( lAlexandrov 



aron s , which are several orders of ma ; 



1996 



Alexandrov and Kornilovitch 



gnitude lighter than 



2002 



Hague et al. 



jipolar ons 



2007ah . 



To illustrate the point let us consider a generic "Coulomb- Frohlich" model (CFM) on 



a lattice, w. 



lich explicit 



range EPI (lAlexandrov 



y includes the finite-range Coulomb 



1996 



r epuls ion and the strong long- 



Alexandrov and Kornilovitch 



2OO2I ). The implicitly present 



(infinite) Hubbard U prohibits double occupancy and removes the need to distinguish the 
fermionic spin, if we are interested in the charge rather than spin excitations. Introducing 
spinless fermion annihilation operators Cn and phonon annihi lation operators dm, the Hamil 
tonia n of CFM is written in the real space representation as (lAlexandrov and Kornilovitch 
20021 ) 



if = ^ T(n - n')4cn' + ^ K(n - n')cJjCnc|^/Cn' + 



(260) 



n^n' 



^0 ^(m- n)(e 

r 



where T(n) is the bare hopping integral in a rigid lattice. 

If we are interested in the non- or near-adiabatic limit and the strong EPI, the kinetic 
energy is a perturbation. Then the model can be grossly simplified using the Lang-Firsov 
canonical transformation. In particular lattice structures like a staggered triangular ladder 
in Fig {TBI the intersite lattice bipolarons tunnel already in the first order in t(n). That 
allows us to average the transformed Hamiltonian over phonons to obtain its polaronic part 
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as Hp = Hq + Hpert, where 

Ho = -Ep ^ ciCn + ^ X] ^'(n - n )cI,Cnicn', 



and 



Hpert = ^ t(n - n')4cn'. 
n^n' 

is a perturbation. E'p is the famihar polaron level shift, 

Ep = t^^5(^(m - n)(e ■ e^^^f, 

m 

which is independent of n. The polaron-polaron interaction is 



(261) 



^;(n - n') = K(n - n') - VpH{n - n'), 



where 



Vphin - n) = 2u;o X^5'(m - n)g{iLn - n )(e ■ en,_„)(e ■ e^-n')- 

m 

The transformed hopping integral is t(n — n') = T(n — n') exp[— 5f^(n — n')] with 
g'^{n-n') = ^ ^((m - n)(e.e^_n) x 

m 

[g{ui - n)(e ■ e^-n) - 5'(m - n')(e ■ e^_n')] 
at T ^ cuq. The mass renormalization exponent can be expressed via Ep and Vph as 



(262) 



(263) 



(264) 



g\n-n') = - 



Ep - -Vph{n - n' 



(265) 



When Vph exceeds Vc the full interaction becomes negative and polarons form pairs. The 
real space representation allows us to elaborate more physics behind the lattice sums in Vph 
( lAlexandrov and Kornilovitch 1 . 120021 ). When a carrier (electron or hole) acts on an ion with 
a force f , it displaces the ion by some vector x = f //c. Here k is the ion's force constant. The 
total energy of the carrier-ion pair is — f^/(2A;). This is precisely the summand in Eq. fl26ip 
expressed via dimensionless coupling constants. Now consider two carriers interacting with 
the same ion. The ion displacement is x = (fi + {2)/k and the energy is — / (2A;) — f|/ (2k) — 
(fi ■ f2)/k. Here the last term should be interpreted as an ion-mediated interaction between 
the two carriers. It depends on the scalar product of fi and ¥2 and consequently on the 
relative positions of the carriers with respect to the ion. If the ion is an isotropic harmonic 
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FIG. 15 Mechanism of the polaron-polaron interaction, (a) Together, two polarons (sohd circles) 
deform the lattice more effectively than separately. An effective attraction occurs when the angle 



between xi and X2 is less than 7r/2. (b) A mixed s ituation: ion 1 results in repu. 



two polarons while ion 2 results in attraction. (After (jAlexandrov and Kornilovitch 



sion between 



2002)) 



oscillator, as we assume here, then the following simple rule applies. If the angle between 
fi and f2 is less than 7r/2 the polaron-polaron interaction will be attractive, otherwise it 
will be repulsive. In general, some ions will generate attraction, and some ions - repulsion 
between polarons. 

The overall sign and magnitude of the interaction is given by the lattice sum in Eg. (12631) . 
One should note that according to Eg. (I265p an attractive EPI reduces the polaron mass 
(and conseguently the bipolaron mass), while repulsive EPI enhances the mass. Thus, the 
long-range EPI serves a double purpose. Firstly, it generates an additional inter-polaron 
attraction because the distant ions have a small angle 0. This additional attraction helps to 
overcome the direct Coulomb repulsion between polarons. And secondly, the Frohlich EPI 
makes lattice bipolarons lighter. 

The many-particle ground state of Hq depends on the sign of the polaron-polaron interac 



tion, the carrier density, and the lattice structure. Following (jAlexandrov and Kornilovitch 



2OO2I ). we consider the staggered ladder, FigJT6| assuming that all sites are isotropic two- 
dimensional harmonic oscillators. For simplicity, we also adopt the nearest-neighbour ap- 
proximation for both interactions, g{\) = g, Vc{n) = Vc, and for the hopping integrals, 
T(m) = T(a) for / = n = m = a, and zero otherwise. Hereafter we set the lattice period 
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a = 1. There are four nearest neighbours in the ladder, z = A. The single-particle polaronic 
Hamiltonian takes the form 



Hp = -Ep 5^(4 c„ + pipn) + (266) 

n 

X][^'(cl+lCn + pi+lPn) + t{piCn + pLl^n) + H.C.], 
n 

where Cn and pn are polaron annihilation operators on the lower and upper legs of the 
ladder, respectively, Fig|T6l Using Eqs. fl261|263ll265l) one obtains Ep = 'ig'^ujQ, t' = 
T{a) exp(— T-Ep/Scuo), and t = T{a) ex\i{—?>Ep / Auq) . 

The Fourier transform of Hp yields two overlapping polaron bands, 

Ep{k) = -Ep + 2t' cos{k) ± 2t cos(A:/2) (267) 

with the effective mass m* = 2/|4t' ± t| near their edges. 

Let us now place two polarons on the ladder. The nearest neighbour interaction is 
V = Vc — Ep/2, if two polarons are on different legs of the ladder, and v = Vc — Ep/A, if both 
polarons are on the same leg. The attractive interaction is provided via the displacement of 
the lattice sites, which are the common nearest neighbours to both polarons. There are two 
such nearest neighbours for the intersite bipolaron of type A or B, FigJT6b. but there is only 
one common nearest neighbour for bipolaron C, FigJTBH. When Vc > Ep/2, there are no 
bound states and the multi-polaron system is a one-dimensional Luttinger liquid. However, 
when Vc < Ep/2 and consequently v < 0, the two polarons are bound into an inter-site 
bipolaron of types A or B. 

Remarkably, bipolarons tunnel in the ladder already in the first order with respect to 
the single-polaron tunnelling amplitude. This case is d ifferent from both ons i te bip olarons 



discussed above, and from intersite chain bipolarons of fiBonca and Trugman I . I200ll ). where 
the intersite bipolaron tunnelling appeared in the second order in t as for the on-site bipo- 
larons. Indeed, the lowest-energy configurations A and B are degenerate. They are coupled 
by Hpert- Neglecting all higher-energy configurations, we can project the Hamiltonian onto 
the subspace containing A, B, and empty sites. The result of such a projection is the 
bipolaronic Hamiltonian 



H, = (v - lEp] J2^AiA^ + BiB^] - t' J2[BlAn + B^A^ + H.c. 



(268) 
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(a) 




(b) 



(c) 



(d) 



polaron 'p 




polaron 'c' 




bipolaron B 




bipolaron C 



FIG. 16 One-dimensional zig-zag ladder, (a) Initial ladder with the bare hopping amplitude 
T(a). (b) Two types of polarons with their respective deformations, (c) Two degenerate bipolaron 
configurations A and B (here t,t' are renormalised (i.e. polaronic) hopping integrals), (d) A 
different bipolaron configuration, C, whose energy is higher than that of A and B. 

where An = CnPn and -B„ = p„c,„+i are intersite bipolaron annihilation operators, and the 
bipolaron-bipolaron interaction is omitted. The Fourier transform of Eg. (12681) yields two 
bipolaron bands, 

E2{k) = Vc- ^Ep ± 2t' cos(A;/2). (269) 

with a combined width 4|t'|. The bipolaron binding energy in zero order with respect to t, t' 
is A = 2Ei(0) - ^2(0) = Ep/2 - 

The bipolaron mass near the bottom of the lowest band, m** = 2/t', is 

E„ 



m 



Am' 



1 + 0.25 exp 



(270) 



The numerical coefficient 1/8 in the exponent ensures that m** remains of the order 
of m* even at sufficiently large Ep up to Ep ^ IOcjq- This fact combined with a 



weaker renormalization of m* prov i des a suyerliqht sma ll bipolaron (lAlexandrov 



1996 



Alexandrov and Kornilovitch 



2002 



Hague et al. 



2007af ) 
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E. Discrete all-coupling Frohlich bipolaron 



The CFM model discussed above is analytically solvable in the strong-coupling 
nonadibatic {ujq > T{a)) limit using the Lang-Firsov transforma tion of the Hamil 



tonian, and projecting it on the inter-site pair Hilbert space (lAlexandrov 



1996 



Alexandrov and Kornilovitch 



2002I ). The the ory has been extende d to the w hole parameter 



space 



2007a 



using CTQMC technique for bipolarons (jHague et al. 



2007a 



Refs. (jHague et al. 



b|) simulated the CFM Hamiltonian on a staggered triangular ladder (ID), triangular 



(2D) and strongly anisotropic hexagonal (3D) lattices including triplet pairing. On such 
lattices, bipolarons are found to move with a crab like motion. FigJTBl which is distinct from 
the crawler motion found on cubic lattices flAlexandrov and Ranninger I . Il981al ). Such bipo- 
larons are small but very light for a wide range of electron-phonon couplings and phonon 
frequencies. EPI has been modeled using the force function in the site-representation as in 
Eq. fll94p . Coulomb repulsion has been screened up to the first nearest neighbors, with on 
site repulsion U and nearest-neighbor repulsion Vc- The dimensionless electron-phonon cou- 
pling constant A is defined as A = f^{0)/2MLj'^zT{a) which is the ratio of the polaron 
binding energy to the kinetic energy of the free electron zT{a), and the lattice constant is 
taken as a = 1. 

Extending the CTQMC algorithm to systems of two particles with strong EPI and 
Coulomb repulsion solved the bipolaron problem on different lattices from weak to strong 
coupling in a realistic parameter range where usual strong and weak-coupling limiting ap- 
proximations fail. Figure [T7] shows the ratio of the polaron to bipolaron masses on the 
staggered ladder as a function of effective coupling and phonon frequency for Vc = 0. The 
bipolaron to polaron mass ratio is about 2 in the weak coupling regime (A ^ 1) as it should 
be for a large bipolaron, section I VI. CI In the strong-coupling, large phonon frequency limit 
the mass ratio approaches 4, in agreement with strong-coupling results, section IVI.DI In a 
wide region of parameter space, the bipolaron/polaron mass ratio has been found between 
2 and 4 and a bipolaron radius similar to the lattice spacing. Figs. [171 Thus the bipolaron 
is small and light at the same time. Taking into account additional intersite Coulomb re- 
pulsion Vc does not change this conclusion. The bipolaron is stable for Vc < 4T(a). As 
Vc increases the bipolaron mass decreases but the radius remains small, at about 2 lattice 
spacings. Importantly, the absolute value of the small bipolaron mass is only about 4 times 
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FIG. 17 Polaron to bipolaron mass ratio for a range of w = ujQ/T{a) and A on the staggered ladder. 
Mobile small bipolarons are seen even in the adiabatic regime u) = 0.5 for couplings A up to 2.5 
(Reproduced from J. P. Hague et al., Phys. Rev. Lett 98, 037002 (2007), (c) American Physical 
Society, 2007). 

of the bare electron mass mo, for X = uj/T{a) = 1 (see Fig. [T7I) . 

Simulations of the bipolaron on an infinite triangular lattice including exchanges and 
large on-site Hubbard repulsion U = 20T(a) also lead to the bipolaron mass of about 
QrriQxy and the bipolaron radius Rbp ~ 2a for a moderate coupling A = 0.5 and a large 
phonon frequency u = T{a) (for the triangular lattice, rriQxy = l/3a^T(a)). Finally, the 
bipolaron in a hexagonal lattice with out-of-plane hopping T' = T(a)/3 has also a light 
in-plane inverse mass, m** ^ 4.5moxy but a small size, Rbp ~ 2.6a for experimentally 
achievable values of the phonon frequency u = T{a) = 200meV and EPI, A = 0.36. Out- 
of-plane m** ^ 70moz is Holstein like, where = l/2cPT', [d is the inter-plane spacing). 
When bipolarons are small and pairs do not overlap, the pairs can form a BEG at Tbec = 

If we choose realistic values for the lattice constants of 0.4 
nm in the plane and 0.8 nm out of the plane, and allow the density of bosons to be nB=0.12 
per lattice site, which easily avoids overlap of pairs, then Tbec ~ 300K. 
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F. Polaronic exciton 



Finally let us briefly mention works on electron-hole bound states coupled with 



phoii o ns. S uch e xcitonic polaron states have been analysed by (jiadonisi and Bassani 



1983 



1984 



19871 ). The binding energies of excitonic states in interaction with LO 



phonons were computed usin g a phonon coherent state and applying a variational method 



(jiadonisi and Bassani 



19831 ). Incomplete relaxation of the lattice is found. The binding en- 



ergies are larger than those obtained with static dielectric screening when the polaron radius 
is comparable to the exciton radius. Analysing the re lative intensities of the one-ph onon and 
zero-phonon lines for a number of semiconductors, (jiadonisi and Bassani . 1983) observed 
that the zero-phonon exciton states are generally much more probable than the phonon 
replicas. 

The binding energy of a core exciton, which depends on the interaction of the conduction 
electron and of the core hole with the valence electrons via a Frohlich-type coupling with the 
electron-hole pairs, i n addition to the Coulomb at traction, was calculated within a functional 
variational method (jiadonisi and Bassani I . Il984f ). When the exciton radius is comparable 
to the polaron radius, the static dielectric screening reduces and the core exciton binding 
energy increases. The excitonic-polaron effective mass renormalization was analysed by 



(jiadonisi and Bassani 



19841 ) using a variational numerical approach. Evidence for this 



renormalization is found from the experimental data on polariton dispersion in hexagonal 
CdS and in CuCl. 

?he angular moment um as a constant of the motion of a Frohlich polaron was introduced 



by (jEvrard et al. 1 . 119701 ). R epresenting the exciton ic polaron problem in angular coordinates 



and solving it numerically, (jiadonisi et al. 



19891 ) showed that exciton states with various 



total angular momenta L are differently affected by the interaction with LO phonons. How- 
ever, when extending this approach to the bipolaron problem, the same authors found that 
the Frohlich interaction is not strong enough to guarantee a stable bipolaron state at least 
in the case of cubic materials. 

Very recently th e polaronic exciton pro blem has been solved using the approximation-free 
DQMC technique (jBurovski et al. 1 . 12008| ). Numerically exact results for the wave function, 
ground-state energy, binding energy and effective mass of this quasiparticle were calculated, 
and the frequently used instantaneous approximation to the retarded interaction due to the 
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exchange of phonons was critically analysed. 



VII. CURRENT STATUS OF POLARONS AND OPEN PROBLEMS 



At present, the basic properties of single polarons are well understood theoretically, and 
to a large extent they are analytically under control at all coupling. It is remarkable how 
the Frohlich continuum polaron, one of the simplest examples of a Quantum Field Theoret- 
ical problem, as it basically consists of a single fermion interacting with a scalar Bose field, 
has resisted full analytical solution at all coupling since ~ 1950, when its Hamiltonian was 
first written. Although a mechani sm for the optical absorption of Frohlich polaro ns was 



already proposed a long time ago (IDevreese et al. 



1972a 



subtle characteristics were onl y clarified very recently ( 



( ICataudella et al 



2007 



i ng numerical DQMC studies (jMishchenko et al. 1. 120031) and improved analytical methods 



Devreese and Klimin 



Kartheuser et al. 



De Filippis et al. 



1969I ). some 



20061 ) by combin- 



20061 ) (IIII.FI) . Of special interest are several 



sum rules derived for the optical conductivity spectra of arbitrary-coupling Frohlich polarons 



( Devreese et al. 



1977 



Lemmens et al. 



19731 ) . A variety of magneto-optical and transport' 



experiments were succes s fully analysed with F r ohlich polaron theory (see e . g. (jPevreese 



2004 



Hodby et al. 



1987: iPeeters and Devreese 



1984 



Peeters and Devreese 



J 



19861) and ref- 



erences therein). 

The charge carriers in a rich variety of systems of reduced dimension and dimensionality 
(submicron- and nanostructures including heterojunctions, quantum wells, quantum wires, 
quant um dots etc.) turn out to b e Frohlich polarons. Several scaling relations were de- 
rived (IPeeters and Devreesd . 119871 ). which connect polaron characteristics (the self-energy, 
the effective mass, the impedance and the mobility) in different dimensions. Quite generally 
confinement enhances EPI and the tendency to polaron formation. A new aspect of the po- 
laron concept has been investigated for semiconductor structures at nanoscale: the exciton- 
phonon states are no t factorisable i nto a n adiabatic product Ansatz, so that a non-adiabatic 



treatment is needed (IFomin et al. 



19981 ). The excito nic polarons in nanost ructures lead to 



the existence of phonon replicas in the luminescence (j Verzelen et al. 1 . 12002| ). Experimental 
evidence of the enhanced phonon-assisted absorption due to effects of non-adiabaticity has 



been provided by the multi-phonon photoluminescence 



PL) spectra observed under s elective 



excitation in self-assembled InAs/GaAs quantum dots (IGarcia-Cristobal et a/.l . Il999l ) and by 
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the photoluminescence excitation (PLE) m easurements on sing le self-assembled InAs/GaAs 
(ILemaitre et a/.l . l200ll ) and InGaAs/GaAs fjZrenner et a/.l . l200ll ) quantum dots. The polaron 



concept was also i nvoked for t 



le exp lanation of the P LE measurements o n self-organised 



lnx.Gai_a;.As/GaAs flHeitz et al 



200 If ) and CdSe/ZnSe flWoggon et a/.l . 120031 ) quantum dots. 
The Frohlich polaron has led to many generali sations. The s t abilit y region of the 



Frohlich la rge bipolaron is now firmly established 



1990 



199ll ) fIVI.Cp . Here the surprise is double (cf. (jPevreese et al. 



Bassani et al. 



19911: 



Verbist et al. 



1995 



Verbist et al. 



199ll )): a) only in a very limited sector of the phase diagram (Coulomb repulsion versus a) 
the bipolaron is stable, b) most traditional Frohlich polaron materials (alkali halides and 
the like) lie completely outside (and "far" from) this bipolaron stability sector, but sev- 
eral cuprate superconductors lie very close and even inside this very restricted area of the 
stability diagram. The stability of a strong-coupling singlet bipolaron was studied in two- 
and three-dimensional parabo lic quantum dots using the Landau-Pekar variational method 



(IMukhopadhyay et al. 



19961 ). It was shown that the confining potential of the quantum dot 



reduces the stability of the bipolaron. A theory of bipolaron states in a spherical parabolic 
potential well was further developed applying the Feynman variational principle. The basic 
parameters of the bipolaron ground state (the binding energy, the number of phonons in the 
bipolaron clou d, and the bipo l aron radius were studied as a function of the radius of the 
potential well (jPokatilov et a/.l . Il999l ). It was found that confinement can enhance the bipo- 
laron binding energy, when the radius of a quantum dot is of the same order of magnitude 
as the polaron radius. A unified insight into the stabil ity criterion for bipo l aron formation 
in low-dimensionally confined media was provided by (jSenger and Ercelebi l2000l ) using an 
adiabatic variational method for a pair of electrons immersed in a reservoir of bulk LO 
phonons and confined within an anisotropic parabolic potential box. Bipolaron formation in 
a two-dimensional lattice with harmonic confinem ent, representing a simplifi ed model for a 
quantum dot, was investigated by means of QMC (IHohenadler et al. U2007al ). This method 
treats all interactions exactly and takes into account quantum lattice fiuctuations. Calcula- 
tions of the bipolaron binding energy reveal that confinement opposes bipolaron formation 
for weak electron-phonon coupling but abets a bound state at intermediate to strong cou- 
pling. We also mention the exciton-polaron formation in nanostructures and quantum-light 
sources studied r ecently by QMC in lattice models with short- or long-range carrier-phonon 
interaction (see f Hohenadler et al. . boOTb ). and references therein) 
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The richness and profundity of Landau-Pekar's polaron concept is further illustrated by 
its extensions to discrete (lattice) polarons. Even the simplest two-site polaron model by 
Holstein (lIV.Ap proved to be very useful for a qualitative understanding of nontrivial fea- 
tures of the polaron proble m, and f or ob t aining some novel analytical and semi-a n alytic al 



results (see, for example, (IBerciu 



2007 



Han et al. 



2002 



Kudinov and Firsov 



19971)). 



The "1/A" expansion technique based on the Lang-Firsov transformation (1IV.B.1|) and un- 
biased numerical analysis of the finite and infinite Holstein and Frohlich models combining 
Lanczos diagonalisations of clusters, density matrix renormalisation group, cluster pertur- 
bation theory techniques, DMFT and different QMC algorithms, allowed for a description 
of properties of a single lattice polaron a nd a lattice bipolaron. The Lang-Firsov canonical 



transformation (ILang and Firsov 



19621 ) was proven particularly instrumental in calcula- 



tion of different kinetic and optical coefficients, which can be represented as expansions in 
powers of the unrenormalised hopping integral t (jV]). Sometimes it is possible to sum the 
expansion and get results, which are valid for arbitrary values of parameters providing the 



underst a nding o 



hopping (IFirsov 



Recent ED (Fehske and 



the c rossover region from the Boltzmann kinetics to thermally activated 



2007h. 



Trugman I . 120071), CTQMC flKornilovitch I . 120071 ) and DQMC 



(IMishchenko and Nagaosal . 120071 ) techniques allow for determination of the ground-state and 



excited states of lattice polarons with arbitrary precision in the thermodynamic limit for any 
dimension and any type of lattices (1IV.E.2I1 . The spectral properties (e.g. photoemission) , 
optical response and thermal transport, as well as the dynamics of polaron formation in the 
Holstein model have been numerically analysed for all EPI strengths and phonon f reque n- 
cies, including the intermediate-coupling regime (see IVl and (IFehske and Trugman 1 . 120071 )). 
CTQMC methods have proven to be powerful and versatile tools providing unbiased results 
for the polaron properties in any lattices for any-range EPI, including Jahn- Teller polarons 
(lIV.Ep . Combining the Lang-Firsov transformation and quantum Monte Carlo simulations 
allows for an exact sampling without autocorrelations, w hich proves to be an enormous ad 
vantag e for small phonon frequencies or low temperatures (IHohenadler and von der Linden 



2007ch . 



Importantly, variational and numerical techniques confirmed that the Frohlich and 
Holstein-Lang-Firsov theories are asymptotically exact in the weak, A ^ 1, and strong- 
coupling, A ^ 1, regimes, respectively, and the polaron formation represents a contin- 
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uous crossover of the ground state (lIV.Ep . The crossover is related to the exponential 
increase of the effective mass, and the band narrowing with a strongly suppressed elec- 
tronic quasiparticle residue and the Drude weight, accompanied by an increase of the 
weight (IV. Cp ■ These fea t ures s trongly depend on the phonon dis- 



persion ( 


Zoh 


Fehske et al. 



2002 



19981), EPI radius (Alexandrov 



Perroni et al. 



more pronounced in higher dimensions (IFehske et al 



2005!) ■ lattice geometry 



1996 



Alexandrov and Kornilovitch 



2002 



(IHague et al. 



Hague et al. 



1999; 



2006al) and are 



2006ah . Remark- 



ably, the unscreened Frohlich EPI provides relatively light lattice pola rons (IIV.CII. which 



are several orders of magnitud e light er than the Holstein small polarons (lAlexandrov 



1996 



Alexandrov and Kornilovitch 



I999I ) at strong coupling. This classification of weak- and 



strong-coupling regimes still leaves room for continuum Frohlich polarons of not only weak 
but also intermediate and - in the theoretical analysis - of strong coupling classified with 
respect to the Frohlich electron-phonon coupling constant a (]II.D|) . 

While the single polaron has been actively researched for a long time and is now well un- 
derstood, the multi-polaron physics has gained particular attention in the last two decades. 
It has been found — unexpectedly for many researchers — that the Migdal-Eliashberg theory 
breaks down already at A ~ 1 for any adiabatic ratio ujq/ Ep. The effective parameter Xujq/ Ep 
becomes large at A > 1 since the bandwidth is narrowe d and the Fermi energy, Ep is renor- 



malised down exponentially (lAlexandrov 



1983 . 



2001al ). Extending the BCS theory towards 



the strong interaction between electrons and i on vibrations, a charged Bose gas o f tightly 



bound on-site small bipolarons was predicted (lAlexandrov and Ranninger 



1981a), with a 



further conclusion that the highest superconducting transition temperature is attained in 
the crossover region of EPI stren gth between the BCS-like polaronic and bipolaronic super- 



conductivity (lAlexandrov 



also two-site bipolarons (^Aubry 



Macridin et al. 



1983 ). Subsequent studies of the Holstein-Hubbard model 



1995 



Bonca et al. 



200( 



La Magna and Pucci 



bund 
I997I : 



20041 ). Taking into account that many advanced materials with low den- 



sity of free carriers and poor mobility are characterized by poor screening of h igh-frequency 



optical phonons, the Coulomb-Froh 



1996 



Alexandrov and Kornilovitch 



i ch lat tice polaron model was introduced (lAlexandrov 



2002) • The large Hubbard U and intersite Coulomb re- 



pulsions and the unscreened Frohlich EPI provide "superlight" but small intersite bipolarons 



dVLDl) . More r ecent CTQMC simu 



ations of intersite small bipolarons in the Coulomb- 



Frohlich model (IHague et al. I . l2007al ) have found such quasiparticles in a wide parameter 
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range with achievable phonon frequencies and couphngs. They could have a superconducting 
transition in excess of room temperature. 

The many-body theory for polarons has been developed for extreme weak and strong 
coupling regimes. It became clear how — in the weak-coupling limit — this prob- 
l em can be reduced to the study of the structure factor of a uniform electron gas 



(ILemmens et al. 



19771). For strong coupling the prob lem is reduce d to an interacting 



Bose gas of on-site (lAlexandrov and Ranninger 



Alexandrov and Kornilovitch 



2002 



Hague et al. 



1981a[) or intersite (lAlexandrov 



1996 



2007al ) small bipolarons in the dilute 



system. 

While correlation effects in transport through metallic quantum dots with repulsive 
electron-electron interactions received considerable attention in the past, and continue to 
be the focus of intense investigations, much less has been known about a role of attractive 
correlations between small polarons mediated by EPI in molecular qu antum dots (MQD) 



and nanowires. In the framework of the negative U Hubbard model (lAlexandrov et al. 



2003! ) it has been found that the attractive electron correlations within the molecule could 



lead to a molecular switching effect where I-V characteristics have two branches with high 
and low current at the same bias voltage. The switching phenomenon has been also pre- 
dicted by a theory of correlated polaron transport with a full account of both the Coulomb 



repu 



2003 



sion and 



Ermakov 



PI in MQD weakly coupled with electrodes (lAlexandrov and Bratkovsky 



20001) . Ref . (lAlexandrov and Bratkovsky I . l2003l ) has shown that while the 



phonon side-bands significantly modify the shape of hysteretic I-V curves in comparison with 
the negative-t/ Hubbard model, switching remains robust. It shows up at sufficiently low 
temperatures when the effective interaction of polarons in MQD is attractive and the molec- 



ular level is multiply degenerate . Importantly, the switching 



and two-fold degenerate MQDs (lAlexandrov and Bratkovsky 



las not been found in non- 



20031). When the polaronic 



energy shift is very large, the effective charging energy of m olecules can 



favoring ground states with even numbers of electrons. Ref.( lKoch et al 



jecom e negative, 



2006! ) has shown 



that charge transport through such molecules is dominated by tunneling of bipolarons which 
coexists with single-electron cotunneling. 

The conductance of deformable molecules with a local magnetic moment has been studied 
in the framework of a two-impurity Anderson model with positive and negative electron- 
electron interactions and in the two-impurity Anderson-Holstein model with a single phonon 
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mode (IZitko and Bonca 



20061 ). It has been shown that the spin and charge Kondo effects 



can occur simultaneously at any coupling strength. At finite bandwidth and strong cou- 
pling the lattice effects lead to a renormalization of the effective Kondo exchange constants; 
nevertheless, universal spin and charge Kondo effects still occur. 

The theory of dense polaronic systems in the intermediate regime remains highly cumber- 
some, in particular when EPI competes with strong electron correlations. The corresponding 
microscopic models contain (extended) Hubbard, Heisenberg or double-exchange terms, and 
maybe also a coupling to orbital degrees of freedom along with strong EPI, so that even 
numerical solutions with the same precision as in the dilute (bi)polaron case are often prob- 
lematic. A number of ED, QMC, DMFT and combined numerical results give strong evi- 
dence that the tendency towards lattice polaron formation is enhanced in strongly correlated 



recent reviews by ( 


Edwards . 


2002 


Fehske and Trueman 


Hohenadler and von der Linden 




>007c 


) for more details) 



2007 



Hohenadler et al. 



2005b 



relations enhance EPI, resulting in polaron fo rmation for moderate 



also EPI strongly enhances spin correlations (ICapone and Ciuchi 



coup l ing strength, but 



2002 



Macridin et al. 



20061 ). Some of these studies show that increasing carrier density could be accompanied 



by a dissociation of po 



arons leading to normal m etallic behavior in the intermediate cou- 



pling adiabatic regime (IHohenadler et al. I . l2005bl ) that is reminiscent of the "overcrowding" 



effect hypothesized by Mott (1995). On the other hand for parameters favoring small po- 
l arons no such density -driven crossover occurs in agreement with simple analytical arguments 



( lAlexandrov 



2001al ). Thorough investigations of these models will definitely be a great 



challenge in the near future. Here we have focused on the single and two-body problems 
leaving theories of strongly-correlated polarons and t heir applications to high-tempera ture 
superconductors and OMR oxides for future reviews (lAlexandrov and Devreese 1 . 120091 ). 

Finally, exactly solvable models might give a rather limited, sometimes misleading, de- 
scription of polarons in real systems. Qualitative inconsistencies can arise when coupling 
is assumed to be just to one phonon mode, often taken as dispersionless, and ad- hoc ap- 
proximations for EPI matrix elements are applied. Moreover electronic nanoscale disorder 



and long-r ange strain fie 



transport fjPhillips et al. 



ds can interweave with the microscopic mechanisms of polaronic 



20031 ). Hence ab-initio calculations of the phonon spectrum, EPI 



and polaron properties beyond the adiabatic Born-Oppenheimer approximation are required 
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in manv cases for which the 



2008 



;heory and experiment can be compared in detail (IBanacky I . 



Shluger and Stoneham 



19931). 
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